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AEROELASTIC ANALYSIS FOR ROTORGRAFT 


IN FLIGHT OR IN A WIND TUNNEL 
Wayne Johnson 

Ames Research Center, NASA 
and 

Ames Directorate, USAAMRDL 


SUMMARY 


An analytical model is developed for the aeroelastic behavior of a rotor- 
craft in flight or in a wind tunnel. A unified development is presented for a 
wide class of rotors, helicopters, and operating conditions. The equations of 
motion for the rotor are derived using an integral Newtonian method, which 
gives considerable physical insight into the blade inertial and aerodynamic 
forces. The rotor model includes coupled flap-lag bending and blade torsion 
degrees of freedom, and is applicable to articulated, hingeless, gimballed, 
and teetering rotors with an arbitrary number of blades. The aerodynamic 
model is valid for both high and low inflow, and for axial and nonaxial flight. 
The rotor rotational speed dynamics, including engine inertia and damping, and 
the perturbation inflow dynamics are included. For a rotor on a wind-tunnel 
support, a normal mode representation of the test module, strut, and balance 
system is used. The aeroelastic analysis for the rotor craft in flight is 
applicable to a general two-rotor aircraft, including single main-rotor and 
tandem helicopter configurations, and side-by-side or tilting proprotor air- 
craft configurations. An arbitrary unaccelerating flight state is considered, 
with the aircraft motion represented by the six rigid body degrees of freedom 
and the elastic free vibration modes of the airframe. The rotor model 
includes rotor-rotor aerodynamic interference and ground effect. The aircraft 
model includes rotor-fuselage-tail aerodynamic interference, a transmission 
and engine dynamics model, and the pilot's controls. A constant-coefficient 
approximation for nonaxial flow and a quasistatic approximation for the low- 
frequency dynamics are also described. The coupled rotorcraft or rotor and 
support dynamics are described by a set of linear differential equations, from 
which the stability and aeroelastic response may be determined. 


1.0 INTRODUCTION 


The testing of rotorcraft in flight or in a wind tunnel requires a con- 
sideration of the coupled aeroelastic stability of the rotor and airframe, or 
the rotor and support system. Even when the primary purpose of the test is to 
measure the rotor performance, experience shows that the question of dynamic 
stability may be ignored only at the risk of catastrophic failure of the air- 
craft. Moreover, in the development of advanced rotor systems, the measure- 
ment and verification of the dynamic stability are themselves major goals of 
the test. Thus it is most desirable to have an analytical model of the 




i£6MiM %c/tdf f 'aetd£lagticlij or ^effco^tVr f iXght f d^a&i£s .> ' • Fur^ermore , an 
^iiftalytlcai mod^l idt th£ : to tor t , ±S ''?4(jt(i'r4cl as vttie Basis for more ex fen- ' 
^‘f^e • inVesfigatidiisf s df the ^aefo^idLstffc "behavior , ; s : Uc : lf " a^‘ autom^iQ' 1 confr df''.*'’ 
system design. *’ Tft e p r ihdip'a 1 " 1 Impair 6rf of" tlie aria" lys ^s’ 'a va liable '' in "t he 
literature is that they are not applicable to a wide class of rotorcraft. 
Typically, aeroelastic stability analyses have been developed in response to e 
concern with some specific dynamic problem, and thus are suitable only for a 
particular type of rotor or a limited range of operating conditions. Often 
the model does not include the entire aircraft or does not consider the rotor 
shaft motion at all. This report presents the unified development of an 
aeroelastic analysis for a wide class of rotors and rotorcraft. A thorough 
documentation of the analytical model is required to interpret the results of 
past and future investigations of rotorcraft dynamic behavior using this 
model. 


The usefulness of an analysis depends on its ability to handle a large 
class of problems; therefore, the scope of the aeroelastic model developed 
here is kept as wide as possible. The rotor model is applicable to articu- 
lated, hingeless, gimballed, and teetering rotors with an arbitrary number of 
blades (including two-bladed rotors). This generality is accomplished by 
using a modal representation for the blade coupled flap and lag motion, with a 
gimbal or teeter hinge included in the hub from the beginning of the analysis. 
Then an articulated or hingeless rotor may be modeled by dropping the gimbal 
degrees of freedom and using the modes of a hinged or cantilever blade, 
respectively. For a gimballed (or teetering) rotor, the gimbal degrees of 
freedom are retained, with cantilever modes for the blade bending motion. 

The description of the blade motion includes rigid pitch deflection due to 
control-system flexibility and elastic torsion modes. The rotor model also 
includes the rotational speed dynamics (with the effects of engine inertia and 
damping) and perturbation inflow dynamics to account for the unsteady aero- 
dynamics of the rotor. 

The aeroelastic analysis of the rotorcraft in flight is applicable to a 
general two-rotor aircraft, including single main-rotor and tandem helicopter 
configurations and side-by-side or tilting proprotor aircraft configurations. 
An arbitrary, unaccelerated equilibrium flight state is considered, with the 
aircraft motion represented by the six rigid body degrees of freedom and the 
elastic free vibration modes of the airframe. The rotor model for the air- 
craft in flight includes rotor-rotor aerodynamic interference and ground 
effect. The aircraft model includes rotor-fuselage-tail aerodynamic inter- 
ference, a transmission and engine dynamics model, and the pilot’s controls. 

In part I, the rotor model is derived and also the model for the coupled 
rotor and wind-tunnel support dynamics. The equations of motion for the rotor 
are developed using an integral Newtonian method rather than the more common 
Lagrangian or differential Newtonian methods. The integral Newtonian approach 
allows greater use of engineering experience in deriving the equations and 
provides considerable physical insight into the inertial and aerodynamic 
forces of the rotor blade. By introducing a vector representation of the 
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coupled flap /lag bending displacement, a very compact form Is obtained for 
the blade bending equations of motion. In part II, the aeroelastic analysis 
for the rotor craft In flight is derived. The coupled rotorcraft or rotor and 
support dynamics are described by a set of linear differential equations, from 
which the stability and aeroelastic response may be determined. 
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PART I. AEROELASTIC ANALYSIS FOR A ROTOR IN A WIND TUNNEL 


The development of the aeroelastic analysis for a helicopter rotor and a 
wind-tunnel support (fig. 1) begins with a consideration of the rotor model in 
section 2. The structural, inertial, and aerodynamic forces on the blade are 
derived, followed by a consideration of the engine dynamics and the rotor 
inflow model. Then the equations of motion for the rotor are presented for a 
three-bladed rotor. Section 3 discusses further some details of the rotor 
model; section 4 extends the analysis to an arbitrary number of blades. In 
section 5, the support equations of motion are presented. Finally, in sec- 
tion 6, the rotor and support equations are combined to construct the equa- 
tions of motion for the coupled system. Note that, although the analysis 
begins with dimensional quantities, in the final equations all parameters are 
dimensionless, based on air density p, rotor rotational speed Q, and rotor 
radius R. 


2. ROTOR MODEL 


This section develops the aeroelastic analysis of the helicopter rotor. 

The rotor motion is represented by the following degrees of freedom: coupled 

flap and lag bending modes, rigid pitch motion (due to control-system flexi- 
bility), blade elastic torsion modes, rotor rotational speed perturbation, and 
gimball or teetering hinge motion (when required) . The six components of the 
rotor shaft linear and angular motion are included, as well as the rotor blade 
pitch control. Three components of aerodynamic gusts are included as external 
disturbances. The rotor hub and root representation includes: precone, droop, 

and sweep; pitch bearing radial offset; feathering axis droop and sweep; and 
gimbal undersling and torque offset. Chordwise offsets of the blade center of 
gravity, aerodynamic center, and tension center are included in the blade 
representation. The undeformed elastic axis of the blade is assumed to be a 
straight line. The rotor aerodynamic model is generally valid for high and 
low inflow and for axial and nonaxial flight. The effects of reverse flow, 
compressibility, and static stall are included. 

The linear differential equations describing the motion of the three- 
bladed rotor are presented in matrix form, together with equations for the 
forces and moments acting on the rotor hub. Two cases are considered: axial 

flow, which is a constant coefficient system, and nonaxial flow, which is a 
periodic coefficient system. Also, in section 2.7, a constant coefficient 
approximation for the nonaxial flow equations, using the mean values of the 
coefficients in the nonrotating frame, is derived. The development of the 
rotor model begins with the analysis of the blade structural moments . 


2 . 1 Structural Analysis 

The structural analysis consists of an engineering beam theory model for 
the coupled flap/lag bending and torsion of a rotor blade with large pitch and 
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twist. A high aspect ratio (of the structural elements) is assumed, so the 
beam model is applicable. The objective is to relate* the bending moments at 
the section, and the torsion moment, to the blade deflection and elastic tor-/ 
sion at that section. The analysis follows the work of references 1 to 3. 

2.1.1 Geometry- The basic assumptions are that an elastic axis exists, 
and the undeformed elastic axis is a straight line; and that the blade has a 
high aspect ratio (of the structural elements) , so engineering beam theory 
applies. Figure 2 shows the geometry of the undeformed blade. The span 
variable r is measured from the center of rotation along the straight elas- 
tic axis. The section coordinates x and z are the principal axes of the 
section, with the origin at the elastic axis. Then, by definition, 
Jsection( xz )dA = 0. Really, the integral is over the tension-carrying ele- 
ments, that is, a modulus weighted integral: / xzE dA = 0. This remark holds 

for all section integrals in the structural analysis. The tension center 
(modulus weighted centroid) is on the x axis, at a distance Xq aft of the 
elastic axis: /x dA = xqA and J*z dA = 0. Again, these are modulus weighted 
integrals. If E is uniform over the section, then x^ is the area cen- 
troid. If the section mass distribution is the same as the E distribution, 
then the tension center coincides with the section center of gravity. 


The angle of the major principal axis (x axis) with respect to the hub 
plane is 0. The existence of the elastic axis means that elastic twist about 
the elastic axis occurs without bending. Generally, the elastic torsion 
deflection will be included in 0. The blade pitch bearing is at the radial 
station rp^. The blade pitch is described by root pitch 0° (rigid pitch 
about the feathering axis, including that due to the elastic distortion of the 
control system), built-in twist 0 tw , an< * elastic torsion about the elastic 
axis 0 e . So 0 = 0° + 0 tw + 0 e , where 0°(if>) is the root pitch, 0(rp^) = 0°; 
0 tw ( r ) is the built-in twist, 0 tw (rpA) = 0; and Q e (r,\p) is the elastic tor- 
sion, 0 e (rp A ,i|>) = 0. There is shear stress in the blade due to 0 e only. It 
is assumed that 0 e is small, but 0° and 0 tw are allowed to be large. 

_ The unit vectors in the rotating hub plane axis system are ip, jg, and 
kg (fig. ^2). T^e unit vectors for the principal axes of the section (x,r,z) 
are i, j, and k; these vectors are for no bending, but include the elastic 
torsion in the pitch angle 0. So the principal unit vectors are rotated by 
0 from the hub plane: 



2.1.2 Description of bending- Now the engineering beam theory assumption 
is introduced: plane sections perpendicular to the elastic axis remain 

so after the blade bends. Figure 3 shows the geometry of the deformed sec- 
tion. The deformation of the blade is described by (a) deflection of the 
elastic axis, x Q , r Q , and z D ; (b) rotation of the section due to bending, by 
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(f> x and <f> z ; and (c) twist about the elastic axis, 0 e , which is implicit in i 
and k. The quantities x Q , r Q , z Q , <|) x , <J> Z , and 0 e are assumed to be small. 

The unit vectors of the unbent cross section are i,^j, and k. The unit 
vectors of the deformed cross section are l^g, iXS> an< * ^XS» where i^g anc ^ 
k^g are the principal axes of the section and j^g is tangent to the 
deformed elastic axis. It follows then that 

-y -y -y 

i™ = ± + + z j . 


XS 

y 

’xs 


->* ->■ 
jvo = j - <t>A + O 


X 


XS 


-> 

= k - <j) j 
Y x J 


Now, by definition, jf xs = dr/ds, where r = x Q i + (r + r Q )y + z Q k and s is 
the arclength along the deformed elastic axis. Hence, to first order, 


= j + (x Q i + z 0 £) ’ 


J XS 


- j + (x^ + z o 0')i + (z' Q - x o 0')k 
It follows that the rotation of the section is 

“♦z = x o + z o 0 ' 

^x = z o " x o 0 ' 


or 


♦x* + ‘t’z^ = " x o^> ' 


The undeflected position of the blade element is r = r + xi + z£, and 


the deflection position is 


r = (r + r + x i + z £ + xi. 


XS 


+ z 




= r J + + r 0 j" + z 0 £ + (x<j) z - z* x )i + xi + z£ 


The first term in the deflected position is the radial station, the next three 
terms are the deflection of the elastic axis, the next term is the rotation of 
the section, and the final two terms are the location of the point on the 
cross section. For now, the elastic extension r Q is neglected. The strain 
analysis is simplified since then, to first order, s = r; r Q just gives a 
uniform strain over the section, which may be reintroduced later. 

2.1.3 Analysis of strain- The fundamental metric tensor g^ of the 
undistorted blade is defined by 
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/ 


(ds) 2 = dr • dr 


■ ft ft ■ ft ft 


= Sim dx m dx n 

where ds is the differential length in the material and ^ are general 
curvilinear coordinates. Similarly, the metric tensor of the deformed 

blade is 

(dS) 2 = dR • dR 
9R 


' ft ft ’ ft ft 


~ dx dx 
mn m n 


Then the strain tensor is defined by the differential length increment: 

2 ^mn dx m dx n = (dS > 2 " < ds > 2 


or 


Y mn 2 ^tnn ^mn^ 


For engineering beam theory, only the axial components of the strain and 
stress are required. (For a full exposition of the analysis of strain, see 
ref. 2.) 

The metric of the undeformed blade (no bending and no torsion, so^ 

0 ? = 9£ w ) °btained from the undistorted position vector r = xi 4 - r;f + z£, 
giving 

g = = 1 + 0’^(x 2 4- z 2 ) 

rr 3r 3r tw 

The metric of the deformed blade, including bending and torsion, is similarly 
obtained from the position vector r = (x + x Q )i 4* (r 4- x<j> z - z<j> x )j 4* (z4z Q )k, 
giving 


3r 9r 


-rr 3 r 3r = (1 + X(f) z " z 4>i> + K + 0?(z + z o> 1 + t z o " 0 ? < x + x o> ^ 

Then the axial component of the strain tensor is 
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Y rr 2 ^ G rr 

= \ | Cl + x(J)^, - z<pp 2 ~ 1 + [x^ + 9' (z + z Q ) ] 2 

- e tw z2 + f z o " 01 (x + x 0 )] 2 - 0^x 2 j 

The linear strain (for small x Q , z Q , d e , <p x , and <J> Z ) is 

Y rr ~ e rr = ~ Z K + 0 tS (xx o + zz o ) + 0 tw [zx o " xz o + 9 e (x2 + z2)] 


The strain due to the blade tension, e^, is a constant such that the ten- 
sion is given by the integral over the blade section: 

T = / Ee rr dA = e T / E dA 

Substituting for e rr and using the results Jz dA = 0, Jx dA = x c A, and 
J*(x 2 + z 2 )dA = Ip = kp 2 A (where kp is the modulus weighted radius of gyra- 
tion about the elastic axis) gives 


e T 


EA 


= <f>z x C + 


0 tw x o x C “ 6 '- z « x ^ + 0,0 ~k 2 + 


tw o C 


enp 


In this expression, the strain due to the blade extension r G has been 
included. It follows that the strain may be written: 


e rr = e T + (x " X C^K - Ox) " Z ^L + 0 tw+ Z ) + 0 tw 0 4< x2 


+ z 2 


- kp) 


2.1.4 Seat-ion moments- To find the moments on the section, the second 
engineering beam theory assumption is introduced: all stresses except 

o rr are negligible. The axial stress is given by a rr = Ee rr . The direction 
of a rr is 

g = 8r/8r 

|3r/3r| 

The moment on the deformed cross section (fig. 4) is $ = M^g + M,. j xs + M z lc xs . 
The moment about the elastic axis due to the elemental force a rr dA on the 
cross section is 

dM - (xi xs + zk xs ) x (a rr a)dA 

*- -Z1 XS + xk XS + 9 tw< x2 + z2 )ixs-* a rr dA 
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Integrating over the blade section yields the result for the total moments due 
to bending and elastic torsion: 

- /section 20 rr ^ 

<Vea ■ /section XO rr dA 

M r - GJ 04 + / sectlon <* 2 + z 2 >S>»° rr JA 

To M r has been added the torsion moment GJ04> due to shear stresses pro- 
duced by elastic torsion. These moments are about the elastic axis. For 
bending, it is more convenient to work with moments about the tension center 

X C : 

" J zo rr dA 
M z = 1 ( x “ x Q^ a rr dA 

Substituting for a rr and integrating yields the following moments: 

“x - E W*i + 9 '* Z > * 0 ' e i EI ZP 
“s - EI xx<*i - «’*X> + «'«; ei xp 

M r - (gj + k|T + e' 2 Ei p p) e; + e; w k|T 

+ 9'[EI x p(+i - 9’* x ) - EI zp (+;+ 6'* 2 )] 

where 

X ZZ " W dA 

X xx = J (x “ X C )2 dA 
Ip = k^A = J(x 2 + Z 2 )dA 

I xp = J(x - x p )(x 2 + Z 2 )dA 

I zp = J*z(x 2 + Z 2 )dA 

!p p = /(x 2 + Z 2 - k£) 2 dA 
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The integrals are all over the tension-carrying elements, of course (i.e., 
modulus weighted). The tension T acts at the tension center hence the 

bending moments about the elastic axis may be obtained from those about the 
tension center by (M Z ) EA = M z + x c T and (M X ) EA » The bending/torsion 

structural coupling is due to EI^p and EI Z p. For a symmetrical section, 

EI zp = 0. 

2.1.5 Vector formulation - Define the section bending moment vector 
and the flap/lag deflection w as: “ 


m^ 2) = + m z £ 

->■ -> 
w = z o i - x o k 


Mg is not^quite the moment on the section because M x and M z are really 
he ljjg and k^g components of the moment.) The derivatives of w are 


< z o 1 ~ x o k) ' = (z o “ x o 0,) ^ “ (x o + z o 0,) ^ 
= + <|) z k 


(z G i - x 0 k)" « (<!»; + e'<t> z )i + (d»; - 0»+ x )S 


Chen the result for the bending and torsion moments may be written: 

~~ ^w 11 + EI xx kk > * ^o 1 ” x o k ^ + ^tw^e^ E ^XP k ” ^ZP 1 ) 

Mj. = (GJ + kpT + Ol pp )f>4 + 6J u k|T + - EI zp t) • (Z 0 t - x 0 ty 

This is the result sought here, namely, the relation between the structural 
moments and the deflections of the rotor blade. 

Writing the bending stiffness dyadic as El = EI zz xi + EI xx lcic, and 
neglecting (for this paragraph only) the bending/torsion coupling terms (EI ZP 
and EIxp) gives 

^ (2) = EIw" 




“r = GJ eff 0 ; + ^ T0 tw 


In this form, our result appears as a simple extension of the engineering beam 
theory result for uncoupled bending and torsion (for 0 £ w = 0). The vector 
form allows a simultaneous treatment of the coupled inplane and out-of-plane 
bending of the blade, with considerable simplification of the equations as a 
consequence. 
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This relation between the moments and deflections is a linearized result. 
Thus the vectors x and k appearing in El and in w are based on the trim 
pitch angle 0=0°+ 0 tw . The perturbations of x and £ due to the elastic 
torsion give second-order moments, which have already been neglected in the 
derivation. The net torsion modulus is 

GJ eff - GJ + k p 2T + e « EI PP 

where T = f2 2 J r 1; pm dp is the centrifugal tension in the blade. For the 
elastic torsion stiffness characteristic of rotor blades, the GJ term 
usually dominates. The kp 2 T term is only important near the root for blades 
that are very soft torsionally. The 9tw EI PP term is important only for very 
soft, highly twisted blades. 


2.2 Inertia Analysis 

This section derives the inertia forces of a helicopter rotor blade. The 
blade motion considered includes coupled flap/lag bending (including the rigid 
modes if the blade is articulated), rigid pitch, elastic torsion, gimbal pitch 
and roll (which are dropped from the model for articulated and hingeless 
rotors), and the rotational speed perturbation. The geometric model of the 
blade and hub includes precone, droop, and sweep; pitch bearing radial offset; 
feathering axis droop and sweep; and torque offset and gimbal undersling. 

2.2.1 Rotor geometry- Consider an N-b laded rotor, rotating at speed £2 
(fig. 5). The mth blade is at the azimuth location: 

i|> = + mAip , m = 1 , . . . , N 

where AtJi = 2ir/N, and Ui = fit is a dimensionless time variable. The S 
coordinate system (ig, jg, kg) is a nonrotating, inertial reference frame. 

The S system coordinates are the rotor shaft axes when there is no hub 
motion. When the shaft moves, however, due to the motion of the helicopter or 
the wind tunnel support, the S system remains fixed in space. The B sys- 
tem (ig, jg, kg) is a coordinate frame rotating with the mth blade. The 
acceleration, angular velocity, and angular acceleration of the hub, and the 
forces and moments exerted by the rotor on the hub are defined in the non- 
rotating frame (S system). Figure 6(a) shows the definition of the linear 
and angular motion of the rotor hub; figure 6(b) shows the definition of the 
rotor forces and moments action on the hub. The rotor blade equations of 
motion are derived in the rotating frame. 

Figure 7 shows the blade hub and root geometry considered (undistorted) . 
The origin of the B and S systems is the location of the gimbal. For artic- 
ulated or hingeless rotors, where there is no gimbal, this is simply the point 
where the shaft motion and hub forces are evaluated. The hub of the rotor is 
a distance Zg A below the gimbal (gimbal undersling, which is not shown in 
fig. 7). The torque offset x-p^ is positive in the -ig direction. The 
azimuth is measured to the feathering axis line (ij|s projection in the 

hub plane), so the feathering axis is parallel to the jg axis and offset 
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X FA from the center of rotation. The precone angle 6p* gives the orienta- 
tion of the blade elastic axis inboard of the pitch bearing with respect to 
the hub plane; Sp^ is positive upward, and is assumed to be a small angle. 
The pitch bearing is offset radially from the center of rotation by rp^. The 
rigid pitch rotation of the blade about the feathering axis occurs at rp^. 

The droop angle dp^ 2 and the sweep angle dp^g occur at rp^, just outboard 
of the pitch bearing; 6p^ 2 an< * ^FA$ give the orientation of the elastic axis 
of the blade outboard of the pitch bearing, with respect to the precone. Both 
<$Pa 2 and dp^ 3 are assumed to be small angles; dp^ 9 is positive downward and 
6FA3 i s positive aft. Feathering axis droop dp^“ and sweep dpAs define the 
orientation of the feathering axis with respect to the precone; ^FA4 i s posi- 
tive downward, 6p^ is positive aft, and both are small angles. If 
dpA 4 = dpAg = 0, then the feathering axis orientation is just given by the 
precone; if dp^ = dp^ 2 and dp^ 5 = dp^ 3 , then the orientation is the same 
as the outboard elastic axis. 


In summary, the blade root is underslung by zp^ and offset by xp^ 
relative to the gimbal. From the root to the pitch bearing, there is a shank, 
of length rpA, which undistorted is a straight line at an angle dp^ to the 
hub plane (small precone). The blade outboard of the pitch bearing at rp^, 
undistorted, has a straight elastic axis, with small droop and sweep (dp^ 2 and 
dpA 3 ) . The feathering axis also has small droop and sweep with respect to the 
precone (dp^ and dp^ g ) . The shank (inboard of the pitch bearing at rp^) and 
the blade (outboard of rp^) are flexible in bending. The shank is assumed to 
be rigid in torsion; the blade outboard of the pitch bearing is flexible in 
torsion as well as bending. There is rigid pitch rotation of the blade about 
the pitch bearing, which takes place about the local direction of the feather- 
ing axis at r F ^, including the bending of the shank. Incorporation of the 
bending flexibility of the blade inboard of the pitch bearing means that gen- 
eral rotor configurations may be considered — an articulated rotor with the 
feathering axis inboard or outboard of the hinges or a cantilever blade with 
or without flexibility inboard of the pitch bearing. The special case of a 
rigid shank can be considered as well, of course. 

Figure 8 shows the undeformed geometry of the blade. The description of 
the blade for the inertial analysis parallels that for the structural analysis 
(see fig. 2 and section 2 . 1 . 1 ). It is assumed that an elastic axis exists, 
that the undeformed elastic axis is a straight line, and that the blade has a 
high aspect ratio, so engineering beam theory and lifting line theory are 
applicable. Here xj is the locus of the section center of gravity, x^ is 
the locus of the section aerodynamic center, and x^ is the locus of the 
section tension center. The distances xj, x^, and xq are positive aft, 
measured from the elastic axis; generally, they are a function of r. The 
corresponding z displacements are neglected. 


The i Q , j Q , and k Q coordinate system is the elastic axis /principal 
axis system of the section. Subscript o refers to the undeformed frame, 
that is, with no elastic torsion in 0, or gimbal or rotor ^peed degrees of 
freedom. The direction of the undeformed elastic axis is j Q ; i 0 and k Q are 
the directions of the local principal axes of the undeformed section. The 
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spanwise variable is r, measured from the center of rotation. This variable 
is dimensionless, so r = 1 at the blade tip. The section coordinates x and 
z are mass principal axes, with origin at the elastic axis. It is assumed 
that the directions of the mass principal axes and the modulus principal axes 
are the same. The CG is at z - 0 and x « xj. The section mass, center 
of gravity position, and section polar moment of inertia (about the elastic 
axis) are, by definition, then as follows: 

dm = m 

section 

. z dm = I ^ . xz dm = 0 
section * section 

. x dm = x_m 
section I 

J (x 2 4 z 2 )dm » I 

•'section 0 


The blade pitch angle is 0 (at this stage in the analysis, the undis- 
torted or mean pitch, denoted by subscript m) . The angle 0 is measured 
from the hub jilane to the section principal axis. It is thus the angle of 
rotation of i Q and k G from the hub plane axes. The undeformed pitch angle 
consists of the collective pitch 0 C oll plus the builtin twist 0 tw (r) : 

0 - 0 m “ B C oll + 9 tw* define B C oll as t ^ ie P itc ^ at r FA> so 
0 tw ( r FA + ) ” 0- The toot pitch is then 0° = 9 co ll* The rotation, by 9 co ll 
is not present inboard of rp^> but there can be pitch of the local principal 
axes with respect to the hub plane, which is included in 0 tw for r < rp^* 
Note that 0 tw (rp^) is not necessarily zero, hence there is a jump in 0 m at 
r-p^, of magnitude: 


e(r FA> - 6(r FA> ' 


coll 


*tw^ r FA^ 


The trim pitch angle is then: 


coll + e tw (r) » 

r 

> r FA 

° = e coll • 

r 

< 

U 

11 

tw (r) 

r 

< r FA 


It is assumed that 0 m is steady (constant in time), independent of ij>. 

Cyclic variations in 0, as may be required to trim the rotor, are included in 
the perturbation to the pitch angle. We shall alow the trim pitch angle to be 
large, hence 9 co n and 0 tw may be large angles. 

The droop and sweep of the blade elastic axis are defined with respect to 
the hub plane axes, so it follows that unless the feathering axis is parallel 
to the outboard elastic axis, these angles vary with the root pitch of the 
blade. Let and <5p^ be the droop and sweep of the blade when the pitch 
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angle at 75-percent radius is zero. Then the following relation can be 
derived from the root geometry: 


S I& 2 S FA 4 + <l5 FA 2 “ 5 FA 4 )cOS 9 75 + ^ 6 FA 3 " 5 FA 5 ^ Sln 6 75 


S FA 3 " S FA 5 - (5 FA 2 - 6 FA 4 )sin 6 75 + < S £a 3 ' S FA 5 >COS 9 75 


where 0 75 = 0° + 0£ W (r = 0.75). The angles and §fAq are fixed geo 

metric constants, so the variation of the droop and sweep due to blade pitch 
perturbations is 


6 fa 2 6 <6 fa 3 5 fa 5 ^ 


s fa 3 - - < ’’«fa 2 - V 


Between the B coordinate system (rotating hub plane axes) and the o 
system (undistorted ^ection axes), there are the following rotations: 

SpAi “ ^FA? a hout ig (small precone and droop) , about kg (small 

sweep), ana then rotation 0 m about jg^ (the large pitch angle). So 

x n = cos 0 m i B - sin 0 m £ R + j T} [(6 r , A - S RA )sin 0 m - 6 RA cos 0J 


m B ’ J B l v FA 


FA, 


m FA, 


'm J 


^o = sin Vb + COS 0 m^B + J B [_(<S FA 1 " 6 FA 2 )c ° S 9 m “ 6 FA 3 sin 9 m ] 


^EA ^ T5 ’ A ^ 


FA 3 B 


FA, 


FA 2 B 


where 6 


FA 2 


and 


6pA 3 are based on 0^ = 0 CO H> and are absent for r < rp^. 
Subscripts'" o and m will be dropped when it is obvious that the undistorted 
geometry is being considered. 


2.2.2 Rotor motion - The rotor blade motion is described by the following 
degrees of freedom: 


(a) Gimball pitch and roll motion of the rotor disk (omitted for articu- 
lated and hingeless rotors) 

(b) Rotor speed perturbation 

(c) Then torsion about the elastic axis, and rigid pitch motion about 
the feathering axis 

(d) Followed by bending deflection of the elastic axis, including rigid 
flap and lag motion if the blade is articulated. 

Figure 9(a) shows the gimbal motion and rotor speed perturbation in the non- 
rotating frame. The gimbal degrees of freedom are Bqq and $ G g — respectively, 
pitch and roll of the rotor disk in the nonrotating frame. The rotor rota- 
tional speed perturbation is ^ s . The degree of freedom i(; s is a rotation 
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about the shaft axis kg, so the azimuth angle of the mth blade is really 
Tj> m + i(j s . Figure 9(b) shows the gimbal motion in the rotating frame. The 
degrees of freedom are g G and 0 G , given by 

$G - e GC cos + e GS sin 
0 G = ” e GC S±n ^m + B GC COS 

The gimbal effects are primarily d^e to the flapwise rotation about the 

x B axis; 0g, the rotation about jg, only introduces a translation of the hub 
due to ZpA and x-p^. The blade pitch 0 is defined with respect to the hub 
plane, so only the blade inboard of the pitch bearing sees the pitch rotation 
due to 0q. 

Figure 3 showed the geometry of the deformed blade. The blade deforma- 
tion is described by twist 0 about the elastic axis, bending deflection x Q 
and z Q of the elastic axis, and rotations of the section <j> x and c|> z due to 
bending^ The pitch angle 0, including perturbations, is implicit in the it, 
and k coordinate system; x and k are the principal axes of the blade with 
no bending, but now include the blade elastic torsion and rigid pitch motion. 
The XS axes (x^g, 3XS» ^Xs) are the section principal axes and elastic axis 
of the deformed blade, including both torsion and bending. The tangent to the 
deformed elastic axis is jxs* From section 2.1.2, the rotation of the cross 
section by <f> x and <j> z is related to the bending as follows: 

(j> x + <j> £ = (z ' - x 0 ' ) x - (x' + z 0')k = (z i - x k) ' 

T x z o o' o o' s o o' 


The blade position, relative to the root, is then 

? = (r + r Q )j + x Q i + z Q t. + xi xs + zt xs 

= (r + r Q + x<|> z - z<l> x );f + (x Q i + zjx) + xi + z£ 

The perturbation of the radial position, r Q + xc() z - z-<J> x , will be neglected 

since it is much smaller than the radial position r. 

The blade pitch angle 0 is the angle of the major principal axis of the 
section (x axis) measured from the hub plane. The pitch is composed of the 
root pitch 0°(ij>) (the blade pitch at the pitch bearing, r = due to con- 

trol commands, control system flexibility, and kinematic coupling); the 
builtin twist 0 tw (r) (where = 0) 5 an < i torsion about the elastic 

axis 0 e (r,i|O (where 0 e (rp^,^) = 0; only 0 e produces shear stress in the 
blade). The blade shank inboard of r does not have the root pitch 0° or 
the elastic torsion 0 0 . Thus the blade pitch is 


e° + e tw + 0 e , 


r > r 


0 ={ 0 ' 


7 tw 


FA 

r = r 

FA 

r < r FA 
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The commanded root pitch angle is defined as 0 C = 0 CO 11 + ®con* Here ®coll 
is the trim value of the collective pitch, which may be large but is assumed 
to be steady in time; 0 con is the perturbation control input (including 
cyclic to trim the rotor), which is time dependent but is assumed to be a 
small angle. The blade root pitch commanded by the control system is 9 C ; 0° 
is the actual root pitch. The difference (0° - 0 C ) is the rigid pitch motion 
due to control-system flexibility or kinematic coupling in the control system. 
Hence the blade pitch may be written as 


0 =< 


(0 .. + e ■ ) + (e° - e c ) + 0 +0 , 

coll tw' con e 

0 ° = 0 „ + ( 0 ° - 0 C ) + 0 , 
coll con 


r > r. 


FA 


r = r. 


FA 


V 0 tw 


r < t 


FA 


The pitch angle 0 may now be separated into trim and perturbation 
contributions: 


0 =< 


0 


0 ° + 0 ' 
m 


m 


r > r. 


FA 


r = r 


FA 


r < r 


FA 


where the trim terms are (as above) 

/ 

0 


6 m=< 


coll + 0 tw 9 

r > 

r 

FA 

coll 

r = 

r FA 

tw 

r < 

r FA 


and the perturbations are 


0 


( 0 ° - 0 ) + 0 + 0 . , 
' con e ’ 

e° = (9° - e c ) + e 

con 


0 


r > r. 


FA 


FA 


r < r 


FA 


The trim value of the pitch 0 m is composed of 0 C oll and e twv it is a large 
steady angle. The perturbation of the pitch angle 0 is composed of the 
blade motion terms (0° - 0 C ) , 0 con , and 0 e ; all are small angles, so 0 is 
small. For the rigid pitch degree of freedom, the notation p Q is used where 


P = 0° = (0° - 0C) + 


con 


17 



(The notation p Q is chosen to be consistent with that for the modal expan- 
sion of the elastic torsion 0 e described below.) Note that p Q is the 
total rigid pitch motion of the blade, including the control angle 0 con . 

2.2.3 Coordinate frames- Table 1 summarizes the coordinate frames used, 
and the axis rotations between them. The unit vectors of the B system are 

i_ = sin ip i_ - cos ip j 0 
B r m S r m J S 

— -f- -y 

i „ = cos ^ + sin ip 

J B m S m J S 

S B ' 6 S 

Between the B system and the blade system, there are the following rota- 
tions: first, + d FAl ~ ^FA 2 a hout ig, and \p s - d FA g about £g; then 0 

about j EA . Hence the unit vectors are 

1 = cos 0t B - Sin 0^ B + 3g[(e G + « FAi - 6 F A 2 > sln 9 + ( * s " 6 FA 3 )coS 9] 


t = sin Ql B + cos 0 £ b + j B [-(B G + S FAi - « FA ^)cos 0 + 0P S - « FA ^)sin 0] 


X 1 


J = ^EA = ^B - ( *s - 6 FA 3 )t B + (P G + 6 FA 3 " ^FA^B 


The unit vectors of the XS are 


-f -f , . *t* 

■xs ■ 1 + V 


J xs -1 - M + ■ J + + z o»' 


^ I "t" 

k XS - k - 


For the undisturbed blade system, the rotations by an< * are 

dropped, and also the pitch perturbations in 0. Hence the unit vectors are 


i - cos 0 m iT> “ sin 0 m k + - 6 . )sin 0 m cos 0 m ] 


m B 


m B 


'B lv FA 


FA„ 


m FA, 


*o = sin Vb + COS 6 m^B + Jb [_( 6 FA 1 " 6 FA 2 )cOS 9 m " 6 FA 3 S±n ®m ] 


j o “ h + 6 FA 3 1 B + (6 FA x “ 6 FA 2 )k B 


18 



TABLE 1*- SUMMARY OF COORDINATE FRAMES 


Coordinate frame 


S system 

nonrotating, hub 
plane frame 


B system 

rotating hub plane 
frame, mth blade 


H system 

hub frame 


FA system 

blade inboard 
elastic axis 


Axis rotations 



EA system 

blade outboard 
elastic axis 


Blade system 

principal axes , 
including torsion 


XS system 

principal axes, 
including torsion 
and bending 


- 90° about k n (shaft rotation) 
m S 


S G about ig (gimbal) 

e Q about jg (gimbal) 

ijj about (rotor speed perturbation) 

SB 


6 about i (precone) 

rA^ n 


-<5 F A2 about ^FA ( dro °P) 

-6 pa 3 about kp^ (sweep) 

0 about J_, A (pitch/torsion) 
BA 

-0 G about J EA (gimbal) 


<j> x about i (bending) 

(j> about k (bending) 
z 


Now since the blade motion (0, and ^ s ) is small, the blade system unit 

vectors can be expanded in terms of those of the undisturbed frame: 



<$FA )]sin 0 + |> s + 0°(6p A 


FA, 


)]cos 0 


k 


= k. 


+ §i„ + 


■I- 


[b g- 


0 ° (6 


FA, 


” 6 FA ^ C0S 9 + + ^ ° ( d TT a - )]sin 


FA 0 


FA, 
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It follows that 


(x 0 i + z Q t) = (x 0 i 0 + z Q t Q ) + 0(z o t o - x/ Q ) + 1 b {[* 8 + 0' 
" ^G " §0(6 FA, ’ 6 FA s )] ^B ) * (x o ? o + Z oK ) 


(6 fa 2 “ 6 fa 4 )]1 b 


which is an expansion of the bending/torsion deflection of the blade in terms 
of the undisturbed axis system. 


2.2.4 Blade position, velocity , and acceleration- The distance from the 
gimbal to a point on the blade section is 


*)■ 

r 



x 1 
X FA X H 


FA 





which may be written: 


r ^B^"" X FA "" Z FA^G "" r FA^FA ^ ” x -rrA^e) 


B FA G FA s 

->■ 4 1 


+ k B (-Z FA + X FA 0 G + r FA 6 FA. } + rj + ^o 1 + z o k) + (xi + zk) 


= ^B ^"" X FA " Z FA 6 G " rlp ~ + (r " 1 


FA FA C 


+ k B [ “ Z FA + X FA 8 G + r(6 G + 8 FA 1 ) - <r - r FA )6 FA 2 ’ 

+ V r + Z FA S G “ ’‘faV + 'V + 2 o k) + + zk) 

The velocity of a point on the blade relative to the rotating frame 
(B system) is 

V r = (dt r ) B = 1 B (_Z FA®G - r ^s ) + j B (z FA®G “ X FA^'s ) + k B (x FA®G + r ^G ) 

" (r “ r FA^®° ^ ^FA _ 6 FA + (6 FA 3 " 6 FA 5 )k B 1 
2 4 3 5 


+ [(x o + x)x + (z o + z)£] 

where 

[(x 0 + x)i + (z 0 + z)£j* -a (x 0 i Q + z 0 ^ 0 )* + 0[(z o + z)i Q - (x Q + x)it Q ] 

Finally, the acceleration of a point on the blade relative to the rotat- 
ing frame and neglecting the squares of velocities is 
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■ (^r) B - V'Vc - r *»> + V Z FA S G - *FA*s> + e B (x FA § G + rS G> 
- (r - r FA )9“[(«FA 2 - ^A.^B + <5 FA 3 - S FA 5 )6 B ] 


+ [(x Q + x)i + (z Q + z)k] 


[ (x n + x)i + (z_ + z)t]" = (x t + zi J + 0[{z_ + z)i - (x + x)£] 


The acceleration of the blade is required with respect to an inertial 
frame, specifically, the S system. The B system rotates at a constant 
angular velocity ft = ftk B with respect to the S frame. The shaft motion is 
composed of linear and angular displacement of the origin of the S frame 
(the gimbal point at the hub center of rotation). The acceleration, angular 
velocity, and angular acceleration of the S system, with respect to the 
nonrotating inertial frame, are 


■+ .. -f •• *f* •• r>“ 

a o - Vs + Vs * z h k s 

-> • ■+ • -f • ry 

"o ■ Vs + Vs + Vs 

% • + “yj S + 

It is assumed that a Q , 0 > o , and 0 ) o are all small quantities. 


The motion of the blade^relative to the B frame was derived pre- 
viously — the acceleration (a r ) and velocity (v r ) of the blade. Now the 
acceleration of a blade point in inertial space is derived in terms of the 
motion of the shaft, the rotation of the rotor, and the blade motion in the 
B frame. From the result for the^acceleration in a rotating coordinate frame 
(the S frame, rotating at rate oj 0 ) , it follows that 


where a r>s and v r>s are the acceleration and velocity relative to the S 
frame. The B system rotates at angular velocity = fikg with respect to 
the S frame. Hence, with ft constant and no angular or linear acceleration 
of the B frame with respect to the S frame, it follows that 


a = a + 2ftxr + S^x(^xr) 
r,s r 


v 


r ,s 


= v 


ftxr 


where a r and v r are the acceleration and velocity relative to the B frame. 
Thus, 
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a Q + a r + 2^xv r + ftx(ftxr) + 2 ( 0 Q xv r + 2w Q x(fixr) + W Q x((o o xr) + (5 Q xr 


To first order in the velocity and angular velocity, this becomes, finally. 



*f 

+ a 4 * 2 fixv 
r r 


ftx(ftxr) + 2(o o x(fixr) + W Q xr 


The six terms in a are, respectively, the acceleration of the origin, the 
relative acceleration in the rotating frame, the relative Coriolis accelera- 
tion, the centrifugal acceleration, the Coriolis acceleration due to the angu- 
lar velocity of the origin, and £he angular acceleration of the origin. In 
dyadic operator form, and with ft = ftkg, the acceleration is 

i = j o + i r + 2Si( j B t B - y B >; r - a > ( y B + y B >? 

+ 2ft(iLr - rlc )w - (rx)m 
B B o o 


To obtain then the total acceleration of the blade, the acceleration is 
multiplied by the density of the blade point (dm dr) and integrated over the 
volume of the blade. 

2.2.5 Force and moment equilibrium- The equations of motion for elastic 
bending, torsion, and rigid pitch of the blade are obtained from equilibrium 
of inertial, aerodynamic, and elastic moments on the portion of the blade 
outboard of r: 


where Mg is the structural moment on the inboard face of the deformed cross 
section (so -Mg is the external force on the outboard face) ; is the 

total aerodynamic moment on the blade surface outboard of r; and Mj is the 
total inertial moment of the blade outboard of r. The structural moment Mg 
is obtained from engineering beam theory for bending and torsion (section 2.1) 
from the control system flexibility for rigid pitch, or from the hub spring 
for gimbal motion. Alternatively, Mg may be viewed as the force or moment on 
the hub due to the rotor (so -Mg is the force on the rotor) ; Mj is the 
inertial moment of the blade outboard of r, about the point r Q (r): 

M j = J J [r(p) - r Q (r)]xa dm dp 

*r " sect ion 

For bending of the blade, engineering beam theory gives 

^ - V + - <«xs + SxA 

- y y y y. . 

Therefore, the operator (iixs + kkxg) applied to Mj and M^ also. For 
bending, the moments about the tension center (x = kq) are required. Then the 
desired partial differential equation for bending is obtained from 3 2 M^ 2 V9r 2 
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For elastic torsion, engineering beam theory gives M t 


•t 

= 


So this 


lr E J XS ' “e* 

same operator is applied to Mj and M^. For torsion, moments about the sec- 
tion elastic axis (x = 0) at r are required; also, elastic torsion involves 
only the blade outboard of rp^. The desired partial differential equation 
for torsion is then obtained from 3M r /3r. The equation of motion for the 
rigid pitch degree of freedom p Q 


about the feathering axis, M- 


is obtained from equilibrium of moments 

, v __ 


FA 


= e 


r = 


the feathering axis (x = 0) at 
feathering axis, including perturbations due to blade bending: 


p^ * M(rp^). Here M is the moment about 
r F ^ , and ep^ is the direction of the 


'FA 


-> 
= 3 


FA 


+ (x 


i + z k) ’ 

3 O 


FA 


X** + 


X ±B 


The elastic restraint from the control-system flexibility gives the restoring 
moment about the feathering axis, completing the desired equation of motion. 

The equations of motion for the gimbal degrees of freedom 6 qq and Sgg 
are otjtain«|d from equilibrium of moments about the gimbal, M x = ig • M and 
My = jg • M, where M is the total moment (from all N blades) about the 
gimbal point, in the nonrotating frame. The equation of motion for the rotor 
speed perturbation degree of freedom is obtained from equilibrium of 

torque moments Q = -M z = kg • M, where, again, M is the total moment about 
the gimbal point. 

The total rotor force and moment on the hub (at 
obtained from a sum over the N blades of J( m ) and 
moment due to the mth blade: 

N f . 

F = f; ? (m) 
m=l 

... , . .,]■ "" N , . 

. s = f; M (m) 

m-1 - 

Since -F^ and -M^ are the forces on the blade, from force and moment 
equilibrium of the entire blade, it follows that 


die gimbal point) are 
M^ 111 ) , the force and 


-t (m) + F a = Fp 
-M (m) + M a = M x 


The hub force and moment are required in the nonrotating hub plane frame (S 
system); the components are defined as follows (see fig. 6): 

->* ->■ 

F = Hi s + Y jg + Tk g 
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fi - M x i s + M y J s - Q E s 


Note that M produces the gimbal and rotor-speed perturbation motion if those 
degrees of freedom are used, but it is also transmitted through the gimbal to 
the helicopter body or support. 


The aerodynamic forces and moments on the blade are obtained from the 
integral over the span of the aerodynamic forces and pitch moments on the 
blade section. The forces acting on the section at the elastic axis are F x , 
F z , and F r (see fig. 10). These are the components of the aerodynamic lift 
and drag forces in the hub plane axis system (B frame) — F x is in the hub 
plane, positive in the drag direction; F z is normal to the hub plane, posi- 
tive upward; and F r is the radial force, positive outward. There are also 
radial components of F x and F z due to the tilt of the section by blade 
bending; here F r is just the radial drag force. Thus the aerodynamic force 
acting on the section at the deformed elastic axis is 


f =F? +f1 - IT • (Fi + F t ) + F t 

aero x B z B :1 B J XS v x B z B' r J XS 

= F i. + F L + F j . 
x B z B r J B 


where 


K ' F r - "A + V X + • <*„? + z o E > ’ > 


- F [-^ + 

x l T s 


FA C 


+ t 


B 


(x i + z k) ' ] 
o o 


Finally, M a is the section moment about the elastic axis, positive nose-up. 
Thus the aerodynamic moment is M aero = M a J^s • 


2.2.6 Bending equation- The equation of motion for blade bending is 
obtained from 

'•£a< 2 > + 

9r^ I! 9r z I 9r z A 

where M is the moment about the tension center (x = xq) at r and 

= (ii xg + S xg )M = [it + ££ - (x Q i + z Q k) ' j"]M 


Considering first the blade outboard of r F ^, the inertia moment is 

-v r 1 r ^ 

M, 


■u 

-£/.{ 


(r - r )xa dm dp 

. . ' pxz 'rx n o 

section C 


(p - r)j + (x q + x)i + (z q + z)£ - [(x q + x c )i + z Q k] | r 


xa dm dp 
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So 


|f = -[j + (x o i + z o k + x c i)']x jf / 
= j”xj*a dm - j^(x o i 


a dm dp 


3 2 S 

Jr 1 


+ z k + x„i ) ' x 
o C 


-Ja 
•]' 


4- zS - x^T)xa dm 


a dm dpi 




-> -> ■> 

(xi + zk - x^i) xa dm 




t • H = f f [(z + z)l - (x + x)£ - (z t - x Jc - xJc) | ] • a dm dp 

J J o o o o u r 


Finally , 


£r Bf > - (It + K) -^r - [ C*t + ') • Sj]" 

= jx ^a dm + £j*(zi - xk + x^k)j • a dmj 


[(z o t - xjc - x^k) ’ J" J" J * a dm dpj 

{< 


- |(x i + z k) f 
■ o o 


a 


[(.z + z)i - (x q + x)k 


- (z q 1 - x Q k - x c k) | r J * a dm dpj 


The last term in this result — [ (x Q ^ + z 0 k)';j * M l]" ~ will be neglected since 
it is order (c/R) 2 smaller than the first term. Including the case r < r-p^, 
which introduces only an effect of droop and sweep, the result is 


3 2 M (2) 
~9r 


v(2) 

Y — = Jx J*a dm + (zi - xk 4 x G k) J • a dmj 1 


T(z i - x k - 

L ° ° 

k 


x c k) ’ 


5(r ~ r FA ) ( 6 FA 2 1 B + 6 FA 3 k B 


J" J 4 * a dm dpj 

’ )/:; 


j • a dm dp 


where 5(r) is the delta function, that is, an impulse at r = 0. The 


acceleration due to the shaft motion (with r = rjg) is 
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**■ 4- 

a = a Q + 2ft(k B r - rk B )oo 0 - rxm 0 

“a 0 + 2 ^r ( y B - 


S 2 Sr( 2 ) 


= 4 B xa = ill - tl)a 0 + 2to(U )w 0 + mr(t t + 


The blade relative acceleration gives 

3 2 S^ 2 ^ _ /_ « , v \ 


- i x SK dm = m pB (z FA®G +r ^s )+ ^B (x FA § G + rS G )+ 


- 0[ (x 0 + Xl )i + z 0 £] - 0° (r - r pA ) - 6^ 


) t B - ( 6 M 2 - X^]} 




The centrifugal acceleration is a = -^ z (igi B + jBJB^ r > so 

Jxt - 3 = t B |-x FA - z FA e G - r FA « FA3 + i B • [ (x 0 + x)t + (z o + z)S] j + n 2 i B r (l 


S G + S FAj - S FA 


j . a = -ft 2 r 


3 2 M (2) 


= -fi 2 ^(z Q i - x Q ic) ' J" pm dpj 


+ mlLlL • (z i - x it) - { [0 (x i + z It 4- x p i) ] T 


-t . -K 


B B O O 


o o C 


pm dpi + [ (x c - x-j-) 6^r mj T - * (x Q i + zt. + Xjt) - 6 (r - r^)©* 


FA 3 ' X 


" r FArFA 2 6 FA 4 


)**■(' 

*.)'■] 


6 FA 6 FA 

^ 2 4. 


Kir 


pm dp + e’m r|6 FAg - « pA > B 


+ V Z FA 0 G _ V rS G 


- fi 2 l[(x„-x T )krm]' -| (x^k)’ j pm dp | " 6 ( r “ r FA ) ( 6 FA ^ + 6 FA k B 


C I 


pm dp + k B m(x FA +r FA 6 FA3 - *j cos 6j 




“>■ ->■ 

The Coriolis acceleration is a = 2Qk-r>xv , so 
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j ♦ a = 2Q[-r^ s - k g • (z q 1 - x o k) ] 


jxa = 2^{kgj" B • v v + 


;(h ' h) 


4 B + 


[- r* s - k B • (z 


-> 4 e I 

i - x k) ]J 

o o I 




For the Coriolis acceleration due to the radial velocity j • v r , it is neces- 
sary to include the effect of the change in the radial position of the blade 
due to bending: 


Ar 


1 J[ [ ( V + z o i * x I d) ' - (♦. - 5 FA 3 ) t B + ( e G + 4 FAj - S FA 2 ) S Bj 2dp 


so 


J * V » - f (z i - xZ) ' • (z i - x k - x T k) ' dp - (z i - x k) 

r ^ o o 001 oo 


( 6 FA x " <5 FA 2 ) i B “ 6 FA 3 k B 


_Z FA + h * (z o^ ‘ X c£ " X I S) 


+ - (r - r FA )6 pA2 


“ ( r “ r FA^ 6 FA 


J 


- # s x FA + £g . (z 1 - xZ - xS) 


Then 




9r d 


(- e b" jf 


(z o d - x o {) ’ ’ • (z j - x o k - x I i ' , ' dp - V®G -Z FA + d B 


(z o i - xt - xfi + rS^ - (r - + rSg 


^g[- 

:] - V*s[ ! 


X FA + k B 


• <V - X o^ - X I K) + Z FA 6 G + r FA S FA, 


■( 


(x c ~ x I )km[rij; s + k g 


(z o i - x o k)*]}’ - {(z o i - xZ - x ty J [p| + Z^ • (zZ - xZ)'] 


X m dpj + 6(r - r FA ) (fi^tg + ^ 

FA 

+ m ( 6 G + { FAj - ‘ x o t)- + r *Ay 


s + V ( V ' x o k) ’ I” dp 
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The structural moment (from section 2.1) is 

3 2 ^ 2) 

3r 2 - [(El^tl + - xj£)"]" 


zz 


XX 


Finally, the aerodynamic moment about the tension center (x = x c ) at r, due 
to the blade loading acting at the elastic axis at p, is 


*. • ft 


•> 


r | - r I IxF dp 

poo 'rx^oj aero 


“ /’(P * r >< E Ai ' F x £ B> d P 


So 


3 2 ftf 2) 


3r E 


- jxf = F t_ - f! 
aero z B x B 


2.2.7 Elastic torsion equation- The equation of motion for elastic tor- 
sion is obtained from 

* ^ m 3 M 3 M 

- — M - — — M = - — — M 

3r r E 3r 3r r^ 

where M is the moment about the elastic axis at r and 


— M = — - T • M - t • — + [(x 1 + z £)' • Ml’ 

3r r 3r J XS 3 3r 1 '“o o' MJ 


The inertia moment is 

.1 


'section 


(r| - r| )xa dm dp 

'pxz ,'roo 


** * l L 

“ J J*[(p “ + (x Q + x) i + (z Q + z)£ - (x Q i + z Q S)| r ]xa dm dp 


So 


_3M 

3r 


-[j+ (x o i+ Z o k) .]x jf J a dm dp - f (xi + zk) 


x a dm 


28 



Thus we have 


* J*(xlc - zi) • a dm - (zt - xjc)" • J (p - r) J*a dm dp — (x^t + Z Q k) ' 
• J*(xic - zl) J • a dm - (x^ + zJc)" * J J t ( z 0 + z )^ “ ( X D + 

- (z Q t - x o it) | r ]J * a dm dp 


The ordinary differential equation for the kth torsion mode of the mth 
blade is obtained by operating with • • )dr, where £k is the elastic 

torsion mode shape. It is most convenient to apply this operator at this 
point in the analysis: 


r 


3M 


•k 3r 


FA. 


dr = J" Jj£; k (xk “ “ J" 5 k ( z Q i ~ x o k)"(r - p)dp| • a dm dr 

r FA ' r FA * 

- J* C k |(x o t + 2 D ^) ' *J"( X ^ “ z ^) j * a dm + (x q J + z Q ic)" 

r VA ' 


FA 

.1 


J* J*[( z o + z)i - (x Q + x)k - (z Q i - X Q k)| r ]J‘ a dm dpjdr 


and the following notation is adopted: 

•r 


\ ■ 5 k x i k 


-! 

FA 


£ k (z 0 i - x o k)"(r - p)dp 


The acceleration due to the shaft motion gives 


z* 1 m r 

r 1 -> 


[ 5 k IT dr ‘ 

/ V 

dr(i B i B + k B k B )l o +2 Q J 

r FA 

r FA 

r FA 


„ • U) 

1 B o 


/ V m dr( V B " *b*B )S o 


FA 


The blade relative acceleration gives 
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/ 


3M 
. r 

’k 9r 


dr 


-r 


FA 


V> dr • (-r FA e G I B + **./«&) + f 

FA r FA 


X k rm dr • (-* S ± B 


+ -f W + / \ * ( V + Z o^" m dr + / [ \ ’ < z 0 d ‘ V 

r FA r FA 

- X I E > + W 18 - dr - T V [( S FA 2 - 

FA 

+ ( 8 FA 3 - s fa 5 ) S b] (i - r FA )§v dr - / E k Sl 6 dr 

r FA 


where Iq = J(x 2 + z 2 )dm is the section pitch moment of inertia about the 
elastic axis. The blade centrifugal acceleration gives 

J f 1 aM r ( A + -1 

I ? k TT dr ■ V dr • W'g - J V” dr • Vg 

fa 1 r FA r FA 

+ J ?; k 01 0 (cos 2 0 - sin 2 0)dr - J* • k^ • (x Q i + zjc 
r FA r FA 


+ Xji)m dr 


-/ 


L FA 


(x i + z k) • kx T m - r(x i + z k)" 
o o I o o 


' .f <Z o t ■ X o f " x / )m dp dr " | \ • £.( 

J FA u 


X FA + r FA 6 FA. 


) 



In the centrifugal acceleration, we have neglected a number of terms due to 
blade torsion and pitch which are of the same order as the propeller moment, 
but which are normally much smaller than the structural moment. 
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9M 


With T = pm dp, the structural moment (from section 2.1) is 

- ^GJ + k p V jf pm dp + e;£EI pp ^9; ' + ^ W V« 2 j[ P” A' 


+ [ °t„‘ EI XP E - EI Z P T) ’ <* 0 t ‘ X o S) " r 


Finally, the aerodynamic moment about the elastic axis at r is 


M t 


r l .> r l 

= J M a^XS dp + J [(p - r) J + ( V + Z o^ ~ ( V + Z o^ lr ]x:? aero dp 


So 


3M a 


JyQ X 


9r a J XS J XS J aero 

*r 


A 


dp 


9M a* 

r A -t 9M A . -t fv n 

ST = J xs * “aF + ( V + Z o k) • m a 


= -M - (z 
a o 


i - xt)" • J 


(p - r) (F x i B + F 2 k B )dp 


and 


9M 


1 0il 1 ] 

/ «k IT dr - - / Va dr + / V (F ~*« + F - 1 ‘« )dr 

r FA r FA k 


'FA 


x B z B 


where 


X A X k 5 k x I k 
k 


2.2.8 Rigid pitch equation - The equation of motion for rigid pitch is 

obtained from Mp^ + Mp^ = Mp^ , w ^ ere 

E X A. 


M « e 
FA FA 


* ■ JfA + (x o* + z <t ' lr FA - 6 FA 4 { B + S Fa/b_ ' 


M 


and M is the moment about the feathering axis at r = r pA* T ^ e i nert: *- a 
moment is 



a . f 1 f(ti - ?i \ 

%a J V 1X2 r FA°7 


■> 

xa dm dr 


= f f ( r - r PA^J + ( x rt + x )* + ( 2 + z )£ - (x 1 + z k) I Ixa dm dr 

J r J L ° ° r FAJ 


So 


'H ‘ £ /( 


Cz„ + z)i - (x + x)k - 1(5-. - 6 


'FA 
x (r - r. 


- [(v - ’faJ^ + (''FA, - s fa 5 ) : 'b] 


FA } ‘ ( V - X o S)l r T , ' <V* - x o iz) ' , r I ,/ r " r FA>} * 2 dm dr 

FA FA / 


( x i + z k)' ~ 5^* i„ + <$_. > c z„ + z)i - (x„ + x)k 


o ' ’ r p A :fa 5 -b fa 4 ~b 



* Jt Si 

FA 

- <V - X o^ )1 r F J + ( 6 FA 3 S " 6 FA 2 ^b) * [(Z Q +Z)t ~ (x o +x) ^ ] }j * * dm dr 


and the following notation is adopted: 


X =-(zi-xk- x T lt) + |(S_. - 6 

o o o I 11 PA 


[k ~ 6 FA 4 ) ? B + ( 6 FA 3 - 6 FA 5 ) k B_ 


(r " r FA ) 


+ (z 1 - x k) I + (z i - x k) ' I (r - r„. ) 
o o r p A O o 'r FA FA 

The acceleration due to the shaft motion gives 

.1 


^A K m dr( V B + ^B } *o - M / 


X rm dr • k-i,, • m 
o B J B o 


FA 


FA 


-Jt 


X o rm dr(k B i B - 


FA 

The blade relative acceleration gives 
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pi ^ 

J *o“ dr • ( - Z FA § G*B + WW 


Inn dr • (-ipl B + e G £ B ) 


,1 ^ 

[ X • (x i + z f) m dr - [ [I • (z t - x k - x T k) + x 2 ] 0 m dr 

J o o o J o-o o I I 


( *o'[(s- d FA>B + ( ( 




(r - r _ . ) 0 °m dr 
FA 


0 I e dr 


The centrifugal acceleration gives 


^FA 


-“’(Cl 




X m dr * i„z„ A 0 _, 4 - 
o B FA G 


X rm dr • k^ 3 ^ 
o B G 


-r 


0 I Q (cos 2 0 - sin 2 0 )dr 


x„. + r„. 6 _. ) + k^r 

FA FA FA 3 y B 


( 6 FAi - 


6 fa Jr dr + 


X C • k B k B 


-> -y -> - 
(x i + z k + x T i)m dr + 
o o I 


• (z i - x £) I j rm dr -I 

O o r„. J J 


FA r T 


- r FA (z o i - *«,£)' | r ^ - r FA [( 6 F A2 - V 


[ 5 FA 5 t B- 5 FA^B +(x o 1 +Z o S >’lr I .J 

+ 2 o k + Xjt) ' j(Z 0 i - X 0 S) |. 

J 1 ® + (x - vH]}" dt J 


(x i 4 - z k 4 * x 
o o 


Next, the aerodynamic moment about the feathering axis at r^ is 


Vxs dr + 


( [(, - 
T?A L 


)j 4 - (x l 4 - z k) . - (x T 4 - z k) *F dr 

v J o o o o 1 r_ A aero 

FA 


V ' / 

A *r T 


M a dr - J (F xS + F z^ B ) • X A dr 

r FA ° 
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where 


-f ->■ 

X. = X - x T k ; 

A o I 
o 

The aerodynamic and inertial moments about the feathering axis are 
reacted by moments due to the deformation of the control system, moments due 
to the commanded pitch angle, and moments due to feedback (mechanical or kine- 
matic) from the blade bending or gimbal motion. The restoring moment acting 
on the blade about the feathering axis is “M con , which is given by the product 
of the elastic deformation in the control system and the control system stiff- 
ness K con . Hence 

M con = K con{ § ° “ 6 con + £ ^ ± + ~ (0 1S COS K “ 0 i C Sin V^s} 


The variables are the bending degrees of freedom (introduced below) , so 

K p is the kinematic pitch/bending coupling due to the control system and 
r i 

blade root geometry. Similarly, Kp is the pitch/flap coupling for the gim- 

G 

bal motion. For the rigid flap motion of the blade, this coupling is usually 
expressed in terms of a delta- three ($ 3 ) angle so that Kp = tan <5 3 . The ipg 
term is the pitch change due to the rotor azimuth perturbation with a fixed 
swashplate. For a rigid control system (K con 00 ) , the rigid pitch equation 

of motion reduces to 

P o * 5 ° = e con ‘ iC V>i " V 8 G + ( 9 1 S C0S *„'f e ic Sln *m ) *s 

1 1 G 

So, in this limit, p Q becomes just the control input, plus the kinematic 
coupling terms. 


Now the control-system stiffness K con is written in terms of the non- 
rotating natural frequency of the rigid pitch motion of the blade, oj 0 , as 



Then the structural pitch moment is 
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2.2.9 Blade force- The net force of the mth blade on the hub is 
p(m) = - Fjj where F is the force due to the blade at the hub. The iner- 

tial force is 


^ - a 


a dm dr 


Then the acceleration due to the shaft motion gives 

.1 ,1 -1 


I = J m dr a Q + J rm dr^ - + J\ rm dr(t^ - 

A “n 


o w o 

The relative acceleration gives 


fi y 1 

* " J rm dr(-t B * g + it B 0 G ) + J (x Q t + z o k)‘ 


m dr 


The Coriolis acceleration gives 


//> 


F ~ 2ft | | k^xv dm dr 

w _ B r 
o 


= 2fi 5 B |" J* rm dr Ig + J* ig * (x Q i + z o ^)*m drj 


The centrifugal acceleration gives 

* = “ fi2 / / ( Vb + Vb^ dm dr 


X FA + (r " r FA )6 FA, m dr + ?B J ** dr “ H J rm dr 

3 J o o 


•4 ll 

+ * B J[ V ( V + Z S + X I^ )m dr j 

Finally, the aerodynamic force is 


ml ^ *1 

J rm dr - t B J 


dr 


B + F 2 k B + F rV dr 
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2.2.10 Blade moment- The net^moment of the mth blade on the hub about 
the gimibal point is M^®) = M a - Mj. The inertial moment is 


rxa dm dr 


The acceleration due to the shaft motion gives 

-1 *1 

S = J rm dr(t B ^ B - + 2 0 ■ J r 2 m dr • Z 

*o o 

+ J r : -m dr(i B t B + 

r o 

The relative acceleration gives 

rm dr( ZFA 0 G t B + x^ig) + J r 2 m dr(^ g + ) + J (2 t - x it)“ 


x rm 


dr - | 0t(x o + x F )t + z o ^]rm dr - 6°^^ - 6 pA Jt 
FA 


fa 3 ^fa 5 /-b 


- ( 5 fa 2 - 6 fa 4 ) E e] i < 

FA 


(r - r FA ) rm dr 


The centrifugal acceleration gives 


B ["' FA + ^ " VS,]™ ^ + l r2m dr 6 G + l ^ 

■ Z * 1 

• (x t + z it + x T it)rm dr + J 0it D • (z i - x k - 
o o I J B o o 


x F k)rm dr 


- |, (x-x)£ 


5/ r FA 


(r - r FA > rm dr 


) 


using the relation 
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-> 

rx 


[fix(3x?)] “■ r X ^ * 1C 


= -Q2^xr • r 


= -H 2 r(-i B kg • r) 


Now the Coriolis acceleration is 


•"> -> 


rxa = 29, |(l B xr) [r^ g + • (z Q i - xJt) ' ] + r^ B • 




where 


-f ->• 

J B xr 


■ 4 
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Finally, the aerodynamic moment is 


M. 


= / ?xf 

J o 

■r 


dr 


<Vl> - VB )r dr 


37 



2.2.11 Gimbal equation- Th_^ equation of motion for the gimbal degrees of 
freedom are obtained from the ig and j s components of the hub moment / r 

TUT y 1 ^(m) 

M = Ca M : 
m 

where M^g is the spring and damper moment at the gimbal, reacting the rotor- 
applied moments. The gimbal spring and damper are assumed to be in the non- 
rotating frame. Hence 

= L S^ K G e GS + C G^GS^ ” ^S^G^GC + C G^GC^ 

Taking the ig and jg components of M, the gimbal equations of motion are 

M y + C G^GC + K G 6 GC " 0 

+ Vos + k g 6 gs - 0 

The gimbal hub spring and damper coefficients may be written: 

K G - I r o a2 (v G 2 - » 

c g - f v c a 

R 

where I Q - S Q r2m dr, anc * V G t ^ ie rotat d n 8 natural frequency of the gim- 
bal flap motion. 


2.2.12 Modal equations- Consider the equilibrium of the elastic, inertial, 
and centrifugal bending moments. From the results in section 2.2.6, these 
terms give the following homogeneous equation for bending of the blade: 


[(EI zz -tt + EI xx ££)(z o 1 - x o £)"]" - ft 2 / pm dp (z 


->• * 
i - x k) 1 
o o 


• (z i - x i) + m(z i - x k) - 0 
' o o o o 


This equation may be solved by the method of separation of variables. Writing 


(z 1 - x it) = n(r)e ivt 


it becomes 


(Elr)")" - ft 2 J pm dp n 1 I - • n - mv 2 n = 0 
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the modal equation for coupled flap/lag bending of the_^rotating blade. It is 
an ordinary differential equation for the mode shape n (r) ; this mode may be 
interpreted as the free vibration of the rotating beam at natural frequency v. 

This modal equation, with the appropriate boundary conditions for a can- 
tilevered or hinged blade, is a proper Sturm-Liouville^eigenvalue problem. It 
follows that there exists a series of eigensolutions of this equation, 

with corresponding eigenvalues v^ 2 . The eigensolutions or modes are orthog- 
onal with weighting function m; if i f k, ' , 


dr = 0 

is possible to expand the rotor 


r 


v 


These modes form a complete series, so it 
blade bending as a series in the modes: 


t*" 

z 1 - x k 
o o 


i=l 


q ± (t) (r) 


The bending modes are normalized to unit amplitude (dimensionless ) at the tip: 
|n(l) | = 1. 

Consider the homogeneous equation for the elastic torsion motion of the 
nonrotating blade, that is, the balance of structural and inertial torsion 
moments. The results _n section 2.2.7 give 

-(gjg;)’ + i Q e e = o , 

The equation for the torsion motion of a rotating blade, including centrifugal 
forces and some additional structural torsion moments, could be used instead. 
For the torsional stiffness typical of rotor blades, these terms have little 
effect, however, and the nonrotating torsion modes are an accurate representa- 
tion of the blade motion. Solving this equation by separation of variables, f 
we write 0 e = £(r)e ia)t , so 

(GJ £ ’ ) ' + I e m 2 C =0 

This equation is a proper Sturm-Liousville eigenvalue problem, from which 
it follows that there exists a series of eigensolutions ^(r) and correspond- 
ing eigenvalues (o^ 2 (k = 1, . . . , oo ) . The modes are orthogonal with weighting 
function Iq, so if i ^ k, 

f We dr ■ 0 

r FA 
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The modes form a complete set, so the elastic torsion of the blade may be 
expanded as a series in the modes: 

00 

0 = £ P i (t)5.(r) 

1=1 1 

These modes are the free vibration shape of the nonrotating blade in torsion, 
at natural frequency The torsion modes are normalized to unity at the 

tip, 5 k (l) = 1. 

2.2.13 Expansion in modes- The bending and torsion motion of the blade is 
now expanded as series in the normal modes. By this means, the partial dif- 
ferential equations for the motion (in r and t) are converted to ordinary 
differential equations (in time only) for the degrees of freedom. 

For the blade bending, we write 


(zi-xk) = (z i - x k) . 
o o o o trim 


i=l 


q i (t)n i (r) 


where are the rotating, coupled bending modes defined above and 

(z Q i - x o^trim t ^ ie tr ^ m bending deflection. These modes are orthogonal 
and satisfy the modal equation given above. The variables qj> are the 
degrees of freedom for the bending motion of the blade. When the substitution 
for the modal expansion is made_. the subscript "trim" will be dropped, as that 
is all that can be meant by (z Q i ~ x Q k) then. 


For the blade elastic torsion, we write 


0 

e 


- £ 

i=l 


p i (t) C ± ( r ) 


where are the nonrotating elastic torsion modes. These modes are orthog- 

onal and satisfy the modal equation given above. The variables p^ (i > 1) 
are the degrees of freedom for the elastic torsion motion of the blade. The 
degree of freedom for rigid pitch motion is p c = §° = (6° - 0 C ) + 0 con . For 
rigid rotation about the feathering axis, the mode shape is simply £ 0 = 1* 
Thus the total blade pitch perturbation is expanded as the series: 


e - £ p.WiiM 

1=0 

The total blade pitch 0 (mean and perturbation) is then 


0=0 + 0 = (0 ,, + 6* ) + 2 P-£- 
m coll tw' — 

x=0 
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Subscript m on the trim pitch angle is dropped when the substitution for the 
modal expansion is made since it is no longer needed to distinguish between the 
trim and perturbation quantities. 


2.2.14 Fourier coordinate transformation - To this point in the analysis, 
the equations of motion and the rotor hub reactions have been obtained in the 
rotating frame, with degrees of freedom describing the motion of each blade 
separately. In fact, however, the rotor responds as a whole to excitation 
from the nonrotating frame — shaft motion, aerodynamic gusts, or control 
inputs. It is desirable to work with degrees of freedom that reflect this 
behavior. Such a representation of the rotor motion simplifies both the analy- 
sis and the understanding of the behavior. 

The appropriate transformation to obtain the degrees of freedom and equa- 
tions of motion in the nonrotating frame is of the Fourier type. There are 
many similarities between this coordinate change and Fourier series, discrete 
Fourier transforms, and Fourier interpolation; the common factor is, of 
course, the periodic nature of the system. A Fourier series representation of 
the blade motion is appropriate for dealing with the steady-state solution. 

Here we are considering the general dynamic behavior, including transient 
motions; hence the Fourier coordinate transformation is required. This coor- 
dinate transformation has been widely used in the classical literature, 
although often with only a heuristic basis. For example, it has been used in 
ground resonance analyses to represent the rotor lag motion (ref. 4) and in 
helicopter stability and control analyses for the rotor flap motion (ref. 5). 
More recently, there have been applications of the Fourier coordinate trans- 
formation with a sounder mathematical basis (e.g. , ref. 6). 

Consider a rotor with N blades equally spaced around the azimuth, at 
= ^ + mAty (where Aip = 27r/N and the blade index m ranges from 1 to N) . 
Here ip = Qt is the dimensionless time variable. Let q' m ' be the degree of 
freedom in the rotating frame for the mth blade, m = 1 to N. The Fourier 
coordinate transformation is a linear transform of the degrees of freedom from 
the rotating to the nonrotating frame. Thus the following new degrees of 
freedom are introduced: 


p - 1 V 
3 o ‘ N 2- 3 

m=l 

o 2 (m) . 

p = T7 V q cos nip 

nc N ^ m 

m=l 

ft _ 2 ^ (m) • i 

$ TT > _ q sin nip 
ns N ^ Y m 

m=l 

ft - I rt (m) t n m 

8 n/2 'n V (_1) 

m=l 


41 



Here $ 0 is a collective mode, and g^ are cyclic modes, and %/2 

the reactionless mode. For example, for the rotor flap motion, g Q is the 
coning degree of freedom, while g^ and 3 1S are the tip-path-plane tilt 
degrees of freedom. The inverse transformation is 

q (m) = 3 o + £ (3 nc cos wp m + (3 ns sin n^) + 6 N/2 (-l) m 


which gives the motion of the individual blades again. The summation over n 
goes from 1 to (N-l)/2 for N odd and from 1 to (N-2)/2 for N even. The 
3 n / 2 degree of freedom appears in the transformation only if N is even. 

The corresponding transformation for the velocity and acceleration are 


i w . e o + 


Sics 


n 


nc 


ng )cos xvp + 
ns r m 


(3 


ns 


- nS nc )sin n* m ] + 6 N/2 <-l) m 


5 W - 6 


+ + “ n 2 g )COS + = “ 2n S 

nc ns nc m ns 

n 


nc 


- n 2 B )sin nip ] 
ns m J 


Note that transformation to the nonrotating frame introduces Coriolis and 
centrifugal terms. 

The variables 3 0 > 3 nc > 3 ng > and Sjj /2 are degrees of freedom, that is, 
functions of time, just as the variables q( m ) are. These degrees of freedom 
describe the rotor motion as a whole, in the nonrotating frame, while q( m ' 
describes the motion of an individual blade in the rotating frame. Thus we 
have a linear, reversible transformation between the N degrees of freedom 
q (m) in the rotating frame (m = 1, . . . , N) and the N degrees of freedom 
(g Q , 3 nc > 3 ns , 3n/ 2^ *- n nonrotating frame. Compare this coordinate 
transformation with a Fourier series representation of the steady-state solu- 
tion. In that case, is a periodic function of so the motions of all 

the blades are identical. It follows that the motion in the rotating frame 
may be represented by a Fourier series, the coefficients of which are steady 
in time but infinite in number. Thus there are similarities between the 
Fourier coordinate transformation and the Fourier series, but they are by no 
means identical. 


This coordinate transform must be accompanied by a conversion of the 
equations of motion for q( m ) from the rotating to the nonrotating frame. 
This conversion is accomplished by operating on the equations of motion with 
the following summation operations: 


^ | £(...)c° S n ^ m , | £(...)sin n^ m , |2(...)(-l) m 

mm m m 
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The result is equations for the 6 0 , 3 nc , 3 ns , and 3 N / 2 degrees of freedom, 
respectively. Note that these are the same operations as are involved in 
transforming the degrees of freedom from the rotating to the nonrotating frame. 
Since the operators are linear, constants may be factored out. Thus with con- 
stant coefficients in the equations of motion, the operators act only on the 
degrees of freedom. By making use of the definitions of the degrees of free- 
dom in the nonrotating frame, and the corresponding results for the time 
derivatives, the conversion of the equations of motion is then straightforward. 
Complexities arise when it is necessary to consider periodic coefficients, 
such as due to the aerodynamics of the rotor in nonaxial flow (see sec- 
tions 2.6.3 and 4.1). 

The total force and moment on the hub have been obtained by summing the 
contributions from the individual blades. The result is operators exactly of 
the form above, for obtaining the total hub reaction in the nonrotating frame 
from the root reaction of the individual blades in the rotating frame. The 
origin of the summation operation is clear, and the sin or cos if/ m 
factors arise when the rotating forces are resolved into the nonrotating 
frame. One may, in fact, view the equation conversion operators in general as 
simply resolving the moments on the individual blades into the nonrotating 
frame. 

The Fourier coordinate transformation is often associated in rotor dynam- 
ics with the generalized Floquet analysis. The latter is a stability analysis 
for linear differential equations with periodic coefficients. Indeed, there 
is a fundamental link between these topics because both are associated with 
the rotation of the system. However, they are, in fact, truly separate sub- 
jects — either can be required in the rotor analysis without the other. For 
example, a rotor in axial flow on a flexible support (or with some other 
relation to the nonrotating frame) requires the Fourier coordinate transforma- 
tion to represent the blade motion, but is then a constant coefficient system. 
Alternatively, for the shaft-fixed dynamics of a rotor in forward flight, a 
single-blade representation in the rotating frame is appropriate, but there 
are periodic coefficients due to the forward flight aerodynamics which require 
the Floquet analysis to determine the system stability. 

For the present investigation, the degrees of freedom to be transformed 
to the nonrotating frame are blade bending, blade pitch, and gimbal motion. 

The nomenclature for the corresponding degrees of freedom in the rotating and 
nonrotating frames are as follows: 


Rotating Nonrotating 
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The notation g^ 1 ) is used for the ith bending mode in the nonrotating 
frame. With the modes ordered according to frequency, g^ 1 ) is thus usually 
the fundamental lag mode, and g( 2 ) the fundamental flap mode. Similarly, 

0'*' is the ith torsion mode, with 0(°) rigid pitch and the remaining 
modes elastic torsion. The collective and cyclic modes (0,1C,1S) are particu- 
larly important because of their fundamental role in the coupled motion of the 
rotor and the nonrotating system. When the transformation of the equations 
and degrees of freedom is accomplished, for axial flow there is a complete 
decoupling of the variables into the following sets : 

(a) the collective and cyclic (0,1C,1S) rotor degrees of freedom 
together with the gimbal tilt and rotor speed degrees of freedom and the rotor 
shaft motion 

(b) the 2C,2S, . . . ,.nc,ns , and N/2 rotor degrees of freedom (as present) 

Thus the rotor motion in the first set is coupled with the fixed system, while 
the second set consists of purely internal rotor motion. Nonaxial flow 
couples, to some extent, all the rotor degrees of freedom and the fixed system 
variables, primarily due to the aerodynamic terms; still the above separation 
of the degrees of freedom remains a dominant feature of the rotor dynamic 
behavior . 

In this section, only the case of a three-b laded rotor is considered; 
thus the collective and cyclic rotor degrees of freedom (0,1C,1S) are the com- 
plete description of the rotor motion. The equations are extended to a 
general number of blades in section 4. With four or more blades, additional 
degrees of freedom are introduced compared to the N = 3 case, while the 
two-bladed rotor requires special consideration. 

2.2.15 Equations of motion- The elements are now available to construct 
the equations of motion for the blade bending and torsion modes in the rotat- 
ing frame and to construct the forces and moments acting on the hub due to one 
blade. The following steps are required: 

(a) Substitute for the expansions of the bending and torsion motion as 
series in the modes. 

,(b) Use the appropriate modal equation to introduce the mode natural 
frequency into the bending or torsion equation, replacing the structural 
stiffness terms (and for bending also some of the centrifugal stiffness terms) . 

(c) For the bending equation, operate with Jo 1 hk* ( • • • )dr to obtain the 
ordinary differential euqation for the kth mode of the mth blade (q^ 
equation) . 

(d) For the torsion equation, operate with / r } ^(..Odr to obtain the 

FA 

ordinary differential equation for the kth mode of the mth blade (p^ 
equation) . 
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The result is the equations of motion and hub reactions in the rotating 
frame. The transformation to the nonrotating frame involves the following 
steps: 


(a) Operate on the hub force and moment with 2( . . . ) ; that is, sum over 
all N blades to obtain the total force and moment on the hub. 

(b) Find the ig, jg, and kg components of the force and moment in the 
nonrotating frame (S system) . 

^ (c) Write the shaft motion a Q , w 0 , and u> 0 in terms of the ig, Jg, and 

kg components in the nonrotating frame (S system). 

(d) Apply the Fourier coordinate transform to the equations of motion 
and rotor degrees of freedom — operate on the equations for bending and tor- 
sion with (1/N)2(...), (2/N)S( . . . ) cos ip m , (2/N)S( . . . ) sin and introduce 

the nonrotating degrees of freedom. 


Names are now given to all the inertial constants. The equations of 
motion, hub forces and moments, and inertia constants are also normalized at 
this point. The inertia constants are divided by the rotor blade character- 
istic inertia 1^ = /o R r 2 m dr, and we introduce the blade Lock number 
y - pacK-Vl^. This normalization of the inertia constants is denoted by 
superscript *. The rotating equations of motion are divided by 1^; the hub 
forces and moments are divided by (N/2)I b , except for the rotor thrust and 
torque, which are divided by NI^. With this particular normalization, the 
forces and moments are obtained in rotor coefficient form. 


The resulting hub forces and moments are as follows. (The inertial 
coefficients are defined in appendix Al.) 
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2 C / 2 C > 
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The gimbal tilt equations of motion are 


Y -7T + I o* C G *^GC + CP" " 1)6 GC “ 0 


" Y “5 T + C c g%s + C ( V - 1)b gs " ° 

Finally, the equations of motion for coupled flap/lag bending and for elastic 
torsion/rigid pitch of the blade in the rotating frame are 
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Note that the structural damping terms have been included in the bending and 
torsion equations, modeled as equivalent viscous damping. The structural 
damping parameter g s (equal to twice the equivalent damping ratio) generally 
is different for each degree of freedom. The bending and torsion equations in 
the nonrotating frame are presented later (section 2.6). (The inertial coeffi- 
cients are defined in appendix A.) 

The aerodynamic forces required are 
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2.3 Aerodynamic Analysis 


In this section, the aerodynamic forces and moments on the rotor blade 
are derived. We shall consider the general case of high or low inflow and 
axial or nonaxial flight. The aerodynamic terms in the rotor equations of 
motion and the hub forces; and moments are obtained for two cases: axial flow, 

which involves constant coefficient equations, and nonaxial flow with periodic 
coefficients. 

The principal assumptions in the aerodynamic analysis are: lifting line 

theory (i.e., strip theory or blade element theory) is used to calculate the 
section loading; the order c (rotor chord) terms in the aerodynamic lift are 
neglected; the order c 3 terms in the aerodynamic moment are neglected; vir- 
tual mass aerodynamic forces and moments are neglected; only first-order 
velocity terms are retained; and aerodynamic interference effects between the 
rotor and support are neglected. The analysis includes reverse-flow and 
large-angle effects. The effects of transient inflow changes on the system 
dynamics are also included, using an elementary model described in section 2.5 

2.3.1 Section aevodunamic forces - A hub plane reference frame is used for 
the aerodynamic forces. All forces and velocities are then resolved in the 
hub plane (i.e., in the B system). The hub plane reference frame is fixed 
with respect to the shaft:; hence it is tilted and displaced by the shaft 
motion. Figure 10 illustrates the forces and velocities of the blade section 
aerodynamics. The velocity of the air seen by the blade, the pitch angle, and 
the angle of attack are defined as: 0 is the blade pitch, measured from the 

reference plane; u T , up, and Up are the components of the air velocity seen 
by the blade, resolved with respect to the reference frame; U = (u T 2 + up 2 ) 1 / 2 
is the resultant air velocity in the plane of the section; cf> = tan~*^ Up/up is 
the induced angle; and a = 0 - $ is the section angle of attack. The veloc- 
ity up is in the hub plane, positive in the blade drag direction; up is in 
the hub plane, positive radially outward along the blade; and up is normal 
to the hub plane, positive down through the rotor disk. The aerodynamic 
forces and moment on the section, at the elastic axis, are defined as: L and 

D are the aerodynamic lift and drag forces on the section, respectively, 
normal and parallel to the resultant velocity U; F z and F x are the compo- 
nents of the total aerodynamic force on the section resolved with respect to 
the hub plane, normal to and in the plane of the rotor; F r is the radial drag 
force on the blade, positive outward (same direction as positive up) ; and 
M a is the section aerodynamic moment about the elastic axis, positive nose-up 
The radial forces due to the tilt of F z and F x are considered separately; 
hence F r consists only of the radial drag force. 

The section lift and drag are 

L - \ pV 2 cc t 

D = ~ pU 2 cc d 
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where U is the resultant velocity at the section, p- is the air density, and 
c is the chord of the blade. The air density is dropped at this point in the 
analysis while the quantities are made dimensionless by use of p, Q, and R. 
The section lift and drag coefficients, C£ = c^Ca, M) and c^ - cj(a, M) are 
functions of the section angle of attack and Mach number: 

a = 0 — cj) = 0“~ tan” 1 u p /u p 
M = ^ Ip U 

where M T pp is the tip Mach number — the rotor tip speed QR divided by the 
speed of sound. The dependence of c^ and c^ on other quantities, such as 
the local yaw angle or unsteady angle of attack changes, is neglected. The 
radial drag force is 

F r - <V U>D ' 1 U V c d 

This radial drag force is based on the assumption that the viscous drag force 
on the section has the same sweep angle as the local section velocity. The 
moment about the elastic axis is 

M a - -x A L + M a(; + Mjjg 
e 

■- X A itS + flVc, +M„ S 

e ac 

where is the distance the aerodynamic center is behind the elastic axis, 

c m is the section moment about the aerodynamic center (positive nose-up) , 
ac 

and Mjjg is the unsteady aerodynamic moment. The effective distance between 
the aerodynamic center and elastic axis is 



Normal flow 
Reverse flow 


For the section aerodynamic moment, it is necessary to include the unsteady 
aerodynamic terms, which, from thin airfoil theory (ref. 7), are 


M us 




sign(V) 


Here w is the upwash velocity normal to the blade surface 

(w = uj sin 0 - up cos 0), B = dw/dx (mainly, the pitch rate 0), and 

V = up cos 0 + up sin 0. For stalled flow, the unsteady moment is set to zero 

(M US = 0) . 
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The aerodynamic forces with respect to the hub plane axes are then 
F z = L cos <j) - D sin <f> = (Lu^ - Du p )/U 

F^ = L sin c|) + D cos <|> = (Liip + Du p )/U 

Substituting for L and D and dividing by a, the two-dimensional section 
lift-curve slope, and by c m , the mean section chord (which enter the Lock 
number y also), we obtain 
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The net rotor aerodynamic forces are obtained by integrating the section 
forces over the span of the blade and then summing over all N blades, 

2.3.2 Perturbation forces- Each component of the velocity seen by the 
blade has a trim term due to operation of the rotor in its trim equilibrium 
state and a perturbation term due to the perturbed motion of the system. The 
latter results from the system degrees of freedom and is assumed to be small 
when the linear differential equations that describe the dynamics are obtained. 
The blade pitch and section velocities are written trim plus perturbation 
terms : 

e => e + 60 


U„ 


U T + 5u T 


Up =*> u p + OUp 

U R " U R + 6u R 

It follows that the perturbations of a, U, and M are 
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6a = 69 - (u T 6Up - Up6u T )/U 2 
6U = (u T Su T + Up6Up) /U 
6M - 

and the perturbations of the aerodynamic coefficients are 

dc z dc z 

■ -sr s “ + lit “ - c i s “ + x SM 

a M 

(similarly for c m and c^) . The perturbations of the section aerodynamic 
forces may then be obtained by carrying out the differential operation on the 
expressions for F z , F x , F r , and M a , using the above results to express the 
perturbations in terms of 60, 6urp, 6up, and 6u^. The coefficients of the 
perturbation quantities are then evaluated at the trim state. The results for 
the perturbation forces are:'' 
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c 
c 

m 

2.3.3 Blade velocity- Now the air velocity seen by the blade section is 
considered: the trim velocity, composed of the forward speed, rotor rotation, 

and rotor- induced velocity terms; the perturbation velocities due to the rotor 
degrees of freedom and the shaft motion, and to the aerodynamic gust velocity. 

The rotor is rotating at constant speed The velocity of the rotor 

with respect to the air is defined in figure 11. The rotor has a steady trim 
velocity V in inertial axes and a trim angle of attack a^p of the rotor 
hub plane with respect to V. The velocity vector is in the rotor x-z plane, 
and otjjp is positive for forward tilt of the disk, producing a component of 
V downward through the disk. The rotor wake-induced velocity v is assumed 
to be uniform over the disk and normal to the hub plane (fig. 11). Following 
standard helicopter practice, the rotor advance ratio y and inflow ratio X 
are defined as 

y = V cos a H p/ftR 
X = (v + V sin ajjp)/f2R 
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The advance and inflow ratios are the dimensionless components of the rotor 
velocity in the hub plane axis system. The inflow ratio X is usually small 
for helicopter operation; the analysis is applicable to large inflow as well, 
however (as would be encountered, e.g. , in proprotor operation). The 
advance ratio y is zero for hover and axial flow, and y > 0 for helicopter 
forward flight. Note that, in body axes, the trim velocity vector is fixed 
with the reference frame and would therefore tilt with it. However, for the 
rotor and wind-tunnel support analysis, an inertial frame is used, so that 
tilt of the rotor by the shaft motion results in a small change in the direc- 
tions of X and y as seen in the reference frame. 

The rotor-induced velocity is obtained using the momentum theory result: 

^T 

X = y +. — ============== 

2 + h 2 /k ^ 2 


where y z = y tan a^p. Empirical correction factors and Kf are included 
for the effects of nonuniform inflow, tip losses, swirl, blockage, etc., in 
hover and forward flight. For the vortex ring and turbulent wake states, this 
momentum theory result is not applicable. Thus, if 

y 2 + (2y + 3X, ) 2 < X 2 

z h ' h 

the following expression is used instead: 


X = y z 4- K h y £ 


0.373y z 2 + 0. 598y 2 

V 


- 1.991 


where X^ = 'j/Cj/2 . 


The shaft motion consists of small linear and angular velocity, with 
components defined in the nonrotating frame: 


v o ■ Vs + Vs + Vs 


-> • -y • • -»• 

u> = ai c , + aj_ + ak l 
o x S y J S 


z S 


The aerodynamic gust velocity has components uq, vq, and wq (longitudi- 
nal, lateral, and vertical, respectively) defined with respect to the shaft 
axes (fig. 11); these components are the velocity seen by the aircraft and are 
assumed to be small compared to &R. The gust components are normalized by 
dividing by QR, not by V as is often the practice for airplane analyses* 

The aerodynamic gust is assumed to be uniform throughout space. Eventually, 
the gust vector will be transformed to wind or tunnel axes (see section 3.1*4), 
but shaft axes are used in the development of the rotor equations of motion. 
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The result for the trim velocity terms is 


u T = r + p sin - p cos 


♦m[ 6 FA 3 - • (z o T - x o £) Kln + C B ‘ x 0 £) trim 


“P ‘ X + h ’ ( V - x o k> irim + rS G . + v cos *m 


+ i • (z i - x k)' . 
Bo o ' trim 


im 


trim 


! fa i " 


6 fa + 6 g 

2 trim 


u R = P COS * m + (x FA + r FA « FA J + £g • .[(« 0 t - x o^.t rim - - x o^trim ] 


- A [>i ■ 


6 FA 2 + g G triin + *B V (z o i " X o k) trim 


l 


f 


+ p sin <1^6^ - k B - (z Q 
and for the trim pitch angle: 


-f , rK f 

x - x k)' . 
i o trim 


e = e ,,+e,. + e - l - S l q. . 

coll tw eye G„ . T. H i . 

J G trim x trim 


Here 0 C y C is the input cyclic pitch required to trim the rotor. For trim 
velocity, the blade bending and gimbal motion is periodic. For axial flight, 
p = 0, the trim velocities are constant; for nonaxial flow, p > 0, these 
velocities are periodic in ip m due to the rotation of the blade with respect 
to the rotor forward velocity. 


The result for the perturbations of the velocity components and the blade 
pitch, due to the rotor and shaft motion and the aerodynamic gust, is then: 

6u T = (Aa + y, + v p )cos ij> + (Aa - x, + u r )sin ip + p cos i}> (a + ^ ) 
i x n tr m y n ^ m mz s 

+ r(a z + \p s ) + 2 * h i ) + P cos if> m 2 q ± (k g • n!) 

6up = (z h - pa y - w G ) + p cos iP m f$ G + r(3 G + a x sin ^ - a y cos i|> m ) 

+ 2 q.,(t B • n ± ) + P cos 2 q ± (i B * h!) 

6u r = -(Aa x + y h + v G )sin ^ + (Aa y - ^ + u G )cos ABg- p sin ^(a z + 

+ 2 q.jicg • (n i - rn* - p sin ifyij) - it B • n^] 

60 = § = 2 P i 5 i 
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Finally, the quantities required for the unsteady pitch moment are 
V — U.J, cos 0 + Up sin 0 

B = 0 + e G + S q.i B . nl 

w + u R w = 2 p^V + u R 2 p^IV - ^g 2u R COS 6 + sin ^m cos 9 

+ (a + <p )2tL sin 0 - (a + ip )y sin i[> sin 0 - 2u 2 q.i • n! 
z s K zs in R i i 

+ S q.(p sin n! - u^i ♦ n/') 

2 . 3. 4 Rotov aevodynarnic fovees- Combining the expansions for the section 
forces and moment in terms of the velocity perturbations, and the velocity in 
terms of the motion of the rotor and shaft, we may obtain the perturbations of 
the aerodynamic forces on the blade* These are the blade forces expanded as 
linear combinations of the degrees of freedom. Giving names to the aerody- 
namic coefficients at this point in the analysis, the results for the required 
aerodynamic forces on the rotating blade are as follows. The aerodynamic 
force for rotor bending is 




— k^ldr = M 


V 


• # 

+ M [(Aa + y, + v_)cos + (Aa - x, 
q k P x h G m y h 


+ u )sin t|>T + M • (a + if» ') + M (a +ip ) 

G m q R 5 z s q R C z s 

+ M \( z h “ ya - w_) + M ♦($ + a v sin V 

q k A h y G q fc 3 G x T m 

- a cos ip) + M + ZM • q. + Zm q. 

y m q k B G q k qi 4 l 


+ I M p 
) . 

1 


Vi 1 


The radial force is 


X I F -IfF F r 

^ “ ! 0 m [ 6 

-S[- 


G + { FAl - 6 FA 2 + V ( V + Z o k) j 


+8 + 4 FA 3 + ( V + Z o k) ']} dr 
= R p [-(Aa x + y h + v G )sin ^ + (Xa y - ^ + u G )cos if-J + ^{(Xa, + y h 

+ v G )cos + (Xa y - ^ + u G )sin *_] + R-(« z + * 8 ) + R ? (a z + * 8 ) 

+ R X (z h' ^y - W G } + V®G + “x sin *m " “y cos *■> + Vg 

+ 2 R* q. + 2 R q + I R p . 
q. l q. i D; l 
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The aerodynamic force for blade torsion and pitch is 

F 


/»! M .1 /F F \ 

I — dr - ( (-t+ — 

J k ac J \ac B ac Bl 

FA r FA ' ' 


X A dr 

V 


- M [(Act + y. 
p fc y x •'h 


v G )cos + (Xa y ' *h + u G )sin *m ] + 


M 


•(a + 

V Z 


V 


+ M n rK + * s > + V x a h " va y - W G ) + M pJ ( *G + \ Sin " “y COS *-> 


V z 


m 


+ M + S M • q. + S M q. + ZM * p t + S M p. 

P k 6 G P^V P k 1i> P k Pi P l P k P ± P l 


Finally, the aerodynamic hub forces and moments are similar to the result for 
the blade bending, but with the following changes in the integrands and 
notation: 


Integrand Coefficient notation 


Flap moment rF^ 


M 



Torque 

rF 

X 

Blade drag force 

F 

X 

Thrust 

F 

z 


Q 

H 

T 


Combining the results for the exapnsion of the aerodynamic forces, and the 
expansions of the velocities, the aerodynamic coefficients can be evaluated. 
The coefficients of the degrees of freedom in the aerodynamic forces are con- 
stant for axial flow, but for nonaxial flow they are periodic functions of 
(The aerodynamic coefficients are defined in appendix A2.) 


2.4 Rotor Speed and Engine Dynamics 

The rotor rotational speed degree of freedom (^ s ) is frequently an impor- 
tant factor in rotorcraft dynamics. With a turboshaft engine, the rotor 
behaves almost as a windmill. For a powered wind-tunnel model with an elec- 
tric motor, the motor inertia and damping can significantly reduce the rotor 
speed perturbations. The equation of motion for \jj Q is given by the rotor 
torque equilibrium. We shall examine the extremes of a windmilling rotor and 
constant rotor speed and then derive a more general model including the engine 
inertia and damping. 

For windmilling or autorotation operation, the rotor is free to turn on 
the shaft. No torque moments are transmitted from the rotor to the shaft and 
no shaft rotational motion is transmitted to the rotor. The equation of motion 
for the rotor speed perturbation (ip s ) is just Q = 0, or yCq/aa - 0. In 
axial flow, there is no spring term in the ip s equation, so the system is 
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first order in t|> s . The rotor azimuth perturbation iJj s is defined with 
respect to the shaft axes, which also have a yaw angle a z ; thus the rotor 
speed perturbation with respect to space is i> s + a z . 

For constant rotor speed, the ip s degree of freedom and equation of 
motion are dropped from the system (i.e., the appropriate row and column are 
eliminated from the coefficient matrices) . The solution for the rotor speed 
perturbation is just \p s = 0, so the rotational speed with respect to the 
shaft axes is constant at Q. 

Now consider a more general case, including the inertia and damping of 
the engine or motor. Any drive train flexibility is neglected since it usually 
does not have a major role in the rotor dynamics. Thus the engine model used 
does not add degrees of freedom to the analysis, it simply includes the engine 
inertia and damping in the rotor torque equilibrium. The engine effects are, 
of course, amplified by the transmission gear ratio. The equation of motion 
for i p s is then 

Q = -Ve 2 *s - Ve 2 *s 


where Ig is the engine rotary inertia, rg is the transmission gear ratio, 
and Qft = dQg/dfig is the engine speed damping coefficient. Normalizing as 
usual, the equation becomes 


— ^ + r 2 I *ip + r 
aa ■ E E E 


2 1n*i = 0 


where Ig* = Ig/NIj> and Q^* = Q^/NI^fi. For the windmilling case, Ig* and 
Q^* are set to zero. For constant rotor speed, the ip s degree of freedom 
and equation are dropped. This model may also treat the engine-out case, for 
which there is no engine damping, by setting Q^* =0. 


The engine damping may be related to the engine trim operating condition 
by 



rotor 

c_ 2 n z / 

E rotor 


where k is a constant depending on the engine type. In coefficient form, 
then, 


* _ 


E 


xy 


fa 

oa 


where Cg is the trim rotor torque or power coefficient. This expression is 
applicable to a wide variety of engines (refs. 8 to 10). The constant takes 
the value k — 1 for a turboshaft engine (refs. 8 and 9) or for a series d.c. 
electric motor (ref. 10). It takes the value k s* 1/(1 - q) for an induction 
electric motor or an armature-controlled shunt d.c. motor (ref. 10; q is the 
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motor efficiency) . For a field- controlled shunt d.c. motor, the only damping 
is mechanical or the damping of the load, so k = 0 (ref. 10). For a synchro- 
nous electric motor, there is a spring on the rotational speed due to the 
motor, so the above model is not applicable (ref. 10). Generally, the inertia 
of the engine or motor is more a factor in the dynamics than the damping. 

A rotor speed governor may be included in the model. For example, inte- 
gral plus proportional feedback of the rotor speed perturbation to collective 
pitch (neglecting the governor dynamics) gives the control equation: 

A0 COn = K t ^ + K v \p 
o I r s P T s 

Note that the integral feedback (Kj ^ 0) adds a spring term to the rotor speed 
dynamics. A governor would be unusual for a wind-tunnel model unless an 
actual helicopter was being tested (in which case, a throttle governor would 
be more likely). Moreover, the governor has little effect on the rotor 
dynamics generally, because it is basically a very- low-frequency feedback 
system. The rotor speed governor has a more important role in helicopter 
flight dynamics. Thus a more detailed model is given in part II (section 9) 
for the governor, as well as for the engine and transmission dynamics. 


2.5 Inflow Dynamics 

The aerodynamic forces on the rotor result in wake-induced inflow veloc- 
ity at the disk, for both the trim and transient loadings. The wake- induced 
velocity perturbations can be a significant factor in the rotor aeroelastic 
behavior; an extreme case is the influence of the shed wake on rotor blade 
flutter (ref. 11). Therefore, the rotor inflow dynamics should be incorpo- 
rated into the aeroelastic analysis. However, the relationship between the 
inflow perturbations and the transient loading is likely more complex even 
than for the steady problem (nonuniform wake-induced inflow calculation) , and 
models for the perturbation inflow dynamics are still under development. In 
the present analysis, an elementary representation of the inflow dynamics is 
used. The basic assumption is that the rotor forces vary slowly enough (com- 
pared to the wake response) that the classical actuator disk results are 
applicable to the perturbation as well as the trim inflow velocities. 

A contribution to the velocity normal to the rotor disk of the following 
form is considered: 


Sil = X + X r cos iL» + X r sin \b 
r x m y m 

Here the perturbation inflow component X is uniform over the rotor disk, 
while the inflow due to A x and varies linearly over the disk. The inflow 

dynamics model must relate these inflow components to the transient aerody- 
namic forces on the rotor, specifically to the thrust C T , pitch moment C™ > 
and roll moment C™ . * 
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2.5.1 Moment-induced velocity perturbations- For hover, the perturbation 
inflow 6v(r,\p) at a point on the rotor disk may be related to the perturba- 
tion of the local disk loading dT/dA by 

^ _ dT/dA 
‘ 2pv 0 


where v c is here the trim value of the induced velocity. This result can be 
derived either by momentum theory for the disk element dA (cf. v Q = T/2pAv c 
for steady state) or by vortex theory (ref. 12). It is applicable only for 
harmonic changes of the blade loading, however, that is, variations occurring 
at a frequency w/ft = n/rev in the rotating frame, where n is a nonzero 
integer. 


Assuming a linear variation of the loading over the disk, the pitch and 
roll moments give 

, T 6M 6M 

dA = " 4 “RA r COS ^ + 4 IT r Sln ^ 

It follows that <$v also has a linear variation. Furthermore, the moments 
involve harmonic (1/rev) loading, so the above result is applicable, and it 
gives in coefficient form: 



r cos ip + 



r sin ip 


where A Q is the trim inflow ratio (see section 2.3.3). This result can be 
extended to forward flight following the usual approach of momentum theory. 
The mass flow through the differential disk area dA is determined by the 
resultant velocity through the disk, so generally m = p(V 2 + y q 2 ) 1 ^ 2 dA. 
Then dT = 2m 6v gives 


6v = 


dT/dA 


2p)/l 


V z + v_ 


Thus the inflow perturbation becomes 

26C, 




26C. 


M 


fiX - - 


A 2 + v 


r cos + 


x 


i ~ \ 


r sin ^ 


For speeds above transition (p above 0.1 or 0.15), this is approximately 

26C M 260, 

M M 

p V X 

<5A = “ r cos ^ H r sin ^ 

U T II 
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which may also be obtained directly from the differential form of the induced 
velocity in forward flight (X.^ = C T /2y). 

2.5.2 Thvust-induoed velocity perturbation- Now consider the inflow 
changes due to rotor thrust transients. The above relation between 6v and 
dT/dA is not applicable in this case. That relationship is based on low- 
frequency variations of a harmonic blade loading (ref. 12). The thrust 
changes correspond to low-frequency variations of the mean blade loading, and 
thus a different approach is required. 


For thrust perturbations, it is possible to simply consider a perturba- 
tion form of the hover momentum theory result for the trim inflow, 

X Q * (Cx/2) 1 / 2 . For low-frequency thrust changes then. 


6X 


3X _ 6C T 
ac T 6C t 4X C 


Note that the harmonic loading result would give an inflow change twice as 
large, 6X = 6Cx/2X 0 . The difference is due to the rotor shed wake. For har- 
monic loading (such as 1/rev variations due to the moments) , there is both 
shed and trailed vorticity in the wake, with half the inflow perturbation 
coming from the shed wake and half from the trailed wake. The rotor thrust 
produces trailed wake vorticity only (i.e. , tip vortices), and hence half the 
wake influence as the rotor hub moments. 

The extension to forward flight is based on the momentum theory result 
for the trim inflow: 


Then 


X Q = u tan a Rp + 


2-)/ P 2 + X 0 2 


6X 


_9X_ 

3C„ 


6C„ 


«<L = 


2|/u 2 + X Q 2 + C X 0 /( v 2 + X() 2) 



(The last approximation is valid for small inflow.) In summary then, the 
result obtained for the inflow perturbation due to the unsteady thrust and 
moment changes is 


6C„ 


26C 


M 


26C 


M 


6X = 


2 ( x o + / m 2 + x 0 2 ) -r* 


2 + X o 2 


r cos \jj + 


y v 2 + x 0 2 


r sin ip 
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2.5.3 Lift deficiency function- Wake effects in unsteady aerodynamic 
theory are often represented by a lift deficiency function. Consider the lift 
deficiency function implied by the above results. The aerodynamic pitch and 
roll moments on the disk may be written: 



where QS means the quasistatic loading, that is, all moment terms except 
those due to the wake. The induced velocity change due to the hub moments is 
given above as 



Substituting for the inflow changes produces 



where the lift deficiency function C is 



Thus all rotor aerodynamic hub moments are reduced by the factor C due to 
the rotor wake influence, which can significantly affect the dynamic behavior. 
For forward flight, the lift deficiency function is 


~ 1 + (oa/8y) 

and, for hover. 


C = 


1 

1 + (oa/8A 0 ) 
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The hover result is, in fact, the same as the low reduced frequency, harmonic 
loading limit of Loewy's lift deficiency function (see refs. 11 and 12). 

Miller shows that it is a good approximation to the real part of Loewy's func- 
tion to at least a reduced frequency of 0.5, although it neglects the phase 
shift entirely, of course (ref. 12). 


Similarly, the aerodynamic thrust changes can be written: 



and the inflow perturbation is 


2 (v + Vv 2 + *o 2 ) 


Hence 



with the lift deficiency function here: 


s(a 0 + y y 2 + x Q 2 ) 

For forward flight, this gives the same function as for the moments, while for 
hover the thrust changes give 


C = 


1 


1 + 


oa 

16\ 0 


Thus the wake effects reduce the aerodynamic thrust forces by the factor C; 
the reduction in hover is not as large as for the moments, however, because of 
the shed wake effects for moments. 

Typical values of the lift deficiency function from the above expressions 
are C = 0.8 for forward flight; C ^ 0.7 for thrust changes in hover; and 
C = 0.5 for moment changes in hover. In practice, it is more convenient to 
incorporate the inflow influence in the aeroelastic analysis using a differ- 
ential equation model rather than a lift deficiency function. The lift 
deficiency function is useful, however, in estimating the magnitude of the 
wake effects. 
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2.5.4 Inflow due to velocity perturbations- Perturbations of the rotor 
velocity will also produce a wake- induced velocity change because of the change 
of the mass flow through the disk. Consider again the momentum theory result 
for the trim induced velocity: 

1 2)/y 2 + (v z + 2 

Then, for low-frequency changes of the rotor inplane and normal velocity com- 
ponents (y and y z ) , the induced velocity perturbation is 

C (y6y + A 0 6y z ) 

6 A - 

2 (y 2 + A q 2 ) 3 / 2 + c t A q 


For hover, this reduces to 



and, for forward flight. 


<SA = — 5 y - — 6y 

2y 2 2y 3 

Including the thrust term, the total perturbation of the uniform induced 
velocity component is then 

C T 

6C T " (y 2 +T 0 2 T (v6y + W 
6A = ° 7 -■= 

2 ( A o + + A o 2 > 

The simplest means of incorporating the rotor velocity perturbation terms in 
the analysis is to treat them as additional terms in the thrust perturbation. 
To complete this expression, it is necessary to identify the velocity pertur- 
bations. Since the velocity components are in the shaft axes, the contribu- 
tions of the shaft motion and aerodynamic gusts are: < 5 y x = -x^ + ug, 

6y y = y h + v G , and 6y z = z h - w G . The shaft angular motion gives no net 
induced velocity change because the mass flow through the rotor disk depends 
only on the magnitude of the resultant velocity. 

2.5.5 Inflow dynamics model- Considering an inflow perturbation consist- 
ing of uniform and linear (1/rev) components, 

6 A = A + A r cos i{> + A r sin 
x y 
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the analysis above has related these terms to the perturbation aerodynamic 
thrust and hub moments on the rotor as: 



This relationship might be generalized to 




where 9X/3L may be a full nine-element matrix. Here we shall consider only 
the diagonal terms obtained from momentum and vortex theory. As a further 
generalization, a time lag in the inflow response to loading changes can be 
included by introducing a first-order time lag term: 


T i + 1 - |f t 

oL 

Following reference 12, let x = k( 9X/3L), where k is a diagonal matrix: 

Po 0 OH 






The values k q = 0.85 and k c “ K s = are from reference 13. Refer- 

ences 13 and 14 show that these values for the time lag correlate fairly well 
with experiment, and also that the lag does not, in fact, have a very important 
role in the dynamics. 


In summary then, the inflow dynamics model consists of a set of first- 
order, linear, differential equations for the inflow variables X, X x , and Xy*. 
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Expressions for the aerodynamic thrust and moment to complete these equations 
may be obtained above. It is also necessary to include the aerodynamic 
forces due to the inflow perturbations in the equations of motion for the 
rotor degrees of freedom. This is an elementary extension of the results of 
section 2.3, for, by comparing terms in 6 up, it is seen that A corresponds 
to z^, A x to -tty, and Ay to a x . Thus, generally, three degrees of free- 
dom and equations have been added to the system that describes the rotor 
dynamics. 

The time lag is not usually an important factor , so the quasistatic model 
for the inflow dynamics is generally sufficient. Dropping the time lag terms, 
the equations for A, A x , and A reduce to linear algebraic equations. Thus, 
in the quasistatic case, the inflow perturbations do not increase the order of 
the system. The wake influence reduces to an algebraic substitution relation, 
which, if incorporated analytically, would lead to lift deficiency functions 
as derived previously (with large-order systems, it is more practical to 
accomplish the substitution numerically) . 

An elementary model has been presented for the rotor inflow dynamics. 

Such a model has shown good correlation with experiment (refs. 13 and 14), and 
it gives the correct low-frequency limit for the lift deficiency function (cf. 
refs. 11 and 12). A more accurate model will probably be necessary for some 
applications, and a more complex model might be constructed, but further 
research in this subject is required before a model becomes available which is 
uniformly more valid than that presented here. 


2.6 Rotor Equations of Motion 

2.6.1 Inertia coefficients- At this point, the linear differential equa- 
tions of motion for the rotor model are constructed. The equations of motion 
are in the nonrotating frame, that is, the Fourier coordinate transformation 
has been applied to the bending and torsion degrees of freedom of the blade. 
For now, only a three-bladed rotor is considered (N = 3); the equations are 
extended to an arbitrary number of blades in section 4. Matrix notation is 
used for the equations. The coefficient matrices are constructed from the 
results of the previous sections (primarily, sections 2.2 and 2.3). The vec- 
tors of the rotor degrees of freedom (x R ) , shaft motion (a), rotor blade pitch 
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input (v R ) , aerodynamic gust input (g s , in shaft axes), and the hub forces and 
moments (F) are defined as: 



(k) C k. ^ 

Note that, in the rotor degrees of freedom x R , the notation g and 0 / 

is intended to cover as many bending and torsion modes as the analysis requires. 


66 



The equations of motion for the rotor, and the hub reactions take the 

form: 

A 2*r + A 1*R + *0^ + V + V * Bv r + “aero 

F ‘ C 2*R + Vr + C 0*R + V + ^1“ + F aero 

Here only the structural and inertial terms are included in the coefficient 
matrices; M aero and F aero are the aerodynamic forcing terms. (These inertial 
matrices are defined in appendix Bl.) 

2.6.2 Aerodynamic coefficients — axial flax- The aerodynamic terms F aer o 
and M aero of the rotor equations of motion and the hub reactions are required 
to complete the differential equations of the rotor model. They are obtained 
by summing over all N blades of the blade aerodynamic forces in the rotating 
frame (section 2.3) and introducing the Fourier coordinate transformation for 
the blade bending and torsion degrees of freedom as required. 

The case of a rotor operating in axial flow (y = 0) is considered first. 
In axial flow, the aerodynamic coefficients for the blade, forces in the rotat- 
ing frame are constants, independent of the blade azimuth angle ip m . The 
coefficients are also then entirely independent of the blade index (m) ; hence 
the summation over the N blades operates only on the system degrees of 
freedom, not on the aerodynamic coefficients themselves (which factor out of 
the summation) . The result for the required aerodynamic forces is 

‘“aero * *1*11 + Vr + A l“ + V ' Vs 

F aero = C l i R + C o x r + + V + Vs 

The coefficient matrices are constant for axial flow (see appendix B2) . 

2.6.3 Aerodynamic coefficients — nonaxial flow- Finally, the aerodynamic 
terms for the rotor operating in nonaxial flow, y > 0, are considered. In 
this case, the aerodynamic coefficients of the rotating blade forces are 
periodic functions of -i|/ m because of the periodically varying aerodynamics of 
the edgewise moving rotor. It follows that the rotor in nonaxial flight is 
described by a system of differential equations with periodic coefficients. 

It is possible to express the aerodynamic coefficients of the rotating 
blade forces as Fourier series, and then to obtain the coefficients of the 
nonrotating equations in terms of these harmonics. The result is rather com- 
plicated, however, and, in the present analysis, it would even be necessary to 
numerically evaluate the harmonics of the Fourier series. The simplest 
approach for numerical work with large-order systems is to leave the coeffi- 
cients of the nonrotating equations in terms of the summation over the N 
blades of the rotor. The summation is easily performed numerically, and it is 
found that this form is also appropriate for a constant coefficient approxima- 
tion to the system. 
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The required aerodynamic forces again take the form: 


= A 1*R + Vr + A l“ + V ~ B G g £ 


= C x x R + C Q x R + C x a + C Q a + D G g £ 


For nonaxial flow, the coefficient matrices are periodic functions of the 
blade azimuth angle = ij/ + mAty (A^ = 2ir/N) . For a three-bladed rotor con- 
sidered here, the period of the equations in the nonrotating frame is 
Aip = (2/3) it (120°). (The coefficient matrices are given in appendix B3.) 


2.7 Constant Coefficient Approximation 

The rotor dynamics in nonaxial flow are described by a set of linear dif- 
ferential equations with periodic coefficients. A constant coefficient 
approximation for nonaxial flow is desirable (if it is demonstrated to be 
accurate enough) because the calculation required to analyze the dynamic 
behavior is reduced considerably compared to that for the periodic coefficient 
equations, and because the powerful techniques for analyzing time- invariant 
linear differential equations are then applicable. However, such a model is 
only an approximation to the correct aeroelastic behavior. The accuracy of 
the approximation must be determined by comparison with the correct periodic 
coefficient solutions. The constant coefficient approximation derived here 
uses the mean values of the periodic coefficients of the differential equa- 
tions in the nonrotating frame. 

To find the mean value of the coefficients, the operator 


_L 

2n 



is applied to the periodic aerodynamic coefficients (given in section 2.6.3 
and appendix B3) , resulting in terms of the form: 
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Here M nc and M ns are the harmonics of a Fourier series representation of 
the rotating blade aerodynamic coefficient M: 


M(^ m ) = M° + 2 MnC cos n ^m + mI1S sin n ^m 
n=l 

In the present case, these harmonics must be evaluated numerically. The aero- 
dynamic coefficient M is calculated at J points, equally spaced around the 
azimuth. Then the harmonics are calculated using the Fourier interpolation 
formulas: 

M° - 4 E M(* ) 

J 3 3 

M 110 = 4 2 

j 3 3 

l^ S =42 M(ip.)sin ni|> 

j 3 3 

where tyi « jAip - j(2tt/J) (j = 1, ..., J). The number of harmonics required 
is n = N-l for N odd and n « N-2 for N even (N is the number of 
blades). Good accuracy from the Fourier interpolation requires at least that 
J = 6n. Using these Fourier interpolation expressions, the required harmonics 
are 



It follows then that the constant coefficient approximation is obtained 
from the periodic coefficient expressions by the simple transformation: 

j2 (••) M <*m> * J £ (•••> HOM 

m=l j=l 
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The summation over N blades (m = 1, . N; Ai|» = 2ir/N) for a periodic coef- 
ficient is replaced by a summation over the rotor azimuth (j =1, . .., J; 

Aip = 2tt/J) for the constant coefficient approximation. This is quite conven- 
ient since the same procedure may be used to evaluate the coefficients for the 
two cases, with simply a change in the azimuth increment. The periodic coef- 
ficients must be evaluated throughout the period of ip = 0 to 2ir/N, of 
course; the constant coefficient approximation (mean values only) is evaluated 
only once. 

With the substitution (1/N) 2 (1/J) 2) , the results given in appen- 

® J 

dix B3 for the periodic coefficient matrices are directly applicable to the 
constant coefficient approximation as well. 


3. ADDITIONAL DETAILS OF ROTOR MODEL 


3.1 Rotor Orientation 

The rotor orientation is defined by a rotation matrix between the shaft 
axes and a tunnel axis system: 



The wind-tunnel axis system used has the x axis directed downstream, the 
z axis positive upward, and the y axis positive to the right (looking 
upstream) . The wind-tunnel velocity is then Viij and the components of the 
velocity in the shaft axes are 



The hub plane angle of attack and yaw angles may be obtained from the follow- 
ing expressions: 

a HP = tan_1 y z/V y x 2 + V 

*HP = tan_1 V y * 

As examples, we shall consider a helicopter, a propeller test rig, and a tilt- 
ing proprotor aircraft in a wind tunnel. 
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3.1.1 Helicopter- For a helicopter rotor in a wind tunnel with a turn- 
table, the shaft axes orientation is given by first yaw to the right by 
then pitch aft by 0^.. Thus the rotation matrix is 


R 


ST 


c.c. 

6 ip 


L s eS 


-c 0 s , 

0 ij; 


~S n S . 
0 ip 


~ s e 

o 

c. 


3.1.2 Propeller test rig- Consider a propeller test rig that may be yawed 
by 'I't (positive to the right), with axial flow for = 0; Then, 


ST 


-S 

0 

-C 


<P 


* 


-C 

0 

s 


<p 


* 


3.1.3 Tilting proprotor aircraft- The nacelle and rotor of a tilting 
proprotor aircraft can be tilted by an angle oip, where ap = 0 for axial 
flow (airplane configuration) and ap = 90° for edgewise flow (helicopter 
configuration) . It is assumed that the nacelle has a cant angle <(>p (posi- 
tive outward in helicopter mode, zero for airplane mode) , and that the aircraft 
has a pitch angle 0^ (positive nose-up) . The resulting rotation matrix is 


ST 


C.S C Q + C S fl 
<p a 0 a 0 

C.S.(1 - C )C a + S A S s. 
cf> <f> a 0 <paQ 

-(cjc + sJ)C Q + css 


<f> a <1> 7 e 


$ a e 


-s.s 

4> a 

c 2 + S. 2 C 
<f> $ a 

-C.S, (1 - c ) 
4> <f> a 


-C.S s. + c C Q 
<p a 0 a 0 

-C.S. (I - C )S. + s.s c 0 
<p 4> a ,0 <j> a 0 

(C^C + S 2 )S„ + C.S c 


(j) a 


C(). 0 (j> u 0J 


'3*1.4 Gust orientation - The aerodynamic gust components are defined with 
respect to the tunnel axes (wind axes) for the analysis of the rotor and wind- 
tunnel support dynamics. The vector of gust velocities seen by the rotor is 
then 



where uq is positive downstream, vq positive from the right, and wq posi- 
tive upward. The rotor aerodynamic forces were derived considering gust com- 
ponents in the shaft axes, g g . The substitution g s = Rg^g into the equa- 
tions of motion then transforms the gust components to the tunnel axes. . 
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3.2 Rotor Trim 


There are two direct requirements in the dynamics analysis for the trim 
solution for the rotor^blade^motion and rotor performance: first, the trim 

bending deflection (x Q i 4* z Q k) is required for the coefficients; second, the 
evaluation of the aerodynamic coefficients requires the lift and drag loading 
of the rotor blade. The result of this role of the trim solution in the 
analysis is that the aeroelastic behavior depends generally on the operating 
state of the rotor. The evaluation of the coefficients of the equations of 
motion must therefore be preceded by a preliminary calculation of the rotor 
trim. The trim solution for the blade motion is periodic in the rotating 
frame for the general case of nonaxial flow; for the axial flow (y = 0) , the 
blade motion is steady in the rotating frame. For the trim blade motion in 
the present analysis, only the bending and gimbal degrees of freedom are con- 
sidered. It is assumed that there is no unsteady shaft motion, that the rotor 
speed is constant, and that there is no torsion/pitch motion (except cyclic 
control and any bending/torsion coupling). 

In the trim solution, it is assumed that all blades have the same 
periodic motion and that the gimbal deflection is constant. The equations of 
motion are solved by a harmonic analysis method, which solves the rotating 
equations of motion directly for the harmonics of a Fourier series expansion 
of the periodic motion. The equations for the blade motion are obtained from 
the above analysis, for the bending and gimbal degrees of freedom: 



The inertia coefficients are defined in appendix Al, and the aerodynamic 
forces are evaluated using the trim velocity components (section 2.3). 


After the blade motion calculation converges, the rotor performance is 
evaluated, including the mean aerodynamic forces and moments the rotor pro- 
duces at the hub (in particular, the rotor thrust and power). In an outer 
loop, there can be an iteration on the controls to trim the rotor to some 
desired operating state. Possible options include adjusting the collective to 
achieve a target thrust or torque, or adjusting collective and cyclic pitch to 
achieve a target thrust and flapping or a target lift and propulsive force. 


For axial flow,(y = 0), the trim solution is independent of ip (assuming 
no cyclic pitch input). The gimbal motion is zero, and the equation for the 
blade bending modal deflection reduces to 


q'k "k 


^ = 1 


V 


+ Y 


.1 /F F \ 

f • (— ±u - — k )dr 
k \ac B ac B \ 
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So, in this case, only the mean blade bending motion is nonzero, and an itera- 
tive solution for the blade motion is not required. Furthermore, of the 
forces and moments on the rotor hub, only the thrust and torque are nonzero. 


3.3 Lateral Velocity 

For the development of the aerodynamic analysis of the rotor in sec- 
tion 2.3, it was assumed that the trim velocity of the rotor was in the x-z 
plane of the shaft axis system. Generally, however, it is also necessary to 
consider a lateral velocity component, that is, by - Js * V. An alternative 
would be to rotate the shaft axes until jg • ^ = 6, but that would imply a 
redefinition of the rotor zero azimuth position for every flight state. Such 
a redefinition of ip is not desirable since it changes the values of param- 
eters such as control— system phasing, and even the definition of the rotor 
degrees of freedom such as tip-path-plane tilt. Hence it is preferable to 
directly incorporate the effects of the lateral velocity in the analysis. The 
velocity of the air seen by the rotor disk has then three components in the 
shaft axes^ y x , positive aft; y y _^ positive from the right; and y z , positive 
downward (V air = y x ig - y y j s - y z Ic s ) . 

The incorporation of y y into the analysis developed in section 2 is 
straightforward. The only influence is on the rotor aerodynamics. In the 
trim induced velocity, y 2 is replaced by y x 2 + y y 2 (section 2.3.3); this 
substitution is also made in the coefficients of the equations for the inflow 
perturbations (section 2.5.4). In the trim velocity of the blade (sec- 
tion 2.3.3), the quantity y cos i|> m is replaced by y x cos ij> m - y y sin <J> m , 
and the quantity y sin <|/ m by y x sin <J; m + y y cos 'J'm* these substitutions are 
also required in the C, 6, and qi aerodynamic derivatives and in the unsteady 
terms of the torsion equation aerodynamic coefficients (appendix A2) . Finally, 
the quantity ya y is replaced by y x a y + y y « x in the perturbation velocity 
6up (section 2.3.3). It follows that lateral velocity y y terms are added to 
the a x columns (fourth columns) of the aerodynamic matrices A 0 and Cq in 
the equations of motion, in a fashion similar to the y x terms in the a y 
columns (fifth columns; see appendices B2 and B3) . 


3.4 Clockwise Rotating Rotor 

The equations of motion for the rotor have been developed assuming coun- 
terclockwise rotation of the blades as viewed from above. The equations for a 
clockwise rotating rotor are obtained by implementing a mirror-image transfor- 
mation consisting of the following sign changes: 

(a) Change the signs of the y h , a x , and a z columns of the A and C 
matrices. 

(b) Change the signs of the Cy, Cq, and C M rows of the C, C, and D G 

matrices. x 

(c) Change the sign of the Vq column of the Bq and Dq matrices. 
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(d) Change the sign of the lateral velocity y . 

The degrees of freedom of the rotor remain defined with respect to the actual 
rotation direction; for example, blade lag motion is still positive opposing 
the blade rotation, and the lateral tip-path-plane or gimbal tilt $is is 
still positive toward the retreating side of the disk. The mirror-image trans 
formation accounts for the reversals of some components of a, F, and g (which 
are defined in the nonrotating system) relative to the counterclockwise or 
clockwise rotating rotors. 


3.5 Blade Bending and Torsion Modes 


3.5.1 Coupled bending modes of a rotating blade- Equilibrium of the 
elastic, inertial, and centrifugal bending moments on the blade gives the 
differential equation for the coupled flap/lag bending of the rotating blade 
(see section 2.2). For free vibration— the homogeneous equation with har- 
monic motion at the natural frequency v — we obtain the modal equation for 
bending of the blade: 


(Ein")" 




mv^n 


0 


Here n (r ^ z Q i - ^x Q k is the bending deflection (mode^shape) , 

El = EI zz ii -I- EI xx kk is the bending stiffness dyadic, Q = ftk^ is the rotor 

rotational speedy v is the natural frequency of the mode, and 

K g = Kpigig + K^kgkg is the hinge spring dyadic. The boundary conditions are 

(a) At the tip (r = R) : Elrj" - (EIn") T = 0. 


(b) At the root (r = e) : r\ = n T = 0 for a cantilever blade; n = 0 

and Eln" = K s n f for an articulated blade. 

The root boundary condition is applied at the offset r = e to allow for 
hinge offset of an articulated rotor or a very stiff hub of a hingeless rotor. 
The hinge springs are assumed to be oriented in the hub frame here, but could 
also include a component that rotates with the blade pitch; Kp is the flap 
spring and K L , the lag spring constant. 

This is an eigenvalue problem, a differential equation in r for the 
mode shapes n and the natural frequencies v. The equation with the appro- 
priate boundary conditions constitutes a proper Sturm-Liouville eigenvalue 
problem. It follows that the solution exists — a series of modes r)i(r) and 
the corresponding natural frequencies — and that the modes are orthogonal 
with weight m. Hence, if i ^ k, then 

n . • Th m dr = 0 
x k 
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The frequencies satisfy the energy balance relation: 



The modal equation is solved by a Galerkin method. The mode shape is 
expanded as a finite series in the functions f^Cr): 

n = 2 c^?^(r) 

We require that each satisfy the boundary conditions on n> then the sum 

automatically does. Since a finite series is required for computation, this 
is an approximate calculation. For best numerical accuracy, the functions 
should then be chosen so that at least the lower-frequency modes can be well 
represented. Substituting this series in the differential equation and 
operating with 

l V<- • -> dr 

reduces the problem (after integration by parts and an application of the 
boundary conditions) to a set of algebraic equations for c = [c^] : 

(A - v 2 B)c =0 


The coefficient matrices are 


\i = *k< e) + 

\i - / dr 

*e 








dr 


The eigenvalues of the matrix (IT^A) are the natural frequencies v 2 of the 
coupled bending vibration of the blade, and the corresponding eigenvectors c 
give the mode shape n- As a final step, the modes are normalized to unity at 
the tip: |n(l) 1=1. 

A convenient set of functions for f^ are the following polynomials 
(ref. 7): 

f _ (n + 2) (n 4- 3) n+1 _ n(n + 3) n+2 n(n + 1) n+3 
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where x = (r - e)/(l - e) . For a hinged blade, fj = x is used. The poly- 
nomials satisfy the required boundary conditions but are not orthogonal func- 
tions. For the hinge spring term in (articulated blades only), note that 

f^(e) = 0, except for the first mode where fj(e) = 1/(1 - e). 

3.5.2 Articulated blade modes- For an articulated rotor blade, the modal 
differential equation need not be solved if the higher bending modes are not 
required. Rigid lag and flap motion about the hinges gives the two lowest 
frequency modes : 



Note that separate hinge offsets may be used for flap and lag motion. The 
natural frequencies are obtained directly from the energy balance relation: 



v 


2 

flap 



3.5.3 Torsion modes of a nonrotating blade- Equilibrium of the elastic 
and inertial torsion moments (see section 2.2) gives the modal equation: 

(GJC')’ + I 0 w 2 5 = 0 

with the boundary conditions £' = 0 at the tip (r = R) , and 5 = 0 at the 
pitch bearing (r = r-^) . The modes are orthogonal with weight Iq. Hence, 
if i / k, then 


/ R 


Ve = 0 


’FA 
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GJ5' 2 dr 


The frequencies satisfy the relation: 

R 

FA 
R 

I 0 5 2 dr 
FA 

These are nonrotating torsion modes, so the solution is independent of 
the rotor speed or collective pitch. The equation is solved by a Galerkin 
method. Writing 5=2 c^f^r), where the functions f^_ satisfy the boundary 
conditions on 5, and operating with f 1 f v (...)dr on the differential equa- 

tion produces a set of algebraic equations for c = [c ± ] : 

(A - m 2 B)c = 0 

The coefficient matrices are 

\i ■ jf f i £ i dr 

FA 

\i - /' Vk f i dr 

r 

FA 

The eigenvalues of the matrix (B -1 A) give the natural frequencies of the tor- 
sion vibration, and the corresponding eigenvectors for c give the modes. 
Finally, the torsion modes are normalized to unit at the tip, 5(1) = 1. 

A convenient set of functions to use for is the solution for the 

torsion modes of a uniform beam: 



These functions satisfy the boundary conditions and will usually be close to 
the actual mode shapes. 



3 . 6 Lag Damper 

For articulated rotors, the mechanical lag damper has an important role 
in the dynamics. A lag damping term is added to the blade bending equation 
of motion (section 2.2): 
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s\ q k 


2 , 

'k V 


v, ‘q, ) + 2 


q k q i 1 


q i + i §lag k B * q k (e) ^B * q l (e > q i + 



M 

q k aer0 

+ I* 

ac q. o 

k 


where gl a g “ C^/lfoft and is the lag damping coefficient. The quantity 

ks * n^(e) is the slope of the kth bending mode in the lagwise direction, 
just outboard of the lag hinge. The manner in which the lag damping enters 
the equation of motion is obtained by a Galerkin or Rayleigh-Ritz analysis. 

The lag damper results in a bending moment at the lag hinge. Thus it is 
necessary to evaluate moments at the blade root by integrating along the span, 
which has, in fact, been our practice. 


3.7 Pitch/Bending Coupling 


The kinematic pitch/bending coupling Kp and the pitch/ gimbal coupling 

i 

Kp have a significant role in the rotor dynamic behavior. The definition of 
G 

Kp # is the rigid pitch motion due to a unit deflection of the ith bending 

mode: K p = - d0/dq-. For an articulated rotor, the first "bending" modes 

£ i 


are rigid lag and flap motion about the hinges. The pitch/ flap coupling is 
often defined in terms of the delta- three angle: Kp = tan <$ 3 . It is possible 
to simply input the kinematic coupling parameters to the stability calcula- 
tions if values are available from either measurements or some other analysis. 
It is also desirable, however, to be able to calculate the coupling from a 
model of the blade root geometry. 


Figure 12 is a schematic of the blade root and control-system geometry 
considered here, which shows the position of the feather bearing, pitch horn, 
and pitch link for no bending deflection of the blade. The radial locations 
of the feather bearing and pitch link are rp A and r pH , respectively; the 
lengths of the pitch horn and pitch line are x pH and x pL . The orientation of 
the pitch horn and pitch link are given by the angles tj> P g + 9 75 and <f>PL* 
Control input produces a vertical motion of the bottom of the pitch link and 
hence a feathering motion of the blade about the pitch axis. The bending 
motion of the blade, with either bending flexibility or an actual hinge 
inboard of the pitch bearing, produces an inplane or out-of-plane deflection 
of the pitch bearing. With the bottom of the pitch link fixed in space, a 
pitch change of the blade results. The vertical and inplane displacements of 
the pitch horn (the end at r pH ) due to bending of the blade in the ith mode 
are 

Az = q 1 i B * " n i ( ' r FA^ r FA “ r PH^ 


Ax 


= " q i k B 


n i (r FA ) (r FA 


- r pH ) ] 
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The kinematic pitch/bending coupling is derived from the geometric constraint 
that the lengths of the pitch horn and pitch link be fixed. The result is 


K_ 


(cos <l>p L i B + sm 9 pL k B ) * " n i^ r FA^ ^ r FA r PH^ 


-x pH cos(<|) pH + e 75 + * pL ) 


Similarly, for a gimballed rotor, the pitch/gimbal coupling is 

V 


- (r pH / Xp H )cos <fr pL 

cos(<j) pH + e 75 + 4> pl ) 


(\\± 


pitch horn horizontal 
cos(<(. pH + 0 75 + <j>p L ) 


3.8 Normalization Parameters 

It has been the practice here to deal with dimensionless quantities based 
on the air density, rotor speed, and rotor radius (p , ft, and R) . In addition, 
the equations of motion inertia coefficients, and aerodynamic coefficients 
have been normalized using the following parameters: I^, a characteristic 

moment of inertia of the blade; c m , the blade mean chord; and a, the blade 
two-dimensional lift-curve slope. These parameters have no effect on the 
numerical problem* It is essential, however, to be consistent in the defini- 
tion of the normalization quantities and the parameters that depend on them. 

In particular, the blade Lock number and rotor solidity are defined as 

pac R 4 
m 


The Lock number may be defined for standard air conditions (p Q ) , and then 
y(p/p 0 ) used in the analysis to account for the actual temperature and alti- 
tude at which the aircraft or wind tunnel is operating. It is convenient to 
use the rotor solidity as the primary parameter and to obtain the mean chord 
from c m /R = ott/N. 

4. ARBITRARY NUMBER OF BLADES 


4.1 Four or More Blades 

In this section, the rotor model is extended to an arbitrary number of 
blades. The equations derived in section 2 are for a three-bladed rotor. We 
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begin with a consideration of the case of four or more blades. Each rotating 
degree of freedom of the blade (i.e., bending or torsion motion) must result 
in N degrees of freedom for the rotor as a whole. Thus increasing the num- 
ber of blades adds degrees of freedom and equations of motion to the rotor 
description. In axial flow, these additional degrees of freedom do not couple 
with the collective and cyclic (0,1C,1S) degrees of freedom of the rotor. 

Hence the equations in section 2.6 remain valid for rotors with N > 3 also, 
and we need be concerned here only with the equations of motion for the addi- 
tional degrees of freedom. For nonaxial flow, however, all rotor degrees of 
freedom are coupled. 

The Fourier coordinate transformation for four or more blades adds the 
following degrees of freedom to the collective and cyclic variables for N =3: 



- E 

N ^ 

m 

q(m) 

cos 2ip m 

6 2S 

- E 

N 

m 

q(n0 

sin 2ip m 


^nc 

= - E 

N ^ 

q (m) 

cos mj;. 


m 



6 ns 

N ^ 

q (m) 

sin nif^ 


m 



N/2 

N ^ 

m 

q 0n) 

(-D m 


So 


q(m) = e o + £<e nc COS wp m + 3 ns sin nip m ) + 3 N/2 (-l) m 
n 

The result for the torsion/pitch modes is similar. The summation over m 
goes from 1 to N; the summation over n goes from 1 to (N - l)/2 for N 
odd and to (N - 2)/2 for N even. The 3^/2 an< * %/2 degrees of freedom 
are present only if N is even. Then there are a total of N nonrotating 
degrees of freedom corresponding to each bending and torsion mode of the blade. 

These additional degrees of freedom are not coupled inertially with the 
shaft or gimbal motion. The bending and torsion equations in the rotating 
frame are thus reduced to 
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The equations of motion in the nonrotating frame are obtained by application 
of the following summation operators: 

| £(. . .)cos , | £(. . Osin nip m , ± £(. . . ) (-l) m 

m m m 


and introduction of the Fourier coordinate transformation as required. 


The additional equations of motion for the rotor with four or more 
blades are then: 


Vr + a i*r + a o x r 



+ M 

aero 


The vectors of the rotor degrees of freedom (x R ) and blade pitch control 
inputs (v R ) are: 


*R 


3 


6 


00 

nc 

5 (k) 

ns 

, (k) 
5 N/2 

5 (k) 

nc 

(k) 

ns 

(k) 

N/2 


0 c °n 

nc 


v. 


R 


0 c °n 

ns 


e 


con 

N/2 


(The matrices of the inertia coefficients are given in appendix B4.) 


To complete the equations of motion, the aerodynamic terms must be eval- 
uated. The aerodynamic forces in axial flow still do not couple the addi- 
tional degrees of freedom for N > 4 with the shaft or gimbal motion. Thus 
the aerodynamic forces for the bending and torsion modes in the rotating frame 
reduce to 
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(The matrices of the aerodynamic coefficients are given in appendix B5.) 

The aerodynamic forces in nonaxial flow (y > 0) couple all degrees of 
freedom of the rotor with each other and with the shaft and gimbal motion. 
Then, not only are additional degrees of freedom and equations of motion 
involved if N > 3, but the number of blades also influences the equations and 
the hub reactions given in section 2.6. Rather than directly presenting the 
aerodynamic matrices for the general case of three or more blades in nonaxial 
flow, the analysis is extended by means of an observed pattern in the coeffi- 
cients. In the nonaxial flow equations in appendix B3, note the repeated 
occurrence of the following submatrices: 




(using the notation S n = sin n^ m and C n = cos n^ m ) . These matrices are a 
direct result of the introduction of the Fourier coordinate transformation 
(columns) and the application of the summation operators to obtain the non- 
rotating equations (rows) . The matrix DP arises from application of the 
Fourier transformation to the time derivatives (q^ or p^) . In the Bq and A 
matrices, only some columns of P and DP appear, while in the C matrices 
only some rows appear. The extension to an arbitrary number of blades (N > 3) 
is then, simply, 
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The constant coefficient approximation for the aerodynamic terms in non- 
axial flow is derived for N > 4 following section 2.7. It is found that the 
approximation is obtained from the periodic coefficient result by the 
transformation: 
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So the periodic coefficient results are still applicable to the constant coef- 
ficient approximation if the summation over the N blades is replaced by a 
summation around the rotor azimuth. 



Two-Bladed Rotor 


Rotors with three or more blades may be analyzed within the same general 
framework, but the two-bladed rotor is a special case. The rotor with N > 3 
has axisymmetric inertia and structural properties and hence the nonrotating 
equations have constant coefficients in axial flow. In contrast, the lack of 
axisymmetry with a two-bladed rotor leads to periodic coefficient differential 
equations, even in the inertial terms and in axial flow. Only in special 
cases (e.g. , shaft fixed, or with an isotropic support — analyzed in the 
rotating frame) are the dynamics of a two-bladed rotor described by constant 
coefficient equations. 

The rotor degrees of freedom in the nonrotating frame are obtained as 
usual from the Fourier coordinate transformation. For N = 2, the bending 
degrees of freedom are coning and teetering type modes: 







] 


m 

The pitch/torsion degrees of freedom 0g and are similarly defined. The 
teetering degree of freedom 0^, which corresponds to the gimbal motion of the 
rotor with three or more blades, is also included for the two-bladed rotor. 

The teetering motion is defined in the rotating frame; hence (see fig. 9), 


3 g = e T (-i) m 



The special characteristics of the two-bladed rotor dynamics are reflected in 
the appearance of the teetering-type degrees of freedom (3j, Oj, and 3x) > 
which take the place of the cyclic (1C and IS) motions for N > 3. 

The bending and torsion equations in the rotating frame, derived in sec- 
tion 2, remain valid for N = 2. The rotor equations of motion in the non- 
rotating frame are obtained by operating on the bending and torsion equations 
with the following summation operators: 1 

N V) (* •■•) an d ^ ^ • • -)(-l) 

m m 
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The equation of motion for the teetering degree of freedom is obtained from 
equilibrium of moments about the teeter hinge (in the rotating frame): 

-2M.J, + C t 3 t + K t 3 t = 0 

The teetering moment Mj- is given by the root flapwise bending moment: 

vJSi, 


a (n) (-»" 


m 


where C T and are the damper and spring constants about the teeter hinge 
in the rotating frame. In terms of the natural frequency and damping coeffi- 
cient, C T = 2I 0 ftCx and % = 2I 0 & 2 (vt 2 - 1). The hub forces and moments are 
obtained by summing the root forces and moments from both blades, as for 
N > 3. The equations are normalized and the inertia coefficients named in a 
manner similar to the N = 3 case (section 2.2). Inflow dynamics and the 
rotor speed dynamics are included as for N > 3 (sections 2.4 and 2.5). 

The vectors of the rotor degrees of freedom (xr) and the rotor blade 
pitch input (vr) are defined as follows for the two-bladed rotor: 


X R = 


6 


(k) 


0 


0 


>(k) 

S 1 

(k) 

0 

(k) 

1 

6. 


A 

A. 


v„ = 


con 

0 

con 


The vectors of the shaft motion (a) and the hub reactions (F) are defined as 
in section 2.6. The equations of motion then take the form: 

A 9 iL + A x + A^x,, + A„oi + A, a = Bv D + M 
2 R l R OR 2 1 R aero 


F - C 2 * R + CjX R + C oXr + C 2 a + CjO + 
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(The matrices of the inertia coefficients are given in appendix B6.) Note 
that there are periodic coefficients in the matrices coupling the rotor and 
shaft motion (A, C, £) . 


The aerodynamic forces M aero and F aero are required to complete the 
equations. The teetering aerodynamic moment is defined as 




aero 

ac 



r dr(-l) m 


The aerodynamic hub forces and moments are defined as for N > 3. The aero- 
dynamic forces on the bending modes are 


M. 


3 0 aero 


N L-i 


M 


— q, aero 
m k 


M Q 

3 x aero 




— q t aero 
iu 


(-1) 


m 


The definitions of the torsion aerodynamic forces Mq aero and Mq aero are 
similar. The aerodynamic forces are then: 0 1 


-M = 4- 4- A, a + A a - B_g 

aero 1 R 0 JR 1 0 G & s 


F = + C x n + G,a 4- C a 4* D 0 g 

aero 1 R OR 1 0 G a s 

The vector of aerodynamic gust components (g s ) is defined as in section 2.6. 
(The matrices of the aerodynamic coefficients are given in appendix B7.) 


A constant coefficient approximation for the two-bladed rotor may be 
obtained in a manner similar to the N > 3 case. For the aerodynamic forces, 
as usual the summation over N blades is simply replaced by an average over 
J points around the azimuth: 



1 


1 

cos 

^ j 

COS 

sin 

3 

sin ij/j 

(-l) m 


0 


M(^) 


Note that all aerodynamic terms involving (-l) m drop in the constant coeffi- 
cient approximation. For the inertia terms, the mean values of the periodic 
coefficients are easily obtained. 

The constant coefficient approximation is not as useful — or as accu- 
rate — for the two-bladed rotor as for N > 3. With three or more blades, the 
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source of the periodic coefficients is nonaxial flow, hence the periodicity ih 
of order u or smaller. At low advance ratio then, the constant coefficient 
approximation may be expected to be a good representation of the correct 
dynamics. The two-bladed rotor has also periodic coefficients due to the 
inherent lack of axisymmetry of the rotor. This periodicity is large even for 
axial flow, and neglecting it in the constant coefficient approximation is 
often a poor representation of the dynamics. 


4.3 Single-Bladed Rotor 


A single-blade analysis is useful for problems not involving the shaft 
motion or other excitation from the nonrotating frame. The only rotor blade 
degrees of freedom involved are the bending and torsion motion. The shaft 
motion, gimbal motion, and the rotor speed perturbation are dropped from the 
system. The hub reactions need not be considered. The single-blade analysis 
is, of course, in the rotating frame. The equations of motion for the bending 
and torsion modes of the blade are then: 


I* (q, + g v, q. + V, 2 q, ) + 2 V I* . q. - T. S* •• p. - £ S* „ P* 
q k VH k & s k H k k H k' ^ q k q i 1 " q k p i q k P i 

- Y/yM . q. - Y yM q. - YyM p. = 
" 1 q n q . i T q. q. l 1 q, p. 1 

4 k M i M k M i n k*i 

+ W k - + 2 foil + 2 W* - 2 - 2 'W* 


2 Y \^ " 2 V* 1 ' 2 1 v / i ' 2 V 1 ' V° 5k<rpA 


:),e 


con 


5. WIND-TUNNEL SUPPORT MODEL 


Consider now the aeroelastic model for the wind-tunnel support system. 

The rotor support is described by a set of linear constant coefficient differ- 
ential equations excited by forces and moments at the rotor hub. The hub 
motion produced by the support degrees of freedom completes the description. 

Let x g be the vector of the support degrees of freedom and v s , the vector 
of input or control variables for the support system. As for the rotor model, 
the vectors of the shaft motion at the hub (a) , the aerodynamic gust components 
(g) , and the rotor forces and moments acting on the hub (F) are defined as: 
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The gust components are here in the tunnel axis system (x downstream, y to 
the right , and z upward) , while the components of a and F are in the shaft 
axis system. The equations of motion for the wind-tunnel support then take 
the following form: 

a 2*s + a l*s + a 0 JC s = bv s + b G g + aF 
and the rotor hub motion is given by the linear transformation: 

a = cx s 

The equations are dimensionless, based on p, ft, and R. With F in rotor 
coefficient form, it is also convenient to normalize the equations by dividing 
by the characteristic inertia (N/2)I b . Note )Lhat the matrix a may always be 
obtained from the matrix c (reciprocity theorem) . 

The sophistication required of the description of the wind-tunnel strut 
and balance system and the aircraft or rotor test module depends on the dynamic 
problem being studied and also, of course, on the available data from which 
the model is to be constructed. Consider, for example, a general model based 
on a normal vibration mode description of the elastic wind-tunnel support 
system. The displacement u(r,t) and rotation ^(r,t) at an arbitrary point 
r are expanded in series of orthogonal vibration modes, with generalized 

coordinates q„, (t) : 
s k 

u(r,t) = 5^ q g (t)f (r) 
k k 

$(r,t) = 5^ q (t)y.(r) 

k s k ' 
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The differential equations for the degrees of freedom q s are then 

k 


"4 


S k + Ws k + \ \ 


)-% 


where is the modal mass and the natural frequency; g s is here the 

structural damping coefficient for the mode and is the generalized force. 


The hub motion is obtained from the mode shapes £ and y at the rotor 
hub, a = c Jq Sk j , where 


C = 


i £ 

■+ 

j S 


-t 

h 

k„ 


k 

K 

K 


• • • 


R ST^k 

R ST Y k 


• • • 


Here £ and y are in the tunnel axis system, so Rg T is the rotation matrix 
to the shaft axes (see section 3.1). The generalized forces due to the rotor 
hub forces and moments are {Q^} - aF, where 








Generally, Q^. may also have additional contributions, including mechanical or 
aerodynamic damping forces of the form 2 C^qg / , support-system control 


inputs of form 2 kki v Sj_> and aerodynamic gust forces on the support of the 
form 2 Making the equations dimensionless and normalizing as appro- 

priate produces the required model for the wind-tunnel support. 


6. COUPLED ROTOR AND SUPPORT MODEL 


The equations of motion have been derived for the rotor and for the wind- 
tunnel support. Now these equations may be combined to construct the set of 
linear differential equations which describes the dynamics of the coupled 
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system. The vectors of the degrees of freedom (x) , the control inputs (v) , 
and the aerodynamic gust components (g) are defined as: 


— ' — 


— — 


— — 

X R 


V R 


U G 


II 

> 


> 8 = 

V G 

x s 


v s 







G 

“ — 


— — 


mm mm 


for the coupled rotor and wind-tunnel support model 
form: 

A^x + A 1 x + AqX « Bv + B^g 

To derive the coupled equations, recall the results for the rotor equa- 
tions of motion and the hub forces and moments from section 2.6: 

A 2 X R + A 1 X R + A 0 X R + V + V + V = BV R + B G 8 S 

F = C 2 X R + C 1*R + C 0 X R + C 2 5 + C l“ + V + Vs 
and the results for the support equations and the hub motion (from section 5) 

a 2 X s + Vs + a 0 X s = bV s + b G S + aF 


The equations of motion 
then take the following 


a - cx s 


The gust components in the rotor equations are transformed to the tunnel 
(wind) axes by the substitution g s = Rgx&* The coupled equations of motion 
are obtained by substituting the hub motion (a) into the rotor forces and 
moments (F) and then substituting for F into the support equations of 
motion. The following coefficient matrices for the complete system may then 
be constructed: 


A^c 


-aC r 


a 2 “ aC 2 C 


A 1 " 


A l c 


-aC, 


a l “ a ^l c 


A 0 


V 


-aC 


a 0 aC o c 
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6.1 


r 

0 


L° 

b „ 


r_ 

b g r st 


b G + ^ D G R ST 


Rigid Control System 

Frequently, the rotor is modeled as having a rigid control system. This 
option requires some restructuring of the equations of motion, for the rotor 
equations have been derived assuming that the blade rigid pitch degrees of 
freedom are present in the model and that the blade pitch control inputs enter 
through these degrees of freedom. A rigid control system is the limit of 
infinite control system and blade torsional stiffness. Thus the rotor blade 
elastic torsion motion is zero, and the solution of the rigid pitch equation 
of motion reduces to 

p c - e c ° n - £ Kj, q + Kj, S G + <e is cos - e lc sin * B >* 8 
i G 

or, in the nonrotating frame, 



(the result for N ^ 3 is similar). The blade pitch motion in this limit 
consists of the control input 0 con , feedback of the bending and gimbal motion 
due to the kinematic coupling, and a pitch change due to the azimuth perturba- 
tion with a fixed swashplate. Thus it is first necessary to account for the 
pitch/bending, pitch/gimbal, and pitch/azimuth coupling, which requires only 
operations on the columns of A Q (as indicated by the above equations) . Next 
the control matrix B is reconstructed from the rigid pitch columns of Ag 
since the blade pitch motion becomes a control variable rather than a degree 
of freedom. Finally, the equations of motion for the rigid pitch degrees of 
freedom may be dropped from the system. 


6.2 Quasistatic Approximation 

A quasistatic approximation is often used in rotor dynamics to reduce the 
order of the system. In the present analysis of the rotor in a wind tunnel, 
the quasistatic option is applicable to the inflow dynamics and sometimes to 
the blade pitch/torsion degrees of freedom. Let us assume that the equations 
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of motion have been reordered so that the quasistatic variables (x Q ) appear 
last in the state vector: 



The quasistatic approximation consists of neglecting the acceleration and 
velocity terms of the quasistatic variables. Thus the equations of motion 
take the form: 



Frequently, the Xj acceleration and velocity terms in the x Q equations 
will also be neglected (A^ 1 = = 0). The quasistatic variables now are no 

longer described by differential equations but rather by linear algebraic 
equations. The solution for Xq then is simply 

x Q = [Ag° ]” 1 [B°v - A^Xj - Aj 1 x 1 - Aq^Xj ] 


Substituting for Xq in the Xj equations of motion gives then the reduced- 
order equations for the quasistatic approximation: 

[A^ 1 - aJ°(aO°)~ 1 A° 1 ]x 1 + [A* 1 - Aq° (A§°)“ 1 A!j 11 ] + x^aJ 1 - A* 0 (Ag 0 )"^ 1 .]xj 

« [B 1 - Aj 0 (Ag°)- 1 B°lv 


The quasistatic approximation retains the low-frequency dynamics of the Xq 
variables. Whether that is a satisfactory representation of the aeroelastic 
behavior should always be verified by comparison with the results of the 
higher order model. 
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PART II. AEROELASTIC ANALYSIS FOR A ROTORCRAFT IN FLIGHT 


An aeroelastic analysis for a general two-rotor aircraft in steady-state 
flight is now developed (fig. 13). The rotorcraft to which the analysis is 
applicable include single main-rotor and tail-rotor helicopters, tandem-rotor 
helicopters, and side-by-side or tilting proprotor aircraft. Section 7 
describes the rotorcraft configuration considered. In section 8, the equa- 
tions of motion for the helicopter fuselage are derived, including both rigid 
body and elastic motions of the aircraft. The aerodynamic forces of the 
aircraft wing-body, horizontal tail, and vertical tail are included. A simple 
model for rotor-fuselage-tail aerodynamic interference (trim and perturbation) 
is developed. The rotor model used is that developed in part I; extensions 
of the model required for the helicopter in flight (principally in the inflow 
and transmission models) are developed in section 9. Finally, in section 10, 
the rotor and aircraft body equations are combined to construct the equations 
of motion for the coupled system. The analysis for the side-by-side or tilt- 
ing proprotor configuration is simplified if complete lateral symmetry is 
assumed in both the aircraft and the flight state. In that case, the longi- 
tudinal and lateral motions separate, allowing the solution of two lower order 
problems. 


7 . ROTORCRAFT CONFIGURATION 


The rotorcraft configuration features that have a major role in the 
dynamic behavior are the aircraft velocity and orientation and the position 
and orientation of the rotors. The aircraft flight path is usually specified, 
and by a trim calculation in which zero net force and moment on the aircraft 
are achieved, the control positions and aircraft orientation are determined. 
The orientation and position of the rotors are fixed geometric parameters. 


7.1 Orientation 

7.1.1 Flight-path and trim Eulev angles— The aircraft flight path is 
specified by the velocity magnitude V and the orientation of the velocity 
vector with respect to earth axes. The velocity vector has a yaw angle i/>pp 
(positive to the right) and a pitch angle 0pp (positive upward). Thus the 
climb and side velocities of the aircraft are V c pi m b = V sin 0pp and 
Vgide = V sin ^pp cos 9pp. The aircraft attitude with respect to earth axes 
is specified by the trim Euler angles, first pitch 0pp (positive nose-up), 
then roll <|>pp (positive to the right). Airplane convention is followed here 
for the coordinate systems — x positive forward, y positive to the right, 
and z positive downward (see ref. 15). Between the earth axes (E system) 
and the velocity axes (V system) are rotations ijjpp and 0pp. Between the earth 
axes and the body axes (F system) are rotations 0pp and <j>pp. Thus the rota- 
tion matrix between the V system and the F system is 
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Ca C C, 

9 FT ®FP ^FP 

-c a s , () 

9 FT ^FP 

C 0 S 0 C * 

9 FT 9 FP ^FP 

+ S 0 S fl 
9 FT 9 FP 


" s 0 c e 
9 FT 9 FP 

S. S A C o C , 

e FT <f>p T e pp ^pp 

-s. s, S, 
9 ft '’’ft ^fp 

S fl S * S fl C ,i, 

9 FT ^FT 9 FP ^FP 

+ C. C. S. 
^FT 9 FP ^FP 

+ C. c, 
^FT ^fp 

+ C. S Q S. 
'’’ft 0 FP ^FP 

- C fl S A S fl 

9 FT ^FT 9 FP 


+ C 0 S d> C 0 

9 ft St 9 fp 

S fi C * C fl C ,i, 

9 FT ^FT 9 FP ^FP 

-S„ c. s, 

9 FT ^FT ^FP 

s„ c. S Q c. 
9 ft St 9 fp Sp 

- S, C Q S, 
^FT 9 FP ^FP 

^FT ^FP 

- s. S Q s, 
St °FP Sp 

" C fl S fl 

9 ft ^ft 9 fp 


+ c fl c C 0 

9 FT St 0 FP 


The trim calculation determines the Euler angles 0pp and <f>pp (and, perhaps, 
the flight-path climb angle 0pp also). 

The velocity of the aircraft is V = Vi^, so the components in the body 
axes are 



v V t v 


The acceleration due to gravity is g = gk or, in body axes. 


E 

/- sin 0 


FT 


g = gk w = g[ cos 0 pT sin <|) 


k FT 


V° s 0 p T cos 4> FT y 


7.1.2 Rotor position and orientation— The rotor hub position is speci- 
fied in the bod^; axes^relative to the aircraft center-of-gravity position, 

r, , = (xi„ + yi„ + zk^), , . The rotor orientation is defined by a rotation 
hub v F /J F F hub 7 

matrix between the shaft axes (S system) and the aircraft body axes (F system) 
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The shaft axis components of the velocity seen by the rotor are then: 



The hub plane angle-of-attack and yaw angle may then be obtained from 



tan“l y 
tan -1 y /y 

y x 



The sign of the lateral velocity y^ must be changed for a clockwise rotating 

rotor (section 3.5), and if the induced velocity is included, the inflow ratio 

is X = y„ + X . . 

z 1 

For a helicopter main rotor, the orientation with respect to the body 
axes is specified by a shaft angle of attack 0 R (positive for tilt forward) 

and a roll angle <J> R (positive to the right). Thus, 

- c e 0 _s e 

" S (j) S 0 C cJ> S <t> C 0 

C ( f ) S © S < J > " C <|> C 0 

The orientation of a tail rotor is specified by a cant angle (j>^ (positive 

upward) and a shaft angle of attack 0 (positive forward). The tail rotor 

thrust is to the right for a counterclockwise rotating main rotor and to the 

left for clockwise rotation. Thus the definition of the tail-rotor shaft axes 

depends on the main rotor rotation direction. Let Q be +1 for a 
. r mr 

counterclockwise rotating main rotor -1 for clockwise rotation. Then the 
rotation matrix for the tail rotor is 
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The nacelle and rotor of a tilting proprotor aircraft can be tilted by an 
angle oip, where = 0 for axial flow (airplane configuration), and 
oip = 90° for edgewise flow (helicopter configuration). The rotor orientation 
is also described by a cant angle <j>^ (positive outward in helicopter mode, 
zero in airplane mode) and a pitch angle 9p (positive nose upward), which is 
the angle of attack of the shaft with respect to the body axes when dp = 0. 
Thus the rotation matrix is 



-c.s c a - c s. 

<j> a 0 a 0 

-s.s 

<j> a 

c.s s. - c c 0 

<j> a 9 a 0 

Rot. = 

SF 

-C.S. (1 - c )c. - s.s s ft 

<J) ij) ct 0 <J> a 0 

cl + sj-c 

(p <f> a 

C.S. (1 - C )S. - s.s c a 
<p <J> a 0 cj) a 0- 


(C*C + s2)C. - c.s s. 
(}) a (j) 0 a 0 


-(C^C + sf)s. - c.s c. 

<t> a <j>' 0 <j> a 0 


The rotor hub location ?jjub for the tilting proprotor aircraft is 
defined by the pivot location ^piyot an< * t * ie mas t height h, so that 


r hub 



7.1.Z Gust orientation— The aerodynamic gust components are defined with 
respect to the wind axes for the analysis of the rotorcraft in flight. The 
vector of gust velocity seen by the aircraft is 



where uq is positive downward, vq is positive from the right, and wq is 
positive upward. The aircraft aerodynamic forces are derived for gust com- 
ponents in the body axes, gp. Hence the transformation gp = Rpyg is 
required. The rotor aerodynamic forces were derived considering gust com- 
ponents in the shaft axes, g s . The substitution g s = Rgg, where 
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r g = 


"-1 0 o ' 

0 1 0 

0 0-1 


R R 
SF FV 


into the equations of motion then transforms the gust components to the wind 
axes . 


7.2 Pilot Controls 

The control variables included in the rotorcraft model are collective 
and cyclic pitch of the two rotors, engine throttle 9 t , and the aircraft 
controls (wing flaperon angle 6f, wing aileron angle 6 a , elevator angle S e , 
and rudder angle <S r ). The control vector is thus: 

V 1 - r(e co XrC'’) (e“ n e“ n e“ n ') tu « s 1 
[V o 1C IS /iV o 1C IS /2 t f e a rj 

The pilot controls, however, consist of collective stick <5 0 (positive upward), 
lateral cyclic stick 6 C (positive to the right), longitudinal cyclic stick 
6 S (positive forward), pedal <Sp (positive yaw right), and throttle lever 6 t : 

'i 

A linear relationship between the control inputs of the pilot and the rotor 
and aircraft control variables is. assumed: 


where v D is the control input with all sticks centered (vp = 0), and Tgp E 
is a transformation matrix defined by the control-system geometry. This 
transformation is required to obtain the aircraft response to the control 
input of the pilot. In addition, it is the pilot controls that must be 
adjusted to trim the rotorcraft. (The control transformation matrices for the 
single main-rotor and tail-rotor helicopter, the tandem main-rotor helicopter, 
and the side-by-side or tilting proprotor configurations are given in 
appendix C . ) 

7.3 Aircraft Trim 

The construction of the equations of motion for the rotorcraft dynamics 
must be preceded by a trim calculation, which determines the aircraft control 
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settings and orientation required for the specified equilibrium flight condi- 
tions. Equilibrium flight requires that the net force and moment on the 
aircraft be zero, which gives six equations to solve for the six trim vari- 
ables, consisting of four pilot controls and the two trim Euler angles (6 Q , 
6 C> <5 S , <$p, 0pT> an d ^ft)' This procedure is for level flight (8pp = 0) or 
a specified climb velocity. If, instead, the power available is specified, 
such as for power-off descent, then an additional trim variable (flight-path 
angle 0pp) and an additional equation (the power required equals the power 
available) must be included in the trim calculation. 

The balance of forces and moments about the aircraft center of gravity 
and the balance of power give the trim equations. The contributions to the 
forces and moments are from aircraft weight, aircraft aerodynamic forces, and 
hub reactions of the two rotors. In helicopter coefficient form, the force, 
moment, and power equations are 




( Y NI b ^/R) 2 /C H C, 

(yNI b ^ 2 /R) 1 \a 1 S a 3 S + a 



rotor 




1 +^kt + \(^t + kt +^lt 

, a S + a J S a k sl + \ hub L S + a J S + a k S 


{ (yNI, a 2 ) / c m c m c a 

— .2 (_X -f _1| 9. £ 

(yNI. fi 2 )i \ a S a J s a S 


rotor 1 


(YNI b ^ 2 /R) 2 
( Y NI b fi 2 /R) i 


I— (- 

I 2aAR \ 


i f \t 

1 \ r hub X j l 


a 1 S + a 3 S + a k S 


M x -t + ^ + /zd + 'V 2 
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q 1 V + q q SJ 
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l q V q V/ 
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( Y NI b fi 3 ) 2 



rotor 2 


"rotor power available 


The components of the force and moment equations are obtained in the body 
axes (F system). Here W is the aircraft weight; and the hub reactions for 
rotor 2 are normalized using the parameters of that rotor, hence the factors 
accounting for the normalization of the coefficients. The aircraft aero- 
dynamic forces are acting on the wing-body (WB) , horizontal tail (HT), and 
vertical tail (VT). Here L, D, and Y are, respectively, the aerodynamic 
lift, drag, and side forces; M x , My, and M z are the roll, pitch, and yaw 
moments on the wing-body and q is the dynamic pressure. 


A consideration of the aerodynamic interference between the rotors, wing, 
and tail is required to accurately calculate the trim state. A simple model 
for this interference is used here. The rotor- induced velocity, together 
with the aircraft velocity, is used to determine the angle of attack at the 
wing and tail. For the horizontal tail, the angle-of-attack change due to the 
wing wake is also included. The rotor- induced velocity is assumed to be 

directed along the rotor shaft (kg). A multiplicative factor on the induced 
velocity is used to account for the fraction of the aerodynamic surface within 
the wake and the fraction of the fully developed wake velocity achieved. A 
further multiplicative factor accounts for the decrease in the wake-induced 
velocity away from the wake boundaries (see, also, the discussion of the 
perturbation aerodynamic interference model, section 8.2). The angle-of-attack 
change at the horizontal tail due to the wing is calculated by 


Aa 


0.45 C 

L WB 

a 2 / s )°- 735 a HT /c )°* 25 
w w HT w 


where C^g is the wing lift coefficient; S w , £ w , and c w are, respectively, 

wing area, span, and chord; and is the tail length (ref. 16). Alterna- 

tively, all the interference effects could be included in the wing-body, 
horizontal tail, and vertical tail ^aerodynamic characteristics. 

The trim equations are nonlinear in the control variables, of course. 

Thus an iterative solution procedure is required in which the control variables 
are incremented in the direction to achieve zero net force and moment, based 
on a set of local partial derivatives obtained at the beginning of the trim 
calculation by making step changes in the individual control variables. The 
solution is considered to have converged to the desired trim state when the 
net force and moment are within a certain tolerance level. 


8. AIRCRAFT MODEL 


The aeroelastic motion of the rotorcraft airframe is described by a set 
of linear, constant coefficient differential equations, excited by the hub 
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reactions of the two rotors. Let x g be the vector of the aircraft degrees 
of freedom, v s the vector of the aircraft control variables, and gp the 
vector of the aerodynamic gust components. The equations of motion for the 
rotorcraft in flight are required in the following form: 

a 2 X s + a i*s + a 0 X s = bV s + b G g F + aF 
and the hub motion is given by 

a = cx 

s 

Here F and a are as usual the rotor hub reactions and hub motion in the 
shaft axes (S system) : 



(For convenience, only the terms for one rotor are shown, but, in fact, the 
interface between the aircraft and the rotor is required at both hubs. The 
parameters of rotor 1 are used to make quantities dimensionless and to 
normalize the aircraft equations of motion.) 

In this section, the aircraft equations of motion are constructed in the 
required form. The aircraft degrees of freedom (x s ) consist of the six rigid- 
body degrees of freedom and the elastic free-vibration modes. The input 
variables (v s ) consist of the aircraft aerodynamic controls — flaperon, 
aileron, elevator, and rudder. An elementary model for rotor-wing-tail 
aerodynamic interference is also developed. 

A body axis coordinate frame with its origin at the aircraft center of 
gravity (F system) is used to describe the motion. Airplane practice is fol- 
lowed for these axes — x forward, y to the right, and z downward (ref. 15). 
The coordinate frame used is not a principal axis system, however. Moreover, 
the airplane practice of aligning the x axis with the trim velocity is not 
followed since, for rotorcraft, the hover case (V = 0) must be considered. 
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Lateral symmetry of the aircraft inertia and of the aerodynamic surfaces is 
assumed; the location and orientation of the two rotors is entirely general, 
however. 


8.1 Aircraft Motion 


8.1.1 degrees of freedom- The linear and angular rigid-body motion of 
the aircraft is defined in the body axes (F system). The linear degrees of 
freedom are xp (positive forward), yp (positive to the right), and zp 
(positive downward). These variables are dimensionless based on the rotor 
radius R. Thus the velocity perturbations are normalized using the rotor 
tip speed ftR rather than the forward speed V as is airplane practice. 

The angular degrees of freedom are the Euler angles ^p (yaw to the right), 

0p (pitch nose-up), and <f>p (roll right). Then the linear and angular velocity 
perturbations are 


where 


U F X F D 'F + ^F^F * Z F^F 
W F = R e(^F^F + ®F J F + V*t) 



’ 1 

0 

-sin 

6 FT 


R = 
e 

0 

cos <j>p T 

sin * FT 

cos 

9 FT 


0 

-sin <(> FT 

cos * rI 

cos 

9 ft 


For the elastic motion of the aircraft in flight, the displacement u 
and rotation ^ at an arbitrary point r are expanded in series of the 
orthogonal free vibration modes: 


u(r,t) 



(t)t k (r) 


e(r,t) = X) q s 
k=l k 


The first six degrees of freedom are the rigid body motions: q Sl ...» q Sg , 

are, respectively, <f>p, 0p, i^p, Xp, yp, and Zp. The generalized coordinates 
q for k = 7 to 00 are the elastic modes of the aircraft. Orthogonality 
implies that the elastic vibration modes produce no net displacement of the 
aircraft center of gravity or rotation of the principal axes. 
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For the rigid-body motions, the mode shapes are simply 

[?! ••• Se] - [<-?x)Rjl] 

[yi • • ■ Ye] - [rJo] 

8.1.2 Hub motion— The linear and angular motion at the rotor hub in the 
shaft axes (S system) is then 



= cx 

s 


The total velocity .of^a point is the sum of the trim and perturbation 
velocities, u = V + y qs^k* in body axes. The rotor equations require the 

velocity components at the hub in an inertial frame, however (S system), and 
the Euler angle rotations between th^ body and inertial axes introduce 
perturbations of the trim velocity V. So the perturbation velocity becomes 
4» “V ^ ->• 

u = ctpXV + i_j q g where, in the S system, 
k 

Xa - y a 
y y z 

-Xa - y a 
X X z 


y a + y a 
x y y x 

It follows that the (x^jy^jZ^) columns of the rotor matrices Aj and Cj con- 
tribute to the (a x ,oty,a z ) columns of and C Q , which exactly cancel the 
existing terms due to the rotation of the inertial axes relative to the 
velocity components y and X (see the matrices in section 2.6). The result 
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is that the net angular hub motion columns of A Q and C Q are zero for the 
Euler angles. 

Furthermore, the acceleration is u = WpxV + 2 w ^ ere the second 

term is the inertial acceleration due to the rotation of the trim velocity 
vector by the body axes angular velocity. The components of this additional 
acceleration in the shaft axes system are 



In summary, the hub motion is a = cx s (where the matrix c is given 
above), with two exceptions. First, for the Euler angles, the net (a x ,a y ,a z ) 
columns of the rotor matrices A Q and C 0 are zero because of the use of body 
axes. Second, in a there are additional linear acceleration terms due to 
the Euler angle velocities, Aa = cx s (where only the upper righthand 3x3 sub- 
matrix of "c is nonzero). 


8.1.3 Equations of motion and hub forces— Following the usual steps of 
airplane flight-dynamics analyses (ref. 15), the rigid-body equations of motion 
are obtained by equating the angular and linear acceleration to the net moments 
and forces on the aircraft: I&p = M and MC&p + tft-FrxV) = T! F. In terms 
of the body axes degrees of freedom then, including the gravitational forc% 
the equations are 





Here M is the aircraft mass (including the rotors) and I, the moment of 
inertia matrix: , : 
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(I = I 
xy yz 


0 since lateral symmetry is assumed). The equations are dimen- 


sionless and are normalized by dividing by the characteristic inertia (N/2)I|j 
(using the parameters of rotor 1). Thus M* = M/ [ (l/2)NIb/R 2 ] and 
1* = 1/ [ (l/2)NIb] . Note that the gravitational constant g is also dimen- 
sionless based on the acceleration ft 2 R. 


For the elastic degrees of freedom, since orthogonal free-vibration modes 
are used, the equations of motion are simply 



+ w. 


+ »kV 


k > 7 


where Mj^ is the generalized mass, the normalized mass is 

Mj| = Mfc/ [ (l/2)NIb/R 2 ) ] , is the natural frequency of the mode, and g s is 
the structural damping coefficient. 

There are two sources for the generalized forces Q^: the rotor hub 

reactions and the aerodynamic forces on the aircraft. The generalized force 
due t£ tlje rotor_^hub reaction is 

Qk 58 tk( r hub) * ^hub + fk^hub) * M h u b* Normalizing Q k by dividing by 
(N/2)Ib gives then {Qj£} = aF, where 



The aircraft aerodynamic forces are obtained in section 8,4, 
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8.2 Aerodynamic Interference 

The interference between the rotor wake and the aircraft aerodynamic sur- 
faces (wing and tail) can be a factor in the dynamic behavior. As a simple 
model of this aerodynamic interference, it is assumed that there is a pertur- 
bation velocity at the wing, the horizontal tail, and the vertical tail, which 
is a linear combination of the perturbation induced velocities from the two 
rotors (Arj and Ar^). Including a first-order time lag, the equations for 

the interference velocities are then 

m ) 2 

T w A w + A w ~ K W 1 C W 1 A R 1 + ^2^2 ( fiR )l Ar 2 

(«R) 2 

t A + A •— K C A .+ K" C„ / 0 n \ ^-d 

H H H Hj rl 2 ^2 (^R) 1 

(OR) 2 

T A + A = K C T7 A + K C Ycyd \ — ^■•n 

v v v Vi Vi Ri V 2 V 2 (^R)i R^ 

A time lag of x = £/V is used, where V is the aircraft velocity and i 
is the distance between the aerodynamic surface and the dominant rotor. 

The first multiplicative factors (K) account for the maximum fraction of 
the aerodynamic surface affected by the wake and the fraction of the fully 
developed wake velocity achieved. Typical values would be K = 1.5 to 1.8 
(or 0 for no interference) . The second multiplicative factors (C) account 
for the cosine of the angle between the wake axis and the normal to the 
aerodynamic surface, and the decrease in the wake-induced velocity away from 
the wake surface. The following expression is used: C = (cosine of angle 

between wake and surface) / (maximum of 1 and 1 + L) , where L is the perpen- 
dicular distance from the aerodynamic surface to the nearest wake boundary 
(L < 0 if the surface is inside the rotor wake cylinder). 


8.3 Aircraft Equations of Motion 


The equations of motion for the aircraft in flight may now be written: 

a 0 x + a.x + a n x = Q . + Q 

2s Is Os ^rotor aero 

The vector of the aircraft degrees of freedom consists of the six angular and 
linear rigid-body motions, the generalized coordinates of the antisymmetric 
and symmetric elastic modes, and the aerodynamic interference velocities: 
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*s- anti 


s k Sym 


w 


H 


The generalized force due to one rotor is Qrotor = ^F, and the hub motion 
for one rotor is a = cx g . There are additional linear acceleration terms 
due to the Euler angle velocities given by Aa = (!x g . 


The matrices c and c are defined in section 8.1.2; a is given in 
section 8.1.3 (note that a can be obtained directly from c). The inertia, 
structural, and gravitation forces are included in the matrices of the equa- 
tions coefficients (appendix Dl). 


8.4 Aerodynamic Forces 

To complete the aircraft equations of motion, the aerodynamic forces 
acting on the wing-body, horizontal tail, and vertical tail are required. 
Helicopter airframe aerodynamics involves complex nonlinear phenomena, 
particularly significant aerodynamic interference such as between the tail and 
the rotor and fuselage. It is difficult to include such effects in any simple 
model. For best results, therefore, experimental data should be used as much 
as possible, but often such data are simply not available. Thus analytical 
expressions for the aerodynamic stability derivatives are required. 

The aerodynamic forces on the wing and tail are calculated by a strip 
theory analysis. The generalized force is obtained by integrating the section 
lift and drag forces and the section moment along the span. Three-dimensional 
effects are accounted for in the integrated aerodynamic characteristics used 
for the wing and tail. A body axis system is used, but with the x axis not 
aligned with the aircraft velocity vector. Otherwise, the analysis follows 
the usual derivation of airplane stability derivatives (see ref. 15). 
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Lateral symmetry is assumed for the aerodynamic forces. Specifically, 
it is assumed that the trim velocity and the center of action of the wing 
and tail are in the x-z plane. Then the symmetric and the antisymmetric modes 
of the airframe are not coupled by the aerodynamic forces. 

The aircraft motion consists of the rigid body and elastic degrees of 
freedom. Consistent with the strip theory analysis, the wing elastic motion 
is described by vertical and chordwise bending and torsion, including wing 
root motion due to the fuselage flexibility. For the kth symmetric or anti- 
symmetric mode of the airframe, the wing motion is thus described by vertical 
deflection zj c (£) (positive upward), chordwise deflection x^Cil) (positive 
aft), and torsion or pitch 0i c (£) (positive nose-up), where SL is the span- 
wise coordinate (J£ = 0 at the root, and SL = ±(1/2 )£ w at the wing tips). 

For symmetric modes, x^, zy., and 9^ are nonzero at the root due to the 
fuselage motion (but x£(0) = z k(0) = 0)* For antisymmetric modes, x£(0) 
and z k (0) are nonzero, while x k (0) = z k (0) = 9^(0) = 0; in addition, the 
fuselage motion gives a lateral reflection of the wing y^ (positive to the 
right) . 

For the tail motion, only rigid linear and angular motion due to the 
fuselage flexibility is considered; bending and torsion of the tail surfaces 
are neglected. Thus the horizontal tail motion for symmetric modes is 
described by vertical deflection z^ (positive upward) , longitudinal deflec- 
tion x^ (positive forward), and pitch Qj. (positive nose-up). The vertical 
tail motion for antisymmetric modes consists of lateral deflection y^ (posi- 1 
tive to the right), roll <j>k (positive roll right), and yaw ^ (positive yaw 
right). There is no vertical tail motion in symmetric modes; the horizontal 
tail motion in antisymmetric modes is just roll <t>k (positive roll right). 

The aircraft controls considered are wing flaperon deflection 6f and 
aileron deflection S a (symmetric and antisymmetric motion of the wing control 
surfaces), horizontal tail elevator deflection S e , and vertical tail rudder 
deflection 6 r . Aerodynamic forces due to the three gust components (in the 
F body axis system) are included. 

The aerodynamic forces on the aircraft, required to complete the equa- 
tions of motion in section 8.3, take the following form: 

Q = -a 0 x - a.x - a„x + bv + b_g_ 

^aero 2 s Is Os s G e F 

The vector of aircraft degrees of freedom x s is defined in section 8.3. 

The vector of the aircraft controls v g and the components of the gust vector 
gp (in the F system) are 
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The matrices of the aerodynamic coefficients are given in appendix D2. 
Expressions for the aerodynamic coefficients required in these matrices are 
given in appendix E. 


9. ROTOR MODEL DETAILS FOR THE FLIGHT CASE 


To treat the general twin-rotor helicopter, a number of extensions of the 
rotor model are required, principally in the models for the inflow dynamics 
and the rotor speed dynamics. Rotor-rotor aerodynamic interference is con- 
sidered, in both the trim- and perturbation- induced velocities. Ground effect 
is also included in the inflow dynamics model. Pitch/mast-bending coupling is 
introduced. A transmission and engine model for two interconnected rotors is 
derived. The drive train dynamics are described by the rotor speed, inter- 
connect shaft torsion, and engine shaft torsion degrees of freedom. The 
throttle control variable is introduced. Finally, a governor with collective 
or throttle feedback of rotor speed is considered. 


9.1 Rotor-Rotor Aerodynamic Interference 


With twin-rotor aircraft, it is necessary to account for the rotor-rotor 
aerodynamic interference in both the trim- and perturbation- induced veloc- 
ities. The model used expressed the induced velocity at each rotor as a 
linear combination of the isolated rotor-induced velocities. Let Aj, and 
A-?„ be the trim-induced velocities of the two isolated rotors (calculated 


12 

as described in section 2.3.3). 


Xi - 


Then the trim inflow ratios are 
(QR) 2 


v 9 + X, 

Zl 1 


+ K 


1 


12 


*2 = + X i2 + *21 


z 2 


m) i 

m)i 

(8R), 


X. 

12 


Here k 12 an d K 2l are the rotor-rotor aerodynamic interference factors. 
Separate values are used for the interference factors in hover and forward 
flight, with a linear variation from p = 0.05 to 0.10. 


Similarly, the rotor-rotor interference is included in the uniform inflow 
perturbation. Recall from section 2.5 that the differential equation for the 
inflow dynamics of the isolated rotor is 


tX + X 


- (y~) 

9L V Y aal 

\ 'aero 


where (8X/9L) 1 = (2y/oa) (X /k 2 + A 2 /k£ + y 2 /i< 2 ) . With two rotors, the 

on on i 

inflow perturbations at one is a combination of the influence of both rotors; 
hence the differential equations become 
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(^R.) g / dX \ / C t\ 
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The interference factors k^ 2 and k 2 i are the same as for the trim- induced 
velocity. 


For the side-by-side or tilting proprotor aircraft configurations, 
lateral symmetry gives <j 2 = < 2J = k. Then the trim-induced velocity is 
A = y z + (1 + k)A^. The differential equation for the inflow perturbation 
becomes 

tA + A. = (1 + k) |~ h ~“) 

\ /aero 

for symmetric dynamics and 


tA + A = (1 - k) 


9A 

3L 



for the antisymmetric dynamics of the aircraft. 


9.2 Ground Effect 

To account for the effect of the ground on the rotorcraft dynamics, it 
is necessary to correct the trim- and perturbation-induced velocities for 
the proximity of the ground plane. Based on reference 17, an approximate 
expression for the ratio of the induced velocities in and out of ground effect 
is 



1 - 


1 


( 42 eff ri 


Here z e ff ~ z/cos e, where z is the altitude of the rotor hub above ground 
level, normalized by the rotor radius, and e is the angle between the ground 
and the rotor wake: 


(vlX - yjo - xfc e ) 


cos e “ 


x S 


y J S 


k. 
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y 1 2 - + y z + A^ 
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|2 
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which thus accounts for the effect of forward speed. Note that ground effect 
is essentially negligible for altitudes greater than the rotor diameter, or 
at forward speeds y > 2(C>j/2) V 2 . This expression compares well with test 
results, down to an altitude of about 1/2 rotor radius (see ref. 17). The 
trim- induced velocity in ground effect (before incorporating the rotor-rotor 
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interference) is thus 


X = y + 
z 


cos^ e 1 
16z 2 , 


The effect of the ground on the inflow dynamics is to add a perturbation 
due to changes in the rotor height above the ground: 


xX + X 


(-S), 


+ 5Z 


where, again based on the results in reference 17, 

X^ cos 2 e 
o A O 


(Actually, the ground-effect term is added to the equations for rotors 1 and 
2, including the rotor-rotor interference terms, as in section 9.1.) The 
rotor height perturbation Sz is due to the rigid body and elastic degrees 
of freedom of the airframe. The vertical component of the displacement at 
the rotor hub gives 

Sz ■ -%s • ( Vs + y Js + Z h 5 S> 

■ ( Z hub C0S 0 FT sIn *FI ‘ y hub C ° s e FT COS ♦ftK 

+ [( Z hub C0S *FT + y hub Sln +Fr) Sl " 6 FI + Vb C ° S 6 Fl] 6 F + (sin 6 FT >X I 

00 

+ (-cos e FT sin <J> FT )y F + (-cos e pT cos <j> pT ) z p + £ * ^q g 


Since 9X/3z > 0, ground effect introduces a positive spring to the rotorcraft 
height dynamics (zp perturbations). A decrease in the rotor height above the 
ground produces a decrease in the induced velocity, hence a rotor thrust 
increase that acts as a spring against the vertical height change. 

For the side-by-side helicopter configuration, the antisymmetric dynamics 
exhibit an unstable roll oscillation due to interaction of the rotor wake and 
the ground. Such behavior can be included in the ground effect model derived 
here by using a negative value for 9X/9z (a negative roll spring), which must 
be obtained from experimental data. 


9.3 Pitch/Mast-Bending Coupling 

Flexibility between the rotor swashplate and hub will produce a blade 
pitch change due to elastic motion of the airframe. This coupling between 
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the rotor pitch and the mast-bending will be accounted for by introducing 
kinematic feedback terms of the following form: 



9.4 Transmission and Engine Dynamics Model 

The rotor rotational speed degree of freedom can be an important factor 
in the rotorcraft flight dynamics. A model is required which accounts for 
the coupling of the two rotors through the flexible drive train, and for the 
engine damping and inertia. The throttle control of engine torque must also 
be introduced. Figure 14(a) is a schematic of the transmission-engine model 
used for the single main-rotor and tail-rotor, and the tandem main-rotor 
helicopter configurations. The two rotors are connected by a shaft, and the 
engine is geared to one rotor (rotor 1 in fig. 14(a)). The torsional flexi- 
bility of the drive train is represented by the rotor shaft springs K]^ and 

%2 5 t * ie interconnect shaft spring Kj, and the engine shaft spring Kg. The 

transmission gear ratios are rg (ratio of the engine speed to rotor 1 speed), 
and rj 1 and rj 2 (ratio of the interconnect shaft speed to rotor speed). Thus 
r I 1 / r I 2 is the ratio of the trim rotational speeds of rotors 2 and 1. 

The^degrees^of freedom are the rotational speed perturbations of the two 
rotors (if s and ^ S2 ) and the engine speed perturbation ($ e ) . The engine shaft 

azimuth perturbation \jj R is defined relative to rotor 1 rotation, so the 
total engine speed perturbation with respect to space is ^e^si + ^e) * With 

coupling of the speeds of the two rotors by the drive system, it is more 
appropriate to use the degrees of freedom: 



Here is the differential azimuth perturbation between the two rotors. 

The degrees of freedom ijj-j- and \jj Q therefore involve elastic torsion in the 
drive train (in the interconnect shaft and engine shaft, respectively) and so 
represent high-frequency modes. The degree of freedom \ p s is the rotational 
speed perturbation of the drive system as a whole — both rotors, the engine, 
and the transmission. 

The engine model includes the inertia, damping, and control torques: 



Here % is the engine speed and Qe is the perturbation torque on the 
engine. The engine rotary inertia is Ig; is the engine speed damping 

coefficient, that is, the torque per unit speed change at constant throttle 
setting (see section 2.4). The variable 0 t is the engine throttle control 
position and Q t is the torque applied due to a throttle change at constant 
speed: 
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Thus Q t and can be obtained from data on the engine power as a function 

of throttle position and engine speed. 


The differential equations of motion for the rotor speed dynamics are 
obtained from equilibrium of the torques on the two rotors and the engine. 
The resulting equations for ip s , \jjj 9 and ifj e are 
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r E K E ( K M 2 + r l2 Kl ) + Kmm 
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Sh - k m 1 k m 2 + + 

The spring constants are normalized by dividing by (NI^ft 2 )jj l| = Ig/CNI^). 
An alternative configuration for the transmission is with the engine by 
rotor 2, instead of by rotor 1 as in figure 14(a). The equations of motion 
for that case are obtained simply by exchanging subscripts 1 and 2 in the 
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three equations above; note that the definitions of the degrees of freedom are 
then: 



The normalized damping and throttle coefficients may be written 

r 2 0 P 
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The approximate expression for 


Qfi 


is discussed in section 2.4. 


The side-by-side or tilting proprotor aircraft requires a different trans- 
mission and engine model due to the lateral symmetry assumed for these con^- 
figurations. Figure 14(b) is a schematic of the model considered. The two 
rotors are connected by a cross-shaft, and there are two engines, one geared 
to each rotor. The degrees of freedom are the. rotor speed perturbation <j> s 
and the engine speed perturbation (defined relative to the rotor speed 

again). The equations of motion for \p s and if) e follow as above: 
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For antisymmetric motion, ip s is the differential azimuth perturbation between 
the two rotors, involving torsion in the, interconnect shaft. For symmetric 
motions, ; there is no torque on the interconnect shaft, so the above equations 
apply with Kj - 0 (so K MR = 1 and K EI = K MI = 0) . 

The case of a rotorcraft in autorotation can be treated with this model 
by dropping the ip e degree of freedom and dropping the engine terms from the 
and equations (helicopters usually have an overrunning clutch to 
disconnect the rotors from the engine at zero torque) . The engine-out case 
(engine and rotors still connected) can be handled by dropping the engine 
damping term. The case of constant rotor speed is modeled by simply dropping 
the and ip e equations and degrees of freedom from the system. Gener- 

ally, the ■ ijij and ij) e degrees of freedom are only involved with high-frequency 
dynamics, and so it is usually sufficient to consider the ip s degree of free- 
dom for flight dynamics analyses. 


9.5 Rotor Speed Governor 


When the rotor rotational speed perturbation is included in the flight 
dynamics analysis, it is usually necessary to also include the rotor speed 
governor in the model for a consistent calculation of the aeroelastic 
behavior. The governor model considered is integral and proportional feed- 
back of the engine speed to throttle and to the collective pitch of rotors 1 
and 2. The control equations are then 


46 1 = Kj, (* s + * e ) + Kj (* e + * e ) 


<A9 > , - 1L, (} + } ) + K t (* + * ) 

o rotor 1 Pj r s Y e 1^ s Y e 


(A6 ) 
o 


rotor 2 



+ K t 4 + 


+e> 


Note that ip s is the rotor speed error and \p s then is the integral of the 
error. Since the governor dynamics are neglected, this model does not add 
degrees of freedom to the aeroelastic analysis. 


10. COUPLED ROTOR AND BODY MODEL 


The equations of motion have been derived for the two rotors and for the 
aircraft body. Now these equations may be combined to construct the set of 
linear differential equations that describes the dynamics of the complete 
rotorcraft system. The equations of motion for the coupled model then have 
the following form: 



Here x is the vector of the degrees of freedom for the enrUe system, v is 
the vector of the individual control inputs, vp is the vector of the pilot's 
control inputs, and g is the aerodynamic gust vector (in velocity axes). 

A 2 x R + AiX R 4- A 0 x r + A 2 a + A]a + A 0 a = Bv R + B G g g 

F = C 2 x r + C jx R + C 0 x R + C 2 a + C^a + C 0 a + D G g g 

for rotors 1 and 2 (see part I; in particular, section 2.6). Recall that the 
vector of rotor degrees of freedom xp consists of the flap/lag bending, 
rigid pitch/elastic torsion, gimbal tilt, rotational speed, and inflow pertur- 
bation variables; the vector of the rotor controls vp consists of the blade 
pitch control inputs; and the aerodynamic gust vector g s is in the shaft 
axes for the rotor model: 



As usual, a and F are, respectively, the hub motion and hub reaction in the 
shaft axes. 

The equation of motion and the hub motion expressions from the aircraft 
model (section 8) are 

a 0 x + a.,x + a n x = bv + b^g- + aF 
2s Is Os s G°F 

a = cx 

s 


Act = cx 

s 

where Act is the linear acceleration due to the rotation of the velocity 
vector in body axes by the Euler angular velocity (section 8.1.2). The air- 
craft degree-of-freedom vector x s consists of the rigid-body angular and 
linear variables, the antisymmetric and symmetric elastic f ree-vibration 
modes of the airframe, and the aerodynamic interference inflow variables; the 
vector of the control inputs v s consists of the aircraft aerodynamic control- 
surface deflections; and the aerodynamic gust vector g R is in the body axes 
for the aircraft model: 
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The gust components for the rotorcraft model must be in velocity axes; 
hence substitutions g s = R^g and gp = Rpyg are required in the rotor and 
body equations (see section/. 1. 3) . The transmission and engine model (sec- 
tion 9.4) replaces the individual rotor speed perturbations and ij/ g 

by the coupled degrees of freedom \p s and iJjj, introduces the engine speed 
degree of freedom i^ e , and adds the engine throttle control 0^ to the model. 
The pilot controls Vp = (6 0 6 C 5 g 6p 5 t )^ are related to the rotor- 
craft input vector v by a linear transformation v = T^pgVp (see sec- 
tion 7.2). Then the state vector x, control vector v, and aerodynamic gust 
vector g fi 
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The coupled equations of motion are obtained by substituting the hub motion 
into the rotor equations and hub reactions, and then the hub reactions into 
the body equations of motion. The resulting coefficient matrices for the 
coupled system are 
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In constructing these matrices, it is necessary to skip the angular hub motion 
(a x , oty, a z ) columns of Aq and Cq for the Euler angles (d> F , 0p, \p-p) since 
body axes are used for the aircraft motion ^see section 8.1.2). Also, the 
linear hub motion (x^, yjp z^) columns of C 2 should be skipped for the body 
degrees of freedom, assuming that the rotor mass is already included in the 
aircraft gross weight and free-vibration generalized masses. 

The rotorcraft equations of motion are normalized based on the parameters 
of rotor 1 (ft, R, N, 1^,, y, a, etc.). The equations for rotor 2 as derived 
are, however, based on the rotor 2 parameters. Therefore, it is necessary to 
multiply the coefficient matrices for the rotor 2 equations of motion and hub 
reactions by appropriate scale factors to account for the differences in nor- 
malization. The degrees of freedom and control variables for rotor 2 will 
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still be normalized based on rotor 2 parameters. Most are angular variables 
anyway, hence inherently dimensionless. The components of the hub motion, 
hub reaction, and gust are based on rotor 1 parameters, however: the linear 

hub displacements in a are based on R* not R 2 ; the gust velocities are 
based on (ftR)j not (J2R) 2 ; and the forces and moments in F are based on 
(NI b fi 2 /R) 1 and (NI b fi 2 )^, respectively. Finally, the scale factors for the 
rotor 2 equations must account for the time scale of the complete system, 
which is based on the trim rotation speed of rotor 1. 

The equations for the rotor inflow dynamics are completed by accounting 
for the rotor-rotor aerodynamic interference (section 9.1) and the effect of 
the ground (section 9.2). The equations for the airframe-rotor aerodynamic 
interference variables (A w , Ajj, A y ) are completed after constructing the 
coupled equations of motion. Note that this aerodynamic interference is the 
only coupling of the rotor and body not taking place through the rotor hub. 
Pitch/mast-bending coupling is accounted for by adding terms for the elastic 
airframe degrees of freedom (q Sk ,k > 7) in the rotor rigid pitch equations 
(section 9.3). The rotor speed governor model is added to the system as 
described in section 9.5. Finally, the unused equations of motion and degrees 
of freedom may be eliminated from the model by deleting the appropriate rows 
and columns from the coefficient matrices. 


10.1 Rigid Control System 

A rigid control system model may be used for either or both rotors. In 
the limit of infinite control system stiffness, the equations of motion for 
the rotor rigid pitch degrees of freedom reduce to the algebraic relations: 



(the result for the number of blades N # 3 is similar). On the basis of 
this equation, the matrices A 0 and B are reconstructed as outlined in sec- 
tion 6.1. 


10.2 Quasistatic Approximation 

It is frequently possible to reduce the order of the system of equations 
describing the rotorcraft dynamics by considering a quasistatic approximation 
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for certain of the degrees of freedom. In the present analysis of the rotor- 
craft Iti flight, the quasistatic approximation is applicable to the inflow v 
dynamics of either or both rotors, to the rotor-body aerodynamic interferences 
variables, to the rigid pitch/elastic torsion degrees of freedom of either, or 
both rotors, to all the degrees of freedom for rotor 1, rotor 2, or both 
rotors, or even to all degrees of freedom except the six rigid-body motions 
of the aircraft. The reduction of the model by eliminating the quasistatic 
variables is described in section 6.2. 

The quasistatic rotor model is frequently useful, and often valid, in the 
analysis of rotorcraft flight dynamics. It Is usually a satisfactory repre- 
sentation for the tail rotor and may also be satisfactory for the main rotor-' 
dynamics for such applications as low-rate stability and control augmentation 
system investigations. Generally, whether the quasistatic model is a satisfac- 
tory representation of the aeroelastic behavior must always be verified 1 for a 
particular application of the analysis by comparison with the results of the 
higher order model. ' ■■■•' -s ■» 


10.3 Side-by-Side or Tilting Proprotor Configuration 

The aeroelastic analysis for the side-by-side or tilting proprotor air- 
craft configuration requires special consideration. Assuming complete lateral 
symmetry of both the aircraft and the flight state, the symmetric and anti- 
symmetric motions are entirely decoupled. Thus the analysis involves the 
solution of two problems of half the order of the whole system. The motions 
of the left and right sides of the aircraft are then given by, respectively; 
the sum and difference of the symmetric and antisymmetric degrees of freedom. 

The symmetry of the flight state requires i^p =0 (no side velocity) 
and (jjpj = 0 (zero trim roll angle). The trim solution has automatically 
6 C 88 6p = = 0, so it is only necessary to solve three equations (vertical 

and longitudinal force, and pitch moment) for three trim variables (<$ 0 , 6 S , 
and ©ft). The construction of the coupled differential equations of motion 
follows basically the steps outlined above for the general two-rotor helicop- 
ter. It is also necessary to obtain the equations of motion for one rotor, 
however, multiplying the hub reactions by a factor of 2 to account for both 
rotors of the aircraft (the renormalization for rotor 2 is not required since 
the two rotors are identical in this case). 


10.4 Two-Bladed Rotor Case 

The two-bladed rotor has special dynamic characteristics compared with 
the case of three or more blades. Generally, the dynamic behavior is described 
by periodic coefficient differential equations, so a Floquet analysis is s 
required except for special cases (such as a shaft-fixed rotor in axial flow; 
see section 4.2). For helicopter flight dynamics, the main concern is with 
the low-frequency impedance of the rotor hub reaction response to shaft 
motion, control inputs , and gusts. * . . 
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The impedance of a linear time- invariant (constant coefficient) dynamic 
system is described by a transfer function H(co): 

F = H(a))a 

that relates the magnitude and phase of the input and output at frequency w. 
The implication of the periodic coefficients of the two-bladed rotor is that 
such a transfer function relation does not exist, for an input at frequency 
a) generally produces a response at all frequencies co ± nfi, n = 0, . . . , 
Then the input-output relation takes the form: 



The flap or teetering response of the two-bladed rotor is found to be primarily 
at frequencies w ± ft. It follows that the low-frequency flap response is at 
±ft, so the low-frequency motion can be written: 

3 = B cos ip + B sin ip 

It is found that the solution for the 3ic an d 3 i s fl a P motion is identical 
to that for the rotor with N > 3 at low frequency. Furthermore, it is found 
that the average of the coefficients of the hub reactions at low frequency is 
the same for the two-bladed rotor as for N > 3. But while this constant 
coefficient result is exact for a rotor with three or more blades in hover, 
due to the rotor inertial and aerodynamic axisymmetry, for the two-bladed 
rotor there really are periodic coefficients in the hub reactions. Specifi- 
cally, there is a large 2/rev variation of the coefficients even in hover due 
to the rotor asymmetry when N =2. 

Difficulties also arise with the quasistatic approximation. As imple- 
mented in section 6.2, the velocity and acceleration terms in an equation are 
dropped, reducing that equation to"an algebraic substitution relation for the 
quasistatic variable. For a rotor with three or more blades, the quasistatic 
approximation applied to the equations in the nonrotating frame produces 
exactly the low-frequency response of the rotor. Note that it is necessary 
to consider both the and 3 i s equations even when only longitudinal or 

lateral dynamics of the helicopter are involved, for the Big and vice versa. 
For the two-bladed rotor, however, the quasistatic approximation does not give 
the low-frequency response because the g]_ equation is really in the rotating 
frame . 

In summary, the two-bladed rotor is indeed a special case. First, the 
description of the dynamics is unique, involving the teetering degree of 
freedom Bi» which is fundamentally in the rotating frame, rather than the 
cyclic degrees of freedom and gjg. The frequency response is not given 
by the common transfer function relation because the system is not time 
invariant. The low-frequency flap response does reduce to a tip-path-plane 
representation, identical to the result for N > 3; but the to = 0 limit, 
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which allows the 3q •- Sic cds f + 3 i s sin $ representation, is a special 
case. 


Second, the equation of motion for the helicopter flight dynamics, while 
the same as for N > 3 if the averaged coefficients are used, in fact involves 
large-amplitude periodic coefficients even for the low-frequency response of 
the hovering rotor. There is a 2/rev variation of the coefficients due to the 
lack of axisymmetry of the two-bladed rotor. For v = 1 (no flap hub spring), 
the effect is mainly on the helicopter pitch and roll damping and cross- 
coupling. 

. Third, the quasistatic approximation as implemented here, when applied 
to the two-bladed rotor, does not give the low-frequency response as it does 
for N > 3. The source of the difficulty is the fact that the 3i equation 
of motion is in the rotating frame still, so the 3i response to low-frequency 
inputs from the nonrotating frame is not at low-frequency also, but rather at 
1/rev. ' & £ ■ " % 

The special characteristics of the two-bladed rotor dynamics pose a 
number of problems for the analysis of the aeroelastic behavior. Generally, 
it is necessary to use the Floquet analysis of the periodic coefficient equa- 
tions more often than for a rotor with three or more blades. In fact, it is 
not possible to use directly the constant coefficient approximation (sec- 
tion 3.2) for flight dynamics since that eliminates the coupling of the rotor 
and the shaft motion. The quasistatic rotor model is very useful for heli- 
copter flight dynamics investigations, for N - 2 as well as N> 3 > Some 
procedure other than that of section 6.2 is required, however, to obtain the 
quasistatic representation of the two-bladed rotor. The simplest procedure 
is to use an equivalent N> 3 model for the rotor. Then the quasistatic 
approximation gives the desired low-frequency , constant-coefficient response 
of the actual two-bladed rotor. For the teetering rotor helicopter, a three- 
bladed gimballed rotor is a good choice for the equivalent model. The funda- 
mental parameters of the rotor (y, a, etc.) must be maintained; hence the 
equivalent rotor will have a chord and mass distribution sealed by a factor 
2/N e quiv A frequent use of such an equivalent model would be to represent a 
two-bladed tail rotor. 

The validity of these approximate analyses of the two-bladed rotor — the 
constant coefficient approximation and the equivalent rotor representation — 
must always be verified for a particular application, of course. While some 
useful range of validity may always be expected, eventually the periodic 
coefficients or high-frequency dynamics become important enough to require a 
more rigorous analysis. 
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11. CONCLUDING REMARKS 
11.1 Applications of the Analysis 

The aeroelastic analysis developed here has been applied in a number of 
investigations of rotorcraft dynamics, both to check the basic features of the 
analysis and to obtain information about the dynamic behavior of specific ' 
rotors and aircraft. References 18 and 21 present some results of these 
investigations. In reference 18, results are given for a number of classical 
problems of shaft-fixed rotor dynamics. The flapping frequency response to 
pitch control inputs is presented, including an examination of the influence 
of the rotor inflow dynamics for articulated and hingeless rotors. A root 
locus of the flapping stability of an articulated rotor in forward flight is 
given, including the influence of the periodic aerodynamic coefficients at 
high advance ratio. Thirdly, reference 18 presents flutter and divergence 
stability boundaries for an articulated rotor in hover. The influence of the 
offset between the center of gravity and the aerodynamic center, of the first 
bending mode, and of forward flight on the flutter boundary is examined. 

The rotor and wind-tunnel support aeroelastic analysis has been applied 
to several configurations. A number of calculations have been made of the 
ground resonance stability of articulated rotors on a test module, strut, and 
balance frame system; reference 18 presents typical results, including a 
comparison with an elementary stability criterion. Reference 19 gives the 
aeroelastic stability calculations for gimballed and hingeless proprotors on 
a cantilever wing. The proprotor and cantilever wing model has also been 
used in an investigation of optimal control designs for gust alleviation 
(ref. 20). Finally, reference 21 presents the predicted dynamic stability for 
a tilting proprotor aircraft in a wind tunnel, including the airframe, strut, 
and balance dynamics. 

The rotorcraft in flight aeroelastic analysis has been used in refer- 
ence 18 to calculate the flight dynamics of four representative helicopters: 
a small articulated rotor helicopter, a large articulated rotor helicopter, 
a soft-inplane hingeless rotor helicopter, and a tandem rotor helicopter. The 
results include an examination of the influence on the flight dynamics of the 
quasistatic rotor model, the rotor lag motion and other degrees of freedom, 
the rotor inflow dynamics, and coupled lateral and longitudinal aircraft 
motion. Finally, reference 21 presents the predicted dynamic characteristics 
of a tilting proprotor aircraft in flight, including trim conditions, flight 
dynamics, gust response, aeroelastic stability, and wing response to control 
inputs . 


11.2 Future Development 

An aeroelastic analysis for a rotorcraft in flight or in a wind tunnel 
has been developed, in which the dynamic behavior is described by a set of 
linear differential equations. From these equations, the dynamic stability, 
flight dynamics, and aeroelastic response of the system may be calculated, and 
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they form the basis for more extensive investigations such as automatic 
control-system design. It is not possible to anticipate all features that 
will be required to model future rotor designs, so it must be expected that 
new applications will often require further development of the model, some- 
times by minor extensions and sometimes by major ones. Thus, in addition to 
its current use in investigations of rotor dynamics, the present analysis also 
provides the basis for the continuing development of models for rotorcraft 
aeroelastic behavior. 
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APPENDIX A. ROTOR INERTIAL AND AERODYNAMIC COEFFICIENTS 


Al. ROTOR INERTIAL COEFFICIENTS 


The inertial coefficients required for the rotor equations of motion 
(see section 2.2) are 
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and for rigid pitch (k = 0)r, 
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A2 . ROTOR AERODYNAMIC COEFFICIENTS 


The aerodynamic coefficients required for the rotor equations of motion 
(see section 2.3) are as follows. Recall that these coefficients are constant 
for axial flow, but are periodic functions of i|) for nonakial flow. The 
coefficients for blade bending are 
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The aerodynamic coefficients for the flap moment are 
M = I F r dr 
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The aerodynamic coefficients for the other hub forces and moments follow 
the pattern of the flap moment, with the following changes in the notation 
and integrands: 
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The radial force coefficients are 
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and F r ^ and F r are similarly defined. Finally, the aerodynamic coefficients 
for the blade pitch and torsion are 
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APPENDIX B. MATRICES OF ROTOR EQUATIONS OF MOTION 


Bl. INERTIAL MATRICES FOR ROTOR EQUATIONS 


The inertial matrices for the rotor equations of motion in the nonrotat- 
ing frame (see section 2.6.1) are given below. For clarity, the superscript * 
denoting the normalization of the inertial coefficients has been omitted. The 
inertial coefficients are defined in appendix Al. 
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B2. AERODYNAMIC MATRICES FOR ROTOR EQUATIONS IN AXIAL FLOW 


The aerodynamic matrices for the rotor equations of motion in axial flow 
(see section 2.6.2) are given below. The aerodynamic coefficients are defined 
in appendix A2. 
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£3. AERODYNAMIC MATRICES FOR ROTOR EQUATIONS IN NONAXIAL FLOW 

The aerodynamic matrices for the rotor equations of motion in nonaxial 
flow (see section 2.6.3) are given below. Note that each matrix is a summa- 
tion over all the blades, that is, m = 1, ..., N. The notation C = cos 
and S = sin is used in these matrices. The aerodynamic coefficients are 
defined in appendix A2. 
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56. INERTIAL MATRICES FOR TWO-BLADED ROTOR 


The inertial matrices for the two-bladed rotor equations of motion (see 
section 4.2) are given below. The notation C = cos i{> m and S = sin ijj m is 
used, where = \(; + mir. 




































































B7. AERODYNAMIC MATRICES FOR TWO-BLADED ROTOR 


The aerodynamic matrices for the two-bladed rotor equations of motion (see 
section 4.2) are given below. The notation C = cos i|> m and S = sin i|j m is 
used, where = ip + mir. 
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APPENDIX C. AIRCRAFT CONTROL TRANSFORMATION MATRICES 


The transformation matrices between the pilot controls and the individual 
rotor controls (T^^, see section 7.2) are given below for the single main- 
rotor and tail-rotor and the tandem main-rotor helicopter configurations. The 
K values are gain factors in the control system, and the Aip values are 
swashplate azimuth lead angles. The main rotor or the front rotor is assumed 
to be rotor 1 and the tail rotor or rear rotor is rotor 2. The parameter Q 
here takes the value +1 for counterclockwise rotation of the rotor and -1 for 
clockwise rotation. 



CFE single 
main 
rotor 







CFE tandem 
main 
rotors 





For the side-by-side or tilting prbpotor configurations,- the lateral 
symmetry of the aircraft allows the control transformation to be separated 
into symmetric and antisymmetric matrices as follows. For a tilting 
propotor aircraft^ th^ gain parameters (K) would generally be' functions of 
the pylon tilt angle a . 
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APPENDIX D. MATRICES OF AIRCRAFT EQUATIONS OF MOTION 


Dl. INERTIAL MATRICES FOR AIRCRAFT EQUATIONS 


The inertial matrices for the aircraft equations of motion are given 
below (see section 8.3). These matrices also include the gravitational forces 
and the structural forces for the elastic body modes. Superscript * denoting 
the normalization of the parameters has been omitted for clarity. 
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D2. AERODYNAMIC MATRICES FOR AIRCRAFT EQUATIONS 


The aerodynamic'' matrices ‘ for-' the aircraft' ‘equations pf motion (see 
tion 8.4) are given below. The aerodynamic coefficients are defined in 
appendix E. 
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APPENDIX E. AIRCRAFT AERODYNAMIC COEFFICIENTS 

The aerodynamic coefficients required for the symmetric and antisynmetric 
aircraft equations of motion (see section 8.4) are defined as follows. 

c ee ' - ™ Vq y{[*\ A + - “«da + v]wb 

+ [ x2k la + z2k du - m(k da + Vim} 

C x§ Oa Vq y P 2K DU ~ xK DA^WB + ^ zK DV ” XK DA^HT^ 

C ae ' - ™ Vq yt <zK LU - ltK LA>WB + (zK lo - *Vht3 
C q lc § oa ^y C -e x^ ZK DU ” XK DA^WB ” ®z^ Z ^LU " ^LA^HB 

+ *k (zK BU ‘ xK da ) ht + z k (zK i» - 
c ei ' - ™ %[ <zK du • =*10)10 ' Vhb 

+ C 1 - f ) <zK du ‘ *Wht] . 

C xx oa V + ^"DU^HT^ 

C zx oa V ^ K lu^WB + ^U^HT ^ ' 

C q k x “ ~ oa V £“ (e x K DU + 6 z K LU^WB + ^\ K DU + Vlu^hJ 

C 0 z oa Vq y K " xK LA^wb “ ^ K MA^WB 

+ f 1 " 9a)^ zK DA “ ^^^LA^HtJ 

C xz oa V [f K DA^WB + ( 1 3 a ^^A^Ht] 
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Subscripts WB, HT, and VT refer to the wing-body, horizontal tail, and 
vertical tail, respectively. The components of the aircraft angular velocity 
introduce the following factors: 

q = cos 4> 

H y FT 


r 

x 


-sin 0 


FT 



= sin <f> FT cos 0 pT 

= COS ()>__ cos 0 m 
y FT FT 


The quantities x and z are the location of the aerodynamic surface center 
of action, in body axes (F system), relative to the aircraft center of 
gravity. The parameter A is the rotor area (ir R 2 ) ; S w , c w , and are, 
respectively, the wing area, chord, and span; and £q are the inboard and 
outboard edges of the wing control surface; \p w and <j> w are the wing sweep and 
dihedral angles (assumed to be small), the horizontal tail and vertical tail 
spans are £ H and Zy, respectively. 


In airplane analyses, it is conventional to use coefficients based on 
the wing area S w . A rotorcraft usually does not have a wing and, generally, 
there is no good reference area for the airframe aerodynamics. Thus, here 
the aerodynamic force characteristics are used in the form L/q, which have 


dimensions of length-squared (q is the dynamic pressure; the moment character- 
istics, M/q, have dimensions length-cubed). This form is appropriate for the 
analysis of a specific vehicle, where the scaling with velocity, but not with 
size, is of concern. The aerodynamic characteristics required for the wing- 
body description are the lift, drag, and pitching moment (L/q, D/q, and H^/q) , 
and their derivatives with respect to angle of attack, flaperon deflection, 
and aileron deflection; and the side force, rolling moment, and yawing moment 
(Y/q, N x /q, and N z /q) and their derivatives with respect to sideslip angle. 


The rolling and yawing moment derivatives N„ and N 

6 3 

alone (no vertical tail contributions), and without the 


are for the wing-body 
sweep and dihedral 


terms already included in C.* and C... The ; vertical and horizontal tail 

<py wy 
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aerodynamic characteristics required are lift and drag and their derivatives 
with respect to angle of attack and control-surface deflection. To account 
for the velocity vector not being aligned with the x axis, the aerodynamic 
characteristics are required in the following combinations for the aerodynamic 
coefficients:. ...... - .. ; 

Sa * i[("f + f) c °s 2 +V + 2 f sl “ z *v + ■(-f + |) sin * v cos t v j 

Sa ‘ A [(f ' q) COs2 *V ' 2 q Si " 2 ' *V ' (f ’ Sl " +V C ° S *v] 

Su* a [ 2 * | c°> 2 *T ■ (‘f ' q) sin! *V -(-f- f) Bin *V 1 , 

V ■ A [ 2 f C0s2 *V + (f + q) Sl " 2 *V ‘(f + q) sln »V C0S +v] 

Sa - S (t cos *V + 2 f *v) 

' ‘ ' r ' , * 

Su - S ( 2 | *v - f sln *v) 

"it * A (■# COS *V + T Sl ” *v) 

t 

» /D L \ .*■ . ‘"i 

St ‘ A VT C ° 8 *V - T Si ” *v) 

he ~ A (v cos +v + -f sin *'v) 

Sc ■ a (t cos *v ' ~i sln *v) . : 

Sw - X [(t + q) ~ »V + (t - q) Sl ” 4 ’ 

Sw • A [(t ' q) cos *V • (t + ?) aln *v] \ :■ 

& V * ' 

218 



Here 


is the angle of attack of the reference axis system, so 


<Py “ tan" 1 V z /V x for the wing-body and horizontal tail and ^ = tan” 1 Vy/V x 
for the vertical tail. 

The wing-induced velocity at the horizontal tail is accounted for by the 
derivative 3e/3a. The following expression is used (from ref. 16): 


0.45 (L 


_3e 

3a 


( £ 2/ s )0.735 ( o /c ) 0 . 25 
w w HT w 


where is the tail length and is the wing lift-curve slope. The 


a 


side-wash velocities at the vertical tail are given by 3o/33, (V/z VT )3a/3p, 
and (V/xyj) 3c/3r for sideslip, rolling, and yawing motion, respectively. 
Typical values are 3o/33 and (V/xyj)3cr/3r near zero and (V/zy^)3a/3p 
approximately 1 (ref. 15). 

Finally, the required integrals of the wing bending and torsion motion 
for the elastic degrees of freedom are: 
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The y^ terms are absent for the symmetric modes; the integrals e^, 
e^, and e Q ^ are required only for antisymmetric modes. 
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(a) Shaft motion. 


T,-Q 

A 


H, M, 


(b) 


Hub reactions. 


Figure 6.- Notation and sign conventions for linear and angular shaft motion 
(displacements in an inertial frame) and forces and moments acting on rotor 
hub (in nonrotating frame) . 
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Rotor 

shaft 


(a) Side view. 


Figure 7.- Schematic of the rotor hub and root geometry showing the pitch bearing radial offset (r^ ) , 

torque offset (x_, A ), precone (5_,. ), droop (6 ), sweep (S_. ), feathering axis droop (6_. ), and 

rA rAj rA 2 rA3 rA4 

feathering axis sweep (6 )• Only a single undistorted blade is shown without the gimbal under- 

1*A5 

sling (z ); gimbal is dropped from the model for articulated and hingeless rotors. 

FA 
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(b) Top view. 
Figure 7.- Concluded. 



Elastic axis (EA) 


Figure 8 


\ 

Locus of section \ 
center of gravity (CG)\ 


Geometry of undeformed blade. 









Figure 11.- Notation and sigr conventions for rotor velocity and orientation. 
(V and a Hp ), induced velocity (v^ , and aerodynamic gust velocity compo- 
nents (u„, v„, w ). 



Pitch 

bearing 
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Pitch 

/^y bearing 


236 



(a) Single main-rotor and tail-rotor helicopter. 

Figure 13.- Typical rotorcraft configurations considered by the aeroelastic 

analysis. , 





(d) Tilting proprotor aircraft. 
Figure 13.- Concluded. 


I'O 

OJ 

VO 


240 


Rotor# I 


Rotor #2 




Rotor #1 
transmission 


Rotor # I 

shaft 

spring 




K i 

■WWW 


Inter -connect 
-spring 


k m- 


cb fl 2 


Rotor #2 

shaft 

spring 



Rotor #2 
transmission 


Engine 

shaft 

spring 


Engine 

inertia 


(a) Single main-rotor and tail-rotor or tandem main-rotor 
helicopter configurations. 


Figure 14.- Schematics of rotor craft transmission and engine dynamics models. 
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(b) Side-by-side or tilting proprotor aircraft configurations. 
■- it-- Figure 14.- Concluded. 
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